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ABSTRACT 


This document contains a detailed description of the computer code for the 
trajectory processor (#TRAJ) of the HP-9825A High Fidelity Relative Motion 
Program (HFRMP). The #TRAJ processor is a 12-degrees-of-freedom trajectory 
integrator (6 degrees of freedom for each of two vehicles) which can be used 
to generate digital and graphical data describing the relative motion of the 
Space Shuttle Orbiter and a free-flying cylindrical payload. The description 
of the processor includes a listing of the code, coding standards and conven- 
tions, detailed flow charts, and discussions of the computational logic. 
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FOREWORD 


This document contains a detailed description of the trajectory processor 
(#TRAJ) of the High Fidelity Relative Motion Program (HFRMP) Version 03T ( 1 229/ 
10N0V79). The HFRMP is coded in Hewlett-Packard Language (HPL) for the 
HP-9825A Desktop Calculator with memory option 002 (23,228 bytes of read/ 
write memory) and the following peripheral devices and read-only memory (ROM) 
modules : 


Required Peripheral Devices 

HP-9885M Flexible Disk Drive 

HP-9872A Plotter 

HP-9866B Thermal Line Printer 

or 

HP-9871A Character Impact Printer 

Required ROM Modules 

HP-9821 OA String-Advanced Programming 

HP-98211 A Matrix Programming 

HP-98217A 9885M Flexible Disk Drive 

HP-9821 6A 9872A Plotter-General I/O-Extended I/O 

or 

HP-9821 5A 9872A Plotter-General I/O 


Applicable HPL programming information is contained, along with operating in- 
structions for the calculator and its peripheral devices, in the following 
HP-9825A manuals published by the Hewlett-Packard Calculator Products Divi- 
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sion: 


HP Part No. 

Title of Manual 

09825-90000 

Operating and Programming 

09825-90020 

String Variable Programming 

09825-90021 

Advanced Programming 

09825-90022 

Matrix Programming 

09825-90024 

General I/O Programming 

09825-90026 

9872A Plotter Programming 

09871-90000 

9871 A Printer .. .Operating Manual 

09885-90000 

Disk Programming 


A basic understanding of the contents of the cited manuals is required for a 
full comprehension of this program document. 


1 . INTRODUCTION 


The HP-9825A High Fidelity Relative Motion Program (HFRMP) is a disk- 
resident flight design software system specifically intended for the analysis 
of Space Transportation System (STS) proximity operations. The mainline logic 

4 * 

flow is illustrated in Figure 1. The code is divided into several primary 
modules that are designated "processors" and identified as such (according to 
HFRMP convention) using the pound-sign (#) symbol as the first character in 
their identification strings (i.e., names). Some processors, being too large 
to fit in the available read/write memory area, are broken down into subsidiary 
links which are identified by the use of the percent (") symbol as the first 
character in their identification strings. Subsidiary links overlay each other 
in a shared memory area, as illustrated in Figure 2. 

Processors (such as #TRAJ) are called into read/write memory from the 
program disk by use of the HPL get command, which erases all variables (numer- 
ical data) that may have been stored in memory. Information transfer between 
processors is accomplished by means of disk-resident data files that are 
stored and retrieved, as necessary, by appropriate commands coded into the 
processors. This includes user-supplied input data, which can be manipulated 
(edited/saved/recalled) by the user with the aid of interactive input proces- 
sors (#DBED and #FPED in Figure 1). 

From the preceding remarks it can be seen that every HFRMP processor rep- 
resents , in a limited sense, a stand-alone program. Therefore it is possible, 


i.e., that which is invoked by a regular production run of the program. Not 
shown are a number of special processors that are used for program mainte- 
nance. 
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TEST I : BODY- FIXED RELATIVE MOTION PLOTS REQUESTED? 
TEST i: MORE FLIGHT SEGMENTS ? 

TEST 5 : DID X PAYLOAO SRM BURN OCCUR ? 

TEST M PLUME IMPINGEMENT DAMAGE ANALYSIS REQUESTED? 
tests: start a new run? 

Figure 1. HFRMP Mainline Logic Flow 
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for instance, to interrupt the normal HFRMP execution sequence r'ter input 
processing has been completed and then later to resume execution by commanding 
get "ITRAJ" from the calculator keyboard and pressing the CffuH) button. By the 
same token, it is not necessary to execute the HFRMP particle impact damage 
integrator ( #P I DI ) immediately after executing #TRAJ. The vehicle ephemeris 
data that #PIDI requires is stored on the program disk by #TRAJ, where it is 
preserved indefinitely until written over by a subsequent execution of #TRAJ. 
Therefore, a user can execute the trajectory processor, shut off the calculator 
and remove his program disk, and then initiate a particle impact damage analy- 
sis of the subject trajectory at any later date simply by inserting his disk, 
commanding get M #PIDI " from the keyboard, and pressing the (ROij) button. 

The remainder of this document will be concerned only with the #TRAJ pro- 
cessor, which contains the HFRMP dynamic models, the trajectory/attitude inte- 
gration logic, and the basic output computations. Detailed program documen- 
tation for the HFRMP input processors and the #PIDI output processor are to be 
published under separate cover. (Reference 1 contains operating instructions 
for the HFRMP, with emphasis on input operations. Reference 2 contains a def- 
inition of the ; • i c equations used in the #PI DI processor.) 
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1.1 GENERAL DESCRIPTION OF #TRAJ 


The fTRAJ processor is a 12-degrees-of-freedom trajectory integrator (6 
degrees of freedom for each of two vehicles) which generates digital and graph 
ical data to describe the relative motion of the Space Shuttle Orbiter and a 
free-flying payload. These data are obtained by differencing the geocentric 
states of the individual vehicles, computed to a numerical precision of 12 
decimal digits, with respect to an oblate earth whose gravitational model 
includes the second harmonic coefficient (J 2 ). The state of the two-vehicle 
system > computed as a function of time by means of a fourth-order Runge- 
Kutta numerical integration scheme which uses unit quaternions (versors) to 
define both the rotational and the translational states of each vehicle. The 
derivation of the versor-based translational state vector is contained in 
Reference 3, which is reproduced for the most part (with some changes in 
coordinate system definition) in Appendix A and Sections 1.2 and 1.3 of this 
report. 

The payload is modeled geometrically as a cylinder whose length and diam- 
eter are specified by program input. The Orbiter and the payload are treated 
as rigid bodies whose individual mass properties (qcoss weight, moments and 
products of inertia, and center of gravity location) are specified by program 
input and are assumed to remain constant f during the HFRMP run. Gravity 
gradient torque is included in the rotational equu ' ; ons of motion for both 
vehicles. At the user's option, aerodynamic torque and drag can also be 


f 

An exception is made when integrating the trajectory of an upper stage during 
a major translational maneuver such as a solid rocket motor burn. In such 
a case, the payload gross weight is decremented during the burn to reflect 
the consumption of propellant. 
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included for either or both vehicles. Aerodynamic forces normal to the rela- 
tive wind vector are ignored. 

The atmosphere is assumed to rotate with the earth, and is modeled as a 
function of geodetic altitude by a curve fit of the 1962 Standard density pro- 
file (Reference 4). The curve fit is valid down to a minimum altitude of 
approximately 400,000 feet. The density profile can be modified by a program 
input factor (K p ) to account for the major effects of solar activity. Aerody- 
namic drag and moment coefficients for the Orbiter are computed as curve-fit 
functions of its attitude with respect to the relative wind vector (References 
5-7). 

Aerodynamically, the payload is modeled as a flat plate whose size and 
shape are determined by projecting its cylindrical outline onto a plane normal 
to the relative wind vector. The payload dr&g coefficient is assumed to be 
2.0, based on its projected frontal area. Aerodynamic effects on the Orbiter 
and the payload can be modified (or cancelled entirely) by means of input fac- 
tors (Ka and ka) which are applied uniformly to all aerodynamic forces and 
torques that are computed for the specified vehicle. The possible effects of 
aerodynamic shadowing (of one vehicle by the other) are not accounted for in 
the internal calculations of #TRAJ. 

Several options are provided for defining the initial state of the two- 
vehicle system. The translational state of the Orbiter can be described 
either in terms of osculating orbit elements referenced to the Mean of 1950.0 
(M50) geocentric equatorial frame, or in terms of invariant orbit elements 
(Reference 8) measured in the Mean of Launch Date (MLD) equatorial frame. The 
initial attitude of the Orbiter is defined by pitch, yaw, and roll angles 
(taken in that order) referenced to the rotating Shuttle-centered local- 
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vertical (SLV) coordinate system. The Orbiter's angular velocity, measured 
in terms of rate components about its body axes, can be defined with respect 
to either the M50 (inertial) frame or the SLV (rotating) frame. 

The initial translational state of the payload is defined by rectangular 
position and velocity components which are measured rel ive to the Orbiter's 
center of gravity (CG). At the user's option, these components can be mea- 
sured in the SLV coordinate system, or in the Shnttle body (SBY) coordinates 
system. The payload's initial pitch, yaw, and roll angles can be referenced 
either to the payload-centered local -vertical (PLV) system, or to the SBY 
system. The payload's angular rate components about its body axes can be 
defined relative to the M50, the PLV, or the SBY frame. 

The initial state of the system can be advanced through up to 40 flight 

profile segments, each of arbitrary length, which are defined by the user in 

a prior execution of the #FPED processor. At the beginning of any segment 

the user may command the application of an impulsive (i.e., instantaneous) 

★ 

increment to the angular rate of either or both vehicles. In this regard, 
the user may specify a particular rate increment (INCR), a desired rate with 
respect to inertial space (IR), or a desired rate with respect to the local- 
vertical frame of the vehicle in question (LVR). In all cases, the components 
of the desired rate or rate increment are measured about the body axes of the 


►f* 

All input data, including the flight profile definition, are saved in disk 
files whence they can be recalled (and edited, if necessary) for use in 
subsequent runs. 

ic 

Commanded angular velocity impulses, and the linear velocity impulses which 
they may induce as a result of RCS translational cross coupling, are the only 
types of state variable discontinuity that are permitted by #TRAJ. These are 
allowed only at the beginning of a flight profile segment. 
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vehicle in question. 

After applying the specified angular velocity increment (if any) at the 
beginning of the flight segment, ^TRAJ then (for each vehicle independently 
of the other) 

1. allows the attitude to drift (P) under the influence of inertia and 
natural torques, or 

2. performs inertial rate hold (IRII) control (i.e., maintains a constant 
angular velocity relative to inertial space), or 

3. performs local-vertical rate hold (LVRH) control (i.e., maintains a 
constant angular velocity relative to the rotating local-vertical 
frame of the vehicle being controlled) 

for the duration of the segment, depending on the user's specifications for 
that segment. 

When the IRH or the LVRII attitude-maintenance option is specified for the 
Orbiter, a simplified RCS/DAP model (Reference 9) is used to compute average 
values for the propellant consumption rates and translational cross-coupling 
accelerations that result from the intermittent thruster firings which are 
required to apply the necessary control torques. The model takes into account 
the mass properties of the Orbiter, the electrical width (an integer multiple 
of the DAP cycle time) and the effective width (the duration of steady-state 
acceleration) of the RCS thruster pulses, and the width of the attitude dead- 
band about each of the Orbiter's body axes. Deadbands can be changed from 
segment to segment in the flight profile, as can the selection of primary or 
vernier thrusters and the mode of cross -coupling compensation. Translational 
cross-coupling accelerations are integrated along with those produced by 
gravity, aerodynamics, and commended translational thrust. They are reflected 
in the output data by the flight path deviations they produce. Propellant 
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consumption rates are also integrated (but not subtracted from the Orbiter 
gross weight), and the accumulated expenditures are tabulated, along with other 
data, at user-specified time intervals. Propellant consumption is broken down 
according to source (forward, aft left, or aft right tank) and function (trans- 
lational or rotational control). 

When the IRH or the LVRH attitude-maintenance option is specified for the 
payload, the magnitudes of the necessary control torques are integrated and the 
accumulated rotational impulse (measured in pound-foot-seconds) in the positive 
and negative direction about each body axis is printed along with the other 
digital output data. Since no specific method of implementation is modeled, 
it is not possible to compute propellant consumption rates or cross-coupling 
effects that may result from payload attitude control. 

Translational thrust acceleration of either or both vehicles can be 
commanded at the beginning of any flight profile segment. Payload transla- 
tional thrust is always applied in the direction of the payload's +X body axis 
and is assumed to be directed through the CG. Once initiated, payload thrust 
acceleration continues until all of the rocket motor propellant is consumed, 
as determined by a table of flow rates versus burn time. 

Translational acceleration of the Orbiter is initiated by commanding 
ignition of either or both of the OMS engines (L, R, or L+R) and/or by firing 
primary RCS thrusters to produce thrust nominally in the positive or negative 
directions of the Orbiter body axes (+X, -X, +Y, -Y, +Z, or -Z). Once initi- 
ated, Orbiter translational acceleration is applied continuously at the 
nominal steady-state level, throughout the duration of the flight profile 
segment. 
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1.2 STATE VECTOR DEFINITIONS 


There are many numerical methods which can be used to integrate the dif- 
ferential equations that govern the motion of a rigid body such as the Orbiter 
or a free-flying payload. A common requirement of such methods is the defi- 
nition of a state vector 1 along with a set of equations which, given the val- 
ues of the state variables at any particular time, permit the calculation of 
the first derivative of every state variable with respect to time. 

A fourth-order Runge-Kutta numerical integration technique is used in the 
#TRAJ processor. However, the state vector definitions and the derivative 
calculations (which comprise the bulk of the computations involved in state 
vector propagation) are not influenced directly by the integration logic 
except in the sense that the integrator determines how frequently the deriv- 
atives have to be evaluated. In #TRAJ, the integration logic and the deriva- 
tive calculations are isolated from each other in separate subroutines, thus 
making it easy to change the numerical integration technique if that should 
become desirable. 

The purpose of this section of the report is to describe the state vari- 
able that have been adopted for use in #TRAJ. The associated derivative 
calculations, described in terms of external forces and torques that act on 
the body under consideration, will be covered in Section 1.3. The detailed 
calculation of these forces and torques will be described in Section 2.3, and 
the Runge-Kutta integration logic will be described in Section 7.2. 


4 . 

A set of variables that define the position and linear velocity of the body's 
CG (in the case of translational motion) with respect to a chosen system of 
reference, and/or the body's attitude and angular velocity (in the case of 
rotational motion about the CG). 
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The #TRAJ processor makes use of versors (unit quaternions) as internal 
state variables to describe the translational as well as the rotational 
motions of the Orbiter and the payload. (The HFRMP user never sees the state 
vectors in quaternion form; input and output data are defined in terms of 
Euler angles.) The computational advantages of using versors to describe 
rotational motion are well known. Certain advantages also accrue from their 
use to describe the translational motion of a satellite, as explained in 
Reference 3. Versors are also used extensively in #TRAJ for coordinate trans- 
formations in general, and a working knowledge of their characteristics and 
rules of manipulation is necessary for a good understanding of the remainder 
of this report. A careful study of Appendix A is recommended at this point 
for the reader who is not accustomed to the use of quaternions. Even for 
those who are familiar with quaternion applications, it probably will be help- 
ful to review the conventions and the system of notation defined in Appendix 
A. All readers should examine the coordinate system definitions in Appendix B. 

1*2.1 Rotational State Vector 

The rotational state of the Shuttle at any given time t is defined by the 

___ T D __ 

variables [qjg,ug ], where q IB is a versor that defines the orientation of the 
Shuttle body-fixed (SBY) coordinate system B with respect to the geocentric 
mean-of-launch-date (MLD) equatorial inertial coordineite system I. The vector 

T D 

symbol tog represents the angular velocity of the Orbiter. The double super- 
script IB identifies the angular velocity as that of system B with respect to 
system I, and the subscript B indicates that the vector is resolved into 
components along the axes of the B system. The rotational state of the payload 
is similarly defined by the variable? [q^,^], whose definitions are the 
same in all respects except that the lower-case symbol b is used to identify 
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the body-fixed coordinate system of the payload. 


1.2.2 Translational State Vector 

The translational state of the Orbiter is defined by the variables 
[qj G ,R,R,fi^]. The symbol q IG represents a versor that defines the orientation 
of the Shuttle local -vertical (SLV) coordinate system G with respect to the 
inertial system I. The symbol R represt nts the distance from the center of 
the earth to the Shuttle CG, and R stands for its derivative dR/dt. The 
symbol represents the scalar magnitude of the instantaneous geocentric 
angular velocity of the Shuttle CG. The subscript K in this case does not 
designate a coordinate system; it merely identifies the angular velocity in 
question as that which is associated with Kepler's law of equal areas. 

The translational state of the payload is similarly composed of the vari- 
ables [qjg S r,r,co^] , which are defined in the same manner as those described in 
the preceding paragraph except for the use of lower-case symbols to repre- 
sent the scalar quantities and the coordinate-system designator (g) of the 
payload local-vertical (PLV) coordinate system. 
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1.3 BASIC EQUATIONS OF MOTION 


The purpose of this section is to define the equations for calculating 
the derivatives of the #TRAJ internal state variables, assuming the forces and 
torques acting on the bodies are known. Equations for calculating the forces 
and torques will be given in Section 2.3. 

The only equations that will be defined explicitly in this section are 
those that apply to the motion of the Shuttle. The payload equations are 
identical in all important respects, and can be obtained simply by substituting 
the lower-case symbols b, g, r, r, and oi^ for their upper-case counterparts 
B, G, R, R, and s^. 

1.3.1 Rotational Motion 

What is required here are the first derivatives of qTjg and with re- 
spect to time. The derivative of the orientation versor is equal to one-half 

T D 

of its quaternion product with Ug , i.e., 
q IB = ( q IB 0 W B V 2- 

TD 

The derivative of Wg is given by the equation 



where the vector Ug represents the external torque acting on the Shuttle. 

The symbol [I]g represents the Shuttle's inertia tensor (a 3x3 matrix whose 
elements consist of the moments and negative products of inertia, referenced 
to the B system), and [I]^ 1 represents its inverse. For the purpose of comput- 
ing the product of a matrix and a vector, the vector is treated as a column 
matrix, e.g. , 
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IB 

“B1 


-IB IB 
“B ~ “B2 



1.3.2 Translational Motion 

We seek now to define the first derivatives of the translational state 
variables q^ G , A, and with respect to time. (The derivative of the state 
variable R is no problem, because it is equal to R). By definition, the 
Shuttle's local -vertical coordinate system G rotates in such a manner that 


R G = - k R (3) 

and 

(RxV) q = - j|RxV| = - j R 2 a K (4) 

at every instant of time, where R and V represent the geocentric position and 

inertial velocity vectors of the Shuttle CG. 

An expression for the inertial velocity can be obtained by differentiat- 
ing Equation (3) with respect to time. This yields 

• T £ A 

V G = - k R - w G x k R, (5) 


where the vector 




+ j to 


IG 

G2 


+ k to 


IG 

G3 


( 6 ) 


(whose components have yet to be evaluated) represents the total angular 
velocity of coordinate system G with respect to the inertial system I. We 
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now differentiate Equation (5) with respect to time and obtain the expression 


A fi = - k R - 2 J q G x k R 


w g G x k R - o)g G xLj 


(“g G * k R ) 


for the linear acceleration of the Shuttle CG, where the vector 


7 -IG : *IG , n . -IG . r -IG 


represents the angular acceleration of coordinate system G. 

T g _^T g 

We now proceed to evaluate the components of wg and nig . First, we 
substitute the expression for wg from Equation ( 6 ) into (5), and expand the 
vector cross product. This results in the equation 

Vg = - i R cog G + j R ujg G - k R. (9) 

Using the expression for Rg and Vg from Equations (3) and (9), we form the 
cross product 

(RxV)g = i R 2 o)g G + j R 2 o)J G . (10) 

Equating the vector components of Equations (4) and (10), we obtain 

4? = 0 OD 


“G2 = ~ n K 


which lead to 


a>AT = 0 
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and 


•IG • 

“G2 " ‘ U K 


04 ) 


Substituting (11), (12), (13), and (14) into Equations (6), (8), and (9), we 
obtain 


-IG _ 5 , . j; IG 

Wg ~ “ 3 K tUf-.*} , 


G3 1 


(15) 


-iG _ '■ % . n * ig 

W G ' " J S “K + k ^63 ’ 


(16) 


and 


Vq = i R - k R. 


(17) 


r p • t p 

Using the expressions for tog and uig from Equations (15) and (16), the 
cross products opearing in Equation (7) can now be expanded to obtain 


Ag = i (2 R p. K + R si K ) + j R oigg + k (R fl* - R) 


(18) 


From Newton's second law, we have 


Fg = M Ag, 


(19) 


where M represents the mass of the Shuttle and the vector 


F G = 1 F G1 + J F G2 + k F G3 


( 20 ) 


represents the external force acting on it. Combining Equations (18), (19), 
and (20), and equating vector components, we obtain the equations 


R = R - Fg 3 /M 


( 21 ) 
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= (Fqi/M ■ 2 R 


( 22 ) 


and 

“G3 = F G2 /(M R fi K^‘ (23) 

Equations (21) and (22) define the derivatives of the translational state 

• =^= 

variables R and The derivative of the remaining state variable q IG is 

given by 

q j q = (Rjq 0 W G (24) 

where the angular velocity vector 

wg G = - j fi K + k F G2 /(M R n K ) (25) 

is obtained by substituting Equation (23) into (15). 

Equations (21), (22), (24), and (25) are valid for all states of motion 
except those where |RxV| = 0, in which case the geocentric flight path would 
be vertical, and the orientation of the X and Y axes in the local -vertical 
coordinate system would be indeterminate. 
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1.4 CONVENTIONS AND CODING STANDARDS 


HP-9825 programs are coded in the Hewlett-Packard Language (HPL), which 
is functionally similar to BASIC but unique in terms of its syntax. To 
facilitate correlation with the #TRAJ code (Appendix C), some of the syntacti- 
cal features of HPL have been carried over into the flow charts (Sections 2-7). 

For instance, consider a computation which involves forming the product 
of two simple variables A and B, where the resulting value is to be assigned 
to a third variable C. In FORTRAN code this operation would be represented by 
writing C = A * B , and in most flow charts by writing C = AB . In HPL code, 
and in the flow charts that follow, the same operation would be described by 
writing AB -» C . (Perhaps the most succinct verbalization of the HPL value- 
assignment symbol is obtained by substituting the phrase "goes to" for the 
right-running arrow "V.) 

Another feature of HPL code and the flow charts, which may be confusing 
to FORTRAN programmers, is the use of square brackets to enclose the index or 
indices of a dimensioned variable. The Jth element of the HPL dimensioned 
variable A[*] is designated by writing A[J], whereas in FORTRAN it would be 
written A(J). In HPL code, the square brackets are always used to enclose 
indexing or dimensioning information about an array (dimensioned) variable, 
and for no other purpose . This is necessary in HPL to distinguish, for 
instance, the simple (scalar) variable P from the array variable P[*] whose 
dimensions might have been declared by a statement dim P[5,10] , which is 
analagous to the FORTRAN statement DIMENSION P(5,10) . In a FORTRAN program 
where such a dimension statement exists, a reference to the variable P is an 
implicit reference to P(1,l), i.e., the first element in the P array. However, 
this is not the case in HPL, where the simple variable P has no connection 
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whatever with the array variable P[*] or any of its elements. For instance, 
the sequence of HPL statements 2 •> P; 3 •» PCI,!]; P P[l,l] •> X would cause the 
number 6 to be assigned to the variable X. 

1.4.1 General Arrangement of Program Files 

The first line of code (line 0) in every HFRMP program file (processor or 
processor link) contains, in the form of a label, the name under which the file 
is stored on the program disk, followed by the date and time of its most recent 
revision. The executable part of this line contains the statement gto "RUN" , 
where "RUN" is the label of the line where the main logic flow begins. The 
next two lines (1 and 2) contain statements that facilitate an automatic list- 
ing of the complete program, and which are never executed in a normal produc- 
tion run. 

All subprograms are located between line 2 and the "RUN" label, and they 
are arranged in the following general order: 

1. Function subprograms 

2. Utility subroutines 

3. Special subroutines. 

The rationale for this order is to make the program listings as insensitive as 
possible to corrections or revisions of the logic. Function subprograms and 
utility subroutines are revised much less frequently than the main (driver) 
logic in a given link. Therefore the code in the top part of each link tends 
to remain comparatively stable. 

To the greatest extent that is practicable, the code in each subroutine 
and in the driver logic of each individual link is arranged so that execution 
control flows from top to bottom, i.e., from smaller line numbers to greater 
ones. 
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1.4.2 Branching Conventions 


For ease of program maintenance, branching to absolute line numbers is 
avoided. Relative addressing (e.g., gto + 4) is used whenever practicable; 
otherwise, the general rule is to branch to a symbolic label (e.g., gto “RUN" ). 
The jmp statement is used only in the rare instances where it is necessary to 
compute a relative address (line number) at execution time. 

To avoid confusion regarding the flow of execution control (and also 
generally to permit the transfer of date through an argument list), the 
HPL sdL statement is used for branching to a subroutine in lieu of the gsb 
statement. 

1.4.3 Program Labels 

No lower-case alphabetic symbols are used in program labels (character 
strings used to identify branching targets, including the names of subpro- 
grams). The purpose of this convention is to provide a contrast between such 
labels and the commands and standard functions of HPL, all of which are spelled 
with lower-case alphabetic symbols. 

1.4.4 Input/Output Mnemonics 

Although numerical option codes are sometimes used internally by #TRAJ, 
it has been adopted as a policy that HFRMP users should not have to memorize 
(or look up) obscure numerical codes in order specify logical options in the 
input data files. Rather, such options are specified by the user in the form 
of mnemonic character strings which, if necessary., are converted to numerical 
codes internally. 

Mnemonic character strings are also used to identify the HFRMP digital 
output data. With regard to all input/output mnemonics, a determined effort 
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has been put forth to make them as meaningful and as consistent as possible. 


1.4.5 Data-Reglster Protocol 

Because of the limitations of the HPL syntax. It has been found necessary 
to adopt a comparatively rigid protocol to govern the use of named variables 
and of the numbered r-registers. Appendix D contains storage allocation infor 
mation for the r-registers. 

1.4. 5.1 Simple Variables 

The simple (scalar) variable names, which are limited in HPL to single 
upper-case alphabetic characters (A thru Z), are used primarily as loop 
counters and for the temporary storage of intermediate computational results. 
Eight of the 26 simple variables are designated as volatile registers, which 
means that they are least rigidly controlled, and that their contents are most 
susceptible to frequent change. These are the registers identified by the 
characters and W,X,Y,Z. 

1.4. 5. 2 r-Registers 

Because of the limited number of usable names for variables, and also 
because array names cannot be passed through the argument lists of HPL sub- 
programs, most of the data that would normally be assigned unique names or 
stored in individual arrays (e.g., as in a FORTRAN program) are stored instead 
in the numbered r-registers. The correlation between the r-numbers and the 
logical symbols, shown in Appendix D, is extremely critical to the understand- 
ing of the #TRAJ code. 

When there is a need for a utility subprogram to perform a standard compu- 
tation involving one or more logical arrays (such as vectors and quaternions), 
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the address of each logical array (i.e. t the r-nuinber of its first element) is 
passed through the argument list (in lieu of an array name, which is not per- 
mitted in HPL). For instance, suppose the components of a vector Rq reside 
in the registers rl through r3, the components of another vector Vq reside in 
r4 through r6, and it is desired to compute the vector H^ = Rq x 7g and store 
it in r7 through r9. Symbolically, in the flow charts, this would be written 
Rq x 7 G -► Rg, or possibly as c£L 'CRSP' (Rg» v g»H G ). In the HPL code, it would 
appear as c jj 'CRSP 1 (1 ,4,7) , where 'CRSP' is the name of the vector cross- 
product routine, and the numbers 1,4,7 are the addresses of the vectors in- 
volved in the operation. 

The first 19 r-registers (rO through rl8) are similar to the simple-vari- 
able registers in that they are used mainly to store the intermediate results 
of array computations. The first ten r-registers (rO through r9) are analo- 
gous to the simple variables H ,1 ,J ,K and W,X,Y,Z in that they are designated 
as volatile, i.e., most frequently used. Preferential use of the lower-num- 
bered registers tends to minimize code-storage requirements , simply because 
they entail writing out fewer digits (each of which requires one byte of code 
storage) to identify the registers. That is to say, the statement c U, 'CRSP ' 
(101 ,104,107) requires 6 more bytes of code storage space than the statement 
c Li 'CRSP' (1,4,7) . 

1.4. 5. 3 Restrictions on Utility Subprograms 

During the execution of an HPL subprogram, the HP-9825 operating system 
allocates temporary storage (only for the duration of subprogram execution) to 
numbered p-variables. These variables are numbered sequentially beginning 
with pO, into which the operating system loads the number of parameters that 
appear in the argument list of the calling routine's reference to the sub- 


program. If there should be three arguments, then at the beginning of sub- 
program execution pO would contain the number 3. The registers pi, p2, and p3 
would contain the values assigned to the arguments of the subprogram by the 
calling routine. The subprogram can make use of as many higher-numbered p-reg- 
isters (p4, p5, etc.) as it may need for temporary storage of intermediate 
computational results. All of the p-registers are de-al located (in effect, 
erased) when execution control is returned to the calling routine. 

The operating system also permits the subprogram to reference (get values 
from or store values into) any of the registers accessible to the main program 
(this includes the simple variables, array variables, string variables, and 
r-registers) . However, to prevent inadvertent modification of the contents of 
registers that the calling routine may be using, certain conventions have been 
adopted that limit the access of utility subprograms to global variables (i.e., 
any other than the p-variables) . 

In general, utility (general-purpose) subprograms are not allowed to make 
literal reference to a"v r-register. They are allowed to i .ke logical refer- 
ence to such registers by means of accesses that may be passed to them 
through their argument lists by the calling routine. For example, the charac- 
ter strings rO or r25 are not permitted in the code of a utility subprogram. 
However, the character strings rpl and r(p3+3 ) are permitted, where pi and p3 
represent addresses (in this case, the first and third numbers in the argument 
list) of logical arrays that are passed to the subprogram by the calling 
routine. 

Function (as distinguished from subroutine) subprograms are not allowed to 
make reference to any simple variable. They are required to use p-variables 
for any temporary storage they might need. The same restriction applies to 
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any utility subroutine whose logical argument list contains only scalars and 
vectors . 

Utility subroutines that perform quaternion and matrix operations are 
permitted to use the volatile simple variables (H,I,J,K and W,X,Y,Z), but no 
others. Special subroutines such as ' DERVIS ' (Section 2.3), whose calculations 
generally are more complicated than those of utility subroutines, have no 
such general restrictions on their access to storage registers. 


2. BASE LINK (#TRAJ) 


The "#TRAJ" program file is the base link of the #TRAJ processor, contain- 
ing all of the general-purpose subprograms, i.e., all that are used by two or 
more subsidiary links. 

2.1 FUNCTION SUBPROGRAMS 

2.1.1 mi 1 

The function subprogram 'ANGT converts an input angle to its equivalent 
value in the range of 0 to 2tt. 

2. 1.1.1 Arguments 

pi = Input angle, measured in radians. 

2. 1.1. 2 Fxample of Usage 

The instruction 1 ANG1 1 (-3n) A in the calling routine would cause the 
value of +7T to be assigned to the variable A. 

2. 1.1. 3 Computations 

The fractional part of pi /2-rr is multiplied by 2ir. If the result is nega- 
tive, 2 tt is added. A flow chart is shown in Figure 3. 
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2.1.2 'ANG2* 


The function subprogram * ANG2 ' converts an input angle to its equivalent 
value in the range of -n to +ir. 

2. 1.2.1 Arguments 

pi = Input angle, measured in radians. 

2. 1.2. 2 Example of Usage 

The instruction 1 ANG2 ' (-3tr) -v A in the calling routine would cause the 
value of -7T to be assigned to the variable A. 

2. 1.2. 3 Computations 

The fractional part of pl/2ir is multiplied by 2n. If the result is great- 
er than n, 2ir is subtracted. A flow chart is shown in Figure 4. 


27 






2.1.3 'ATNl 1 / * ATN2 1 


This multiple entrypoint function subprogram computes the arctangent of 
the ratio of two input quantities, whose individual signs determine the quad- 
rant of the output angle. The 1 ATN2 ' entrypoint is analogous to the FORTRAN 
ATAN2 function, in that the output angle lies in the range of -it to +ir. When 
the 'ATNl' entrypoint is used, the output angle lies in the range of 0 to Pit. 

2.1.3. 1 Examples of Usage 

The instruction 'ATNl 1 (-1 ,1 ) > A would cause the variable A to be as- 
signed the value of 7ir/4, while the instruction 1 ATN2 1 (-1 ,1 ) ■->• A would cause 
it to be assigned the value of -ir/4. 

2. 1.3. 2 Computations 

Tests are made to avoid division by zero. Depending on whether pi or p2 
has the greater magnitude, a base value of the output angle is computed from 
the expression ir/2 - atn(p2/abs(pl ) ) or the expression atn(abs(p!/p2) ) , where 
at]i and abs_ represent the HPL arctangent and absolute value functions. Tests 
on the signs of pi and p2 then are made so that the base value can be rotated 
into the proper quadrant. If pi = p2 = 0, the output angle is set equal to 
zero. A flow chart is shown in Figure 5. 
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2.1.4 1 HM . S 1 


The 1 HM. S 1 function converts time from seconds to hours, minutes, and 
seconds. The output is a single number in the form of HHMM.SS, where HH 
denotes the hours digits, MM denotes the minutes digits, and SS denotes the 
seconds digits. 

2. 1.4.1 Arguments 

pi = Input time, measured in seconds. 

2. 1.4. 2 Example of Usage 

The instruction 'HM.S 1 (36385.8) -» H in the calling routine would cause 
the variable H to be assigned the value 1006.258 (10 hours, 6 minutes, and 
25.8 seconds) . 

2. 1.4. 3 Computations 
See Figure 6. 
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Figure 6. 'HM.S' Function Subprogram Logic Flow 





2.1.5 'SECS' 


The 'SECS' function is the inverse of the ' HM. S ' function described in 
Section 2.1.4. It converts time from hours, minutes, and seconds (expressed 
in the HHMM.SS format) into seconds. 

2. 1.5.1 Arguments 

pi = Input time, in the HHMM.SS format. 

2. 1.5. 2 Example of Usage 

The instruction 'SECS' (1006.258) -» S in the calling routine would cause 
the variable S co be assigned the value 36385.8 (the number of seconds in 10 
hours, 6 minutes, and 25.8 seconds). 

2. 1.5. 3 Computations 
See Figure 7. 
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2.1.6 ' DOTP 1 


The 'DOTP' function subprogram computes the dot (scalar) product of two 
vectors . 

2. 1.6.1 Arguments 

pi = Address of first vector. 
p2 = Address of second vector. 

2. 1.6. 2 Example of Usage 

If vector ft resides in registers rl , r2, r3, and vector Ef resides in 
registers r4, r5, r6, then the instruction 'DOTP' (1 ,4) -> D in the calling 
routine would cause the value of A • Ef to be assigned to the variable D. 

2. 1.6. 3 Computations 
See Figure 8. 
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2.2 UTILITY SUBROUTINES 


2.2.1 'SXV' 

The subroutine 'SXV' calculates the product of a scalar with a vector. 

2. 2. 1.1 Arguments 

pi = Value of scalar. 

p2 = Address of input vector. 

p3 = Address of output vector. 

2. 2. 1.2 Examples of Usage 

The instruction ell 'SXV'(6. 8,0,3) in the calling routine would cause the 
values 6.6 rO, 6.8 rl , 6.8 r2 to be stored in r3, r4, r5. If vector J, resides 
in rll , rl2, rl3, and B resides in r7, r8, r9, the instruction c ISt 1 SXV ' ( ' DOTP ' 
(11 ,7) ,11 ,7) would cause the value of (A • B) ft to be stored in r7, r8, r9. 

2. 2. 1.3 Computations 


See Figure 9. 




Figure 9. 


l 

ip\ rj>x-*r{>5 
^>i r(^i+0^r(|D34O 
r(}aa4:0-*r(|D34a) 

~~r~ 

( ) 


'SXV' Subroutine Logic Flow 




2.2.2 1 VADD 1 / 1 VSUB ' 


Depending on which entrypoint is used, this subroutine will calculate 
either the sum or the difference of two vectors. 

2. 2. 2.1 Arguments 

pi = Address of first input vector. 
p2 = Address of second input vector. 
p3 = Address of output vector. 

2. 2. 2. 2 Example of Usage 

If vector A resides in rl , r2, r3 and vector F resides in r4, r5, r6, the 
instruction c jj 'VADD' (1 ,4,7) in the calling routine would cause the value of 
A + F to be stored in r7, r8, r9. The instruction cjj 'VSUB' (1 ,4,4) would 
cause the value of A - F to be stored in r4, r5, r6. The instruction 
c H 'VADD' (4,4,4) would cause the value of 2F to replace F in r4, r5, r6. 

The instruction c 11 'VSUB 1 (1 ,1 ,1) would cause zeros to be stored in rl , r2, r3. 

2. 2. 2. 3 Computations 
See Figure 10. 
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Figure 10. 'VADD 1 / 1 VSUB 1 Subroutine Logic Flow 






2.2.3 ' CRSP ' 


The 'CRSP' subroutine computes the cross product of two input vectors. 

2.2.3. 1 Arguments 

pi = Address of first input vector. 
p2 = Address of second input vector. 
p3 = Address of output vector. 

2. 2. 3. 2 Example of Usage 

If the vector A resides in rl , r2, r3 and the vector B" resides in r4, r5, 
r6, the instruction cit 'CRSP 1 (1 ,4,7) in the calling routine would cause the 
value of A x B to be stored in r7, r8, r9. The instruction c ll 'CRSP 1 (4,1 ,4) 
would cause the value of B x A to replace B in r4, r5, r6. 

2. 2. 3. 3 Computations 
See Figure 1 1 . 
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Figure 11. 'CRSP 1 Subroutine Logic Flow 
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2.2.4 1 ROT 7*1 ROT ' 


This multiple entrypoint subroutine is used for transforming the com- 
ponents of a vector between two Cartesian coordinate systems that are angularly 
displaced from one another by a rotation about their common origin. Let the 
two coordinate systems be designated F and G, and let the versor (unit quater- 
nion) qpg define the orientation of 6 with respect to F (see Appendix A). 

The 'ROT 1 entrypoint performs the computation symbolized by cfpg o Vp o qfp G -► V Q . 
The symbol Vp represents the value of a vector V in coordinate system F, and Vg 
represents the value of the same vector in system G. The I ROT entrypoint 
performs the inverse transformation, symbolized by qpg o Vg o fpg Vp. 

2.2.4. 1 Arguments 

pi = Address of input vector. 

p2 = Address of versor. 

p3 = Address of output vector. 

2. 2. 4. 2 Examples of Usage 

If the versor qpg resides in rO, rl , r2, r3 and the vector Vp resides 
in r4, r5, r6, the instruction ell l R0T'(4,0,7) in the calling routine would 
cause the vector Vg to be stored in r7, r8, r9. A subsequent instruction 
cM * I ROT 1 (7,0,7) would cause a copy of the vector Vp to replace Vg in r7, r8, 
r9. The instruction cl£ 'ROT* (4,0,1) would cause Vg to be stored in rl, r2, 
r3 (thus destroying the last 3 components of the versor qpg). 

2. 2. 4. 3 Computations 

The quaternion operation qpg o Vp o qpg ■+ Vg is equivalent to the matrix 
operation [R] Vp Vg, where 
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(HH+II-JJ-KK) 

2 ( I J+HK) 

2(1 K-HJ ) 


[R] ■ 

2 ( I J-HK) 

(HH-II+JJ-KK) 

2( JK+HI. ) 

, (26) 


2 ( I K+H J ) 

2( JK-HI ) 

(HH-II-JJ+KK) 



and where q FG =H+i I+jJ+kK. The inverse transformation 
qpG ° Vg o cfpg Vp is equivalent to the matrix operation [R]^ Vg ■> 7p, where 
[R] T is the transpose of [R]. By inspecting Equation (26), it can be seen that 
[R] T is obtained simply by reversing the sign of H, as reflected in the flow 
chart that is shown in Figure 12. 


44 


r^i+x 
r(Joi+rt** y 
r(jpi4a.W 
rUpi+\)->i 
r(|aa.4a>T 
r(|53L43>K 

(HM4lI-7T-kK > )K4A((lT4MK^'Y 4 (XK- -* r|j3 

(HH-1I 4 JJ-toOY 4a.((TK4HT > >2 4 (XT- HK^X > )-> r(Jo34l) 
(WH-II-TJ+klc^?+x((Xk4HX)X4aK-Hl^Y > )^ r(f>34;0 







2.2.5 * QDOT 1 / 1 QXQ * / 3 QXQC 1 / 1 QCXQ 1 


This multiple entrypolnt subroutine performs a variety of related quater- 
nion operations. As in Section 2.2.4, let qtpg represent a versor (unit 
quaternion) which defines the orientation of a Cartesian coordinate system 
Q with respect to another system F that has the same origin. Let the vector 
represent the angular velocity (measured in system 6) of system G with 
respect to system F. The 'QDOT' entrypoint is used to perform the operation 
symbolized by (q pG o u> G )/2 -* qp G , where q pG represents the derivative of q pG 
with respect to time. 

Let F and Q represent any two quaternions (not necessarily versors), and 
let p and Q represent their conjugates. The entrypoints 'QXQ', 'QXQC', and 
'QCXQ' (respectively) perform the operations symbolized by ¥ o Q + F, 

P o (I + R, and ? o Q -> R. 

2. 2. 5.1 Arguments 
'QDOT' Entrypoint : 

pi = Address of the orientation versor qp G . 

p2 = Address of the angular velocity vector ai^ . 

p3 = Address of the output derivative q"p G . 

'QXQ 7 'QXQC 7 'QCXQ' Entrypoints : 

pi = Address of the first input quaternion P\ 
p2 = Address of the second input quaternion Q\ 

p3 = Address of the output quaternion FT. 

2. 2. 5. 2 Examples of Usage 

Assume the variable A has the value of 25, that the versor qp G resides in 

rp 

r50, r51 , r52, r53, and that the vector Wq resides in r54, r55, r56. Then 
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the statement c ll 1 QDOT 1 (A+25,A+29,A+125) will cause the value of qp Q to be 
stored in rl50, rl51, rl52, rl53. If the versor resides in rO, rl , r2, r3, 
the statement c ll 1 QXQ 1 (A+25,0,4) wil I cause q^ to be stored in r4, r5, r6, 
r7. If this were followed by cU 'QXQC' (0,4,8) , the versor 

q GH 0 ^FH = ^GH 0 ^HF = ^GF wou ^ be storec ^ 1n r8 » r9 > rlO, rll. Then the 
statement c/l 1 QCXQ 1 (8,0,0) would cause a copy of the versor 

^GF 0 9 GH = 9 FG 0 q GH = q FH ( a ^ reac *y residing in r4, r5, r6, r7) to be stored 
in rO , rl , r2 , r3 . 

2. 2. 5. 3 Computations 

The computational sequence depicted by the flow chart in Figure 13 results 
from a straightforward application of the rules of quaternion algebra that are 
defined in Appendix A. 
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2.2.6 1 Q231 * / *Q21 3 * 


The digits 1,2,3 in the 'Q231 ' and 'Q213' entrypoint names identify coor- 
dinate-system rotations about the X.Y.Z axes (respectively) in two different 
Euler-angle sequences. Generally, rotation about the X axis is designated as 
roll, rotation about the Y axis is designated as pitch, and rotation about the 
Z axis is designated as yaw. Given an orientation versor q^g, the 'Q23V entry- 
point is used to calculate the values of the pitch, yaw, roll angles (taken in 
that sequence) that are required to rotate the axes of coordinate system F into 
coincidence with those of system G. The 'Q213' entrypoint is used to calculate 
the corresponding Euler angle set for a pitch, roll, yaw sequence. 

2.2.6. 1 Arguments 

pi = Address of the orientation versor. 

p2 = Address of the output array of Euler angles (measured in radians, and 
stored in the same order as the logical rotation sequence). 

2. 2. 6. 2 Example of Usage 

If the orientation versor under consideration resides in rO, rl , r2, r3, 
the statement c U 1 Q231 1 (0,4) would cause the first rotation (pitch) angle to 
be stored in r4, the second rotation (yaw) angle to be stored in r5, and the 
third rotation (roll) angle to be stored in r6. The statement cil '0213(0,0) 
would cause the rotation angles of a pitch, roll, yaw sequence to be stored in 
rO, rl , r2 (thus wiping out the first three conponents of the orientation 
versor) . 

2. 2. 6. 3 Computations 

The output of the ' Q231 1 / 1 Q21 3 ' subroutine is an ordered set of Euler angles 
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a, B , y. For a pitch, yaw, roll sequence the Euler angles are 

related to the 

, — . a A 

orientation versor q FG = H + i I + j J + k K by the equations 

(see Reference 

10) 


H = C C , C - S S, S , 

a ii y a p y 

(27) 

I = c c, s + s S , C , 

a ^ y a h y 

(28) 

J = s c,c + c s , s , 

a [3 y a y 

(29) 

and 


K = C S, C - S C„ S , 

a (a y a 3 Y 

(30) 

where = cos ' 2 a, = sin ^ •x, C = cos ' 2 3, etc. Solving equations (27) 

(30) for a , 3» y results in 


a (a + y) = tan 1 [(J + I )/ (H + K)], 

(31) 

: 2 (o - y) = tan' 1 [(J - I )/ (H - K)], 

(32) 

3 = tan' 1 [2 (HK + Uj/A 2 + Y 2 ], 

(33) 

where 


X = H 2 - I 2 + J 2 - K 2 , 

(34) 

and 


Y = 2 (JK - HI). 

(35) 


As reflected by the flow chart shown in Figure 14, Equations (31) - (35) are 
used to compute the Euler angles when the 'Q231' entrypoint is used. 

For a pitch, yaw, roll sequence, the orientation versor and the Euler 
angles are related by the equations 
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( 36 ) 


H = C 

C, 

c 

+ 

s 

S, 

s , 

a 


Y 


01 


Y 

I = C 

s, 

c 

+ 

s 

c, 

s , 

Ol 


i 


a 


Y 

J = S 

c. 

c 


C 

s. 

s , 

a 


Y 


a 


Y 

K = C 

c, 

s 


S 

s. 

C . 

01 


“y 


a 


Y 


The Euler-angle solutions in this case are 

(« - y) = tan- 1 [(J - K)/(H + I)], 
\s (ot + y) = tan 1 [(J + K)/(H - I)], 

p = tan" 1 [2 (HI - JK)/ + V 2 ], 

where 

X = H 2 - I 2 + J 2 - K 2 
and 

Y = 2 (IJ + HK). 


(37) 

(38) 

(39) 

(40) 

(41) 

(42) 

(43) 

(44) 


It can be seen that Equations (31) - (35) become identical with (40) - (44) 
when 

(a) -K is substituted for I, 

(b) I is substituted for K, and 

(c) the sign of y is reversed in the first set of equations. 

This characteristic of the two Euler sequences allows common logic to be used 
for the major portion of the computations, as illustrated in Figure 14. 
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Figure 14. ' Q231 ' / 1 Q21 3 1 Subroutine Logic Flow 
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2.2.7 '0313* 


Given an orientation versor cfpg, the ' Q31 3 ' subroutine calculates a set of 
Euler angles (a,3,y) that corresponds to a rotation of coordinate system F 
into coincidence with system G, assuming F is rotated sequentially about its 
Z, X, Z axes. 

2.2.7. 1 Arguments 

pi = Address of the orientation versor. 

p2 = Address of the output array of Euler angles (measured in radians, and 
stored in the same order as the logical rotation sequence). 

2.2.7 .2 Example of Usage 

Suppose that coordinate system F is a geocentric inertial system defined 
such that the first point in Aries lies on the Xp axis, and the earth's 
angular momentum vector lies along the Zp axis. Let G be another geocentric 
system defined such that the instantaneous position of some satellite lies 
on the Xg axis, and such that the angular momentum vector of the satellite 
CG lies along the Z Q axis. Let the versor qp G reside in rO, rl , r2, r3. Then 
the statement cjj 1 Q31 3 1 (0,4) would cause the right ascension of the ascending 
node of the satellite orbit to be stored in r4, the inclination of the orbit 
to be stored in r5, and the argument of latitude corresponding to the satellite 
position to be stored in r6. 

2. 2. 7. 3 Computations 

- AAA 

The Euler angles are related to the versor q FG =H+i I+jJ+kKby 
the equations 
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H = COS l 2 (a + y) COS g, 

(45) 

I = cos l 2 (a - y) sin »2 6, 

(46) 

J = sin l 2 (a - y) sin ; 2 P, 

(47) 

and 


K = sin x 2 (a + y) cos p 

(48) 

(see Reference 10). It follows that 


h* = tan 1 (v/II+JJ, v/RIMk), 

(49) 

l 2 (a + y) = tan" 1 (K,H) , 

(50) 


and 


l 2 (ot - y) = tan" 1 (J,I). 

A flow chart of the 'Q313' logic is shown in Figure 15. 
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Figure 15. 'Q313' Subroutine Logic Flow 



2.2.8 * 231Q V 1 213Q 


The '231Q7'213Q' subroutine performs the inverse function of the 'Q231'/ 
'Q213' subroutine described in Section 2.2.6. That is to say, given an order- 
ed set of Euler angles a, p, y, it calculates the associated orientation versor 
qp G . The '231Q' entrypoint is used when the rotation sequence is pitch, yaw, 
roll, and the ' 21 3Q ' entrypoint is used when the sequence is pitch, roll, yaw. 

2.2.8. 1 Arguments 

pi = Address of the input array of Euler angles (measured in raoians, and 
stored in the same order as the logical rotation sequence). 

p2 - Address of the orientation versor. 

2. 2. 8. 2 Example of Usage 

If rl , r2, r3 contain the pitch, yaw, and roll angles that define a 
certain orientation of the Orbiter's body axes B relative to its local-vertical 
frame G, then the statement ell * 231 Q 1 ( 1 ,0) in the calling routine would cause 
the orientation versor qf GB to be stored in rO, rl , r2, r3 (thus destroying the 
input values of the Euler angles, plus whatever value may have been stored in 
rO before calling the subroutine). If the variable A has been assigned the 
value 100, and if r4, r5, r6 contain the pitch, roll, yaw angles that define 
the orientation of system G with respect to another system F, the statement 
M l 213Q , (4,A+10) would cause the versor q^ to be stored in rllO, rill, rl 1 2 , 
rl 1 3 . 


2. 2. 8. 3 Computations 

The computational sequence defined by the flow chart in Figure 16 is a 
straightforward imp! mentation of Equations (27) - (30) and (36) - (39) from 
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Figure 16. * 231 Q ’ / ‘ 21 3Q 1 
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2.2.9 1 3 1 3Q 1 


The '313Q' subroutine performs the inverse function of the 'Q313* subroutine 
described in Section 2.2.7. Given an ordered set of angles a, 8, y that 
defines the orientation of a coordinate system G with respect to another system 
F in terms of sequential rotations of F about its Z, X, Z axes, the * 31 3Q ' sub- 
routine computes the associated versor cjpg 

2.2.9. 1 Arguments 

pi = Address of the input array of Euler angles (measured in radians, and 
stored in the same order as the rotation sequence). 
p2 = Address of orientation versor. 

2. 2. 9. 2 Example of Usage 

Let the F and G coordinate systems be defined as in Section 2.2.7. Let 
the right ascension of the satellite orbit reside in r4, the inclination of the 
orbit in r5, and the satellite's argument of latitude in r6. Then the state- 
ment c ll 1 31 3Q 1 (4 ,0) would cause the orientation versor q"pg to be stored in 
rO, rl , r2, r3. 

2. 2. 9. 3 Computations 

The computational sequence defined by the flow chart in Figure 17 is a 
straightforward implementation of Equations (45) - (48) from Section 2.2.7. 
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Figure 17. ' 31 3Q 1 Subroutine Logic Flow 
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2.2.10 1 IMATQ 1 


Given the inverse 


[R]' 1 = DO T - [T] 


11 

T 12 

T 13 

21 

T 22 

T 23 

31 

T 32 

T 33 


( 52 ) 


of the coordinate transformation matrix [R] defined by Equation (26) in Section 
2.2.4, v.ie 'IMATQ' subroutine computes the associated orientation versor 


A A A 

^FG q 0 + 1 ^1 + "i ^2 + ^ ^3' 


(53) 


2.2.10.1 Arguments 

pi = Address of the inverse transformation matrix, stored (columnwise) in 

the order T^, T 21 » T 31 , T 12 > T 22 # T 32* T 13’ T 23’ T 33’ 

p2 = Address of the output versor. 

2.2.10.2 Example of Usage 

Let F represent an ECI coordinate system, defined such that the first 
point in Aries lies on the Xp axis, and the earth's angular momentum vector is 
aligned with the Zp axis. Let G represent another geocentric coordinate 
system defined such that the X Q axis is aligned with the geocentric position 
vector R of a satellite, and such that the Zg axis is aligned with 7x7, where 
7 is the geocentric inertial velocity vector of the satellite. Suppose that 
7 and 7 are known in terms of their components in F, and that it is desired to 
evaluate qpg. 

The inverse transformation matrix can be expressed as 


where U^ X = R p / 1 R p | , Uj? Z = (R p x V p )/|R p x V p | , and uj? Y = U^ Z x U^ X . If R p re- 
sides in r!3, rl4, rl5, the statement c£l 'SXV' (1AT( l D0TP l (13J3)) ,13,0) will 

— GX — 

cause U p to be stored in rO, rl , r2. If V p resides in r!6, rl7, r!8, the 

statement c// 'CRSP' (13,16.6) ; c l£ * SXV * ( 1 A/~( 'DOTP' (6,6)) ,6,6) will cause U{? Z to 

be stored in r6, r7, r8, and the statement cM 'CRSP' (6,0,3) will then cause 
— TY 

Up to be stored in r3, r4, r5, thus completing the columnwise storage of [T] 
in rO through r8. The statement c il 1 IMATQ 1 (0,9) will then cause to be 
stored in r9, rlO, rll, rl2. 


2.2.10.3 Computations 

By comparing Equations (26) and (52), it can be seen that [T] can be 
expressed, in terms of the versor components appearing in Equation (53), as 


2 2 2 2 

( WW 2 ( q l q 2~ q 0 q 3^ ^ q l q 3 +q 0 q 2^ 


[T] = 


2 ^ q l q 2 +q 0 q 3^ ^ q 0"" q l +q 2~ q 3^ ^ q ? q 3" q 0 q l^ 

^ q l q 3~ q 0 q 2^ ^ q 2 q 3 +q 0 q l ) ^ q 0~ q l~ q 2 +q 3^ 


(54) 


2 2 2 2 

Taking cognizance of the facts that qg + q-j + + q^ = 1 , and that one (but 

only one) of the versor components can be given an arbitrary sign, four 
different solutions for the values of q^fk = 0,1 ,2,3) can be found in terms of 
the T..(i,j = 1,2,3). These solutions are given by the equations 

* O 


q 0 = v/l + T n + T 22 + T 33^ 2 


(54a) 


q l = ^ T 32 ‘ T 23^ 4 q 0 


(54b) 
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or 


or 


or 


r 


q 2 ' ^ T 13 “ T 31^ /4 

q 3 = < T 21 - T l 2 )/4 q 0 
q i = v/TT ~ r n - t^t^/2 

q 2 = (T 21 + T 12^ 4 q] 

q 3 = (T 13 + T 31 )/4 
% = ( T 32 " T 23^ /4 q] 

q 2 = 

q 3 = (T 32 + T 23^ 4 q 2 

q 0 = ( T i3 ■ T 3i^ /4 q 2 

q l = (T 21 + T 12 )/4 q 2 

q 3 = s/T^h ? r ^ 2 7 lr ^/2 
q o = ^ T 2 i " t i 2 )/ 4 q 3 
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(54c) 

(54d) 


(55a) 

(55b) 

(55c) 

(55d) 

(56a) 

(56b) 

(56c) 

(56d) 

(57a) 

(57b) 


i 


q l = ^ T 13 + T 31^ 4 q 3 


(57c) 


q 2 = {T 32 + T 23 )/4 q 3- (57d) 

Since any versor component (or as many as three of them) may turn out to 

be zero, numerical tests are necessary to determine which one of the alternate 

computational sequences defined by the preceding equations can be used in a 

2 2 2 2 

given case. It is known from the versor identity q^ + q^ + q 2 + q^ = 1 that 

at least one component must have a magnitude as great as and that none can 

have a magnitude greater than unity. As defined by the flow chart shown in 

Figure 18, the ' I MATQ ' logic tests the magnitude of each versor component in 

sequence until one is found that has a magnitude at least as great as k. As 

2 

soon as a component q^ is found that satisfies the inequa ty 4 q^ >_ 1 , it is 
arbitrarily given a positive sign, and the remainder of tr ... components are 
calculated by using the equations that contain 4 q^ in the denominator. 
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Figure 18. ' I MAT Q ' Subroutine Logic Flow 
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2.3 STATE VECTOR DERIVATIVES SUBROUTINE ('DERI VS') 

The 'DERI VS 1 subroutine computes the first derivative (with respect to time) 
of each quantity contained in an extended array of state variables that in- 
cludes, in addition to the translational and rotational state vectors defined 
in Section 1.2, such things as propellant expenditures and control torque 
integrals. The memory allocation table contained in Appendix D lists the 
specific quantities that are included in the Shuttle state variable array 
(registers r25-r49) and the payload state variable array (registers r50-r74). 

Containing mathematical models of the orbital flight environment and of 
all vehicle systems affecting the dynamics of the Shuttle and the payload, the 
' DERI VS ' subroutine represents the heart of the #TRAJ processor. 


2.3.1 Input Data 

2.3.1 .1 Argument List 

The ' DERI VS ' argument list consists of a single parameter (pi) which is 
assigned a value of 0 or 1 by the calling routine according to whether state 
variable derivatives are to be computed for the Shuttle or the payload. 


2. 3. 1.2 r-Registers 

The contents of the following r-registers (see Appendix D ) are used as 
input quantities in the ' DERI VS 1 calculations: 


Shuttle Computations 


ill 


Compi 


r23 , r24 
r25-r39 
r75-r86 
r91-r99, rl24 


Payload Computations 
(Pi = 1) 

r23, r24 

r50-r64, r74 

rlOO-rlll 

rl!6-rl24 
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2. 3. 1.3 Array Variables 


The appropriate 12x12 RCS/DAP response matrix (see Section 3) must reside 
in the HPL array variable R[*] (represented symbolically as [R] in the following 
flow charts) before calling ' DERI VS ' for Shuttle derivatives (i.e., with pi = 0). 

When SRM thrusting is to be simulated, the name of the disk file containing 
the appropriate SRM thrust table must have been assigned to file no. 3 and the 
array variable A[*] (see Figure 19d for contents) must have been initialized 
before calling DERIVS for payload derivatives (i.e., with pi = 1). The contents 
of the A[*] array (symbolically, [A]) are updated as required in 'DERIVS' by 
reading new thrust-profile coordinates from file no. 3 as the burn progresses. 

2.3.2 Output Data 

As indicated in Appendix D , the r-registers have been allocated so that 
the address of every state derivative can be found simply by adding 100 to the 
address of the corresponding state variable. Thus, Shuttle state derivatives 
are stored in r!25-rl49 when pi = 0, and payload state derivatives are stored 
in registers rl 50-rl 74 when pi = 1 . When pi = 0, the contents of r285 
(normally zero) will be set equal to 1 if 'DERIVS' finds the Shuttle RCS 
inadequate to maintain the commanded attitude. 

Sometimes the calling routine requires, in addition to the state variable 
derivatives, the angular velocity vector (wq G or ^) 0 f the local-vertical 
coordinate system. This „-ctor resides in the volatile registers rl-r3 when 
execution control is returned from 'DERIVS' to the calling routine. Another 
vector sometimes needed by the calling routine is the linear acceleration 
vector (Ag or a^), which resides in r7-r9 upon the return from 'DERIVS'. 
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2.3.3 Computations 


The computations performed by 'DERI VS' are defined by the flow chart shown 
in Figure 19a through 19g. For the most part, the computations are described 
in terms of logical symbols (as opposed to the HPL variable names and addresses 
that appear in the code.) The major portion of the computational logic applies 
equally to Shuttle and to payload calculations. To minimize confusion in 
following the flow of the logic, Shuttle symbols (e.g., Aq rather than a^) are 
used almost exclusively throughout the flow chart, even in those sections of 
the logic that apply uniquely to the payload. The only exceptions to this 
general rule occur when there is no Shuttle-related equivalent of a payload- 
related variable (e.g., i^ + , i 2+ , etc.). 

Insofar as practicable, the order of computation in the HPL code (Section 
C.l) is arranged to follow that shown in the flow chart. Even so. Appendix D , 
which shows the correlation between logical symbols and r-register numbers, is 
indispensable to the understanding of the code. 

Following is a summary of all the data registers vhat will or may be 
modified by an execution of ' DERI VS 1 : 

Shuttle Computations Payload Computations 

(Pi = 0) (Pi = 1) 


r0-rl8 

r0-r!8 

r87-r90 

rll2-rll5 

rl25-r!49 

r!50-rl74 

r285 

A[5]-A[8] 

A-Z 

A-Z 
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Figure 19a. ' DERI VS ' Subroutine Logic Flow 
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Figure 19b. ' DERIVS ' Subroutine Logic Flow (cont.) 
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Figure 19c. ' DERIVS ' Subroutine Logic Flow (cont.) 
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Figure 19e. ' DERI VS ' Subroutine Logic Flow (cont.) 
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Figure 19f. 1 DERI VS * Subroutine Logic Flow (cont.) 
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2.4 MAIN LOGIC FLOW 


The "main logic" of the #TRAJ link consists of nothing more than an 
instruction that causes the %RMAT link (Section 3) to be appended to #TRAJ. 


3. RESPONSE MATRIX COMPUTATION LINK («RMAT) 


The function of the %RMAT link is to compute matrices that describe the 
response of the Shuttle to DAP/RCS commands, and to store the matrices in 
designated disk files for subsequent use during trajectory integration. 

3.1 INPUT DATA 

3.1.1 Jet Force Table 

Basic inputs to the response matrix computations include the data shown in 
Table 1, which reside in the disk files named "$JFT" and "$JFTM". The "$JFTM" 
file contains the mnemonic identification symbols (shown in the second column 
of Table 1) for the Orbiter's 44 RCS jets, as defined in Figure 4. 3. 2. -2 of 
Reference 11. The "$JFT" file contains the body-fixed thrust components and 
the station coordinates of the thrust application point (which together define 
a torque component normal to the thrust line) for each thruster, along with a 
factor C (in the last column of Table 1) whose product with the thrust vector 
defines a component of torque parallel to the thrust line. This latter com- 
ponent of torque (parallel to the thrust line) is the result of RCS jet plume 
impingement on the exterior surfaces of the Crbiter. Table 1 results from 
combining the basic RCS jet data shown in Table 2 with the plume impingement 
data shown in Table 3, as explained in Reference 9. 

3.1.2 Jet Select Tables 

For each flight profile segment, the HFRMP user must specify which one of 
three basic combinations of thrusters is to be used for attitude and/or trans- 
lational control of the Shuttle. The available options are designated V 
(vernier jets), P (primary jets), and PZI (primary jets with +Z thrusters 
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Table 2. Basic RCS Thrust Data (Without Plume Impingement) 


1 ' 

f 


THRUSTER 

NO. 

THRUSTER 

ID. 

Hi 

wm 

HI 

■ 


RE 

(IN) 


F2F 

-879.4 

-26.2 

119.9 

306.72 

14.65 

392.96 


F3F 

-879.5 

0.0 

122.7 

306.72 

0.00 

394.45 


FI F 

-879,4 

26.2 

119.9 

306.72 

-14.65 

392.96 


F1L 

-26.3 

8?3.6 

18.2 

362.67 

-69.50 

373.73 


F3L 

-21.0 

870.3 

0.5 

364.71 

-71.65 

359.25 


F2R 

-26.3 

-873.6 

18.2 

362.67 

69 . 50 

373.73 

7 

F4R 

-21 .0 

-870.3 

0.5 

364.71 

71.65 

359.25 

8 

F2U 

-32.3 

-11.7 

874.4 

350.93 

14.39 

413.46 

9 

F3U 

-31.9 

0.0 

873.5 

350.92 

0.00 

414.52 


F1U 

-32.3 

11.7 

874.4 

350.93 

-14.39 

413.46 


F2D 

-28.0 

-616.4 

-639.5 

333.84 

61.42 

356.95 


TID 

-28.0 

616.4 

-639.5 

333 . 84 

-61.42 

356.95 

13 

F4D 

-24.8 

-612.6 

-639.4 

348.44 

66.23 

358.44 

14 

F3D 

-24.8 

612.6 

-639.4 

348.44 

-66,23 

358.44 

15 

R3A 

856.8 

0.0 

151.1 

1555.29 

137.00 

473,06 

15 

R1A 

856.8 

0.0 

151.1 

1555.29 

124.00 

473,06 

17 

L3A 

856.8 

0.0 

151.1 

1555.29 

-137.00 

473.06 

18 

LI A 

856.8 

0.0 

151.1 

1555.29 

-124.00 

473.06 

19 

L4L 

0.0 

870.5 

-22.4 

1516.00 

-149.87 

459.00 

20 

L2L 

0.0 

870.5 

-22.4 

1529.00 

-149.87 

459.00 

21 

L3L 

0.0 

870.5 

-22.4 

1542.00 

-149.87 

459 . OO 

22 

L1L 

0.0 

870.5 

-22.4 

1555.00 

-149.87 

459.0;" 

23 

R4R 

0.0 

-870.5 

-22.4 

1516.00 

1 49 . 87 

459.0i 

24 

R2R 

0.0 

-870.5 

-22.4 

1529.00 

149.87 

459.00 

25 

R3R 

0.0 

-870.5 

-22.4 

1542.00 

1 49 . 87 

459,00 

26 

R1R 

0.0 

-870.5 

-22.4 

1555.00 

1 49 . 87 

459.00 

27 

L4U 

0.0 

0.0 

870.0 

1516.00 

-132.00 

480.50 

28 

L2U 

0.0 

0.0 

870.0 

1529.00 

-132.00 

480.50 

29 

LIU 

0.0 

0.0 

870.0 

1542.00 

-132.00 

480.50 

30 

R4U 

0.0 

0.0 

870.0 

1516.00 

132.00 

480.50 

31 

R2U 

0.0 

0.0 

870.0 

1529.00 

1 32 . 00 

480.50 

32 

R1U 

0.0 

0.0 

870.0 

1542.00 

1 12.00 

480.50 

33 

L4D 

170.4 

291.8 

r801 .7 

1516.00 

-111.95 

437.40 

34 

L2D 

170.4 

291.8 

-801.7 

1529.00 

- 111.00 

440 00 

35 

L3D 

170.4 

291 .8 

-801.7 

1542.00 

-110.06 

442.60 

3b 

R4D 

170.4 

-291 .8 

-801 .7 

1516,00 

111 .95 

437.40 

37 

R2D 

170.4 

-291.8 

-801.7 

1529.ro 

111. 00 

440.00 

38 

R3D 

170.4 

-291.8 

-801.7 

1542.00 

1 10.06 

442.60 

39 

F5R 

-0.8 

-17.0 

-17.6 

324..,., 

59.70 

350,1.. 

40 

F5L 

-0.8 

17.0 

-17.6 

324.35 

-69.70 

350.12 

41 

R5R 

0.0 

-24.0 

-0.6 

1565.00 

149.87 

459.00 

42 

L5L 

0.0 

24.0 

-0.6 

1565.00 

- 1 49 . 87 

459.00 

43 

R5D 

0.0 

0.0 

-24.0 

1565.00 

118.00 

455.44 

44 

L5D 

0.0 

0.0 

-24.0 

1565.00 

-118.00 

455.44 
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Table 3. Force and Moment Increments Due to Plume Impingement 
(CG @ STA 1076.7, BL 0, WL 375.0) 


THRUSTER 

ID 

A FORCE (LB) 

A MOMENT (FT-LB) 

F x 

f y 

mtm 


l y 


L4L 

3.34 

-1 .66 

27.17 

-795.30 

762.42 

144.46 

L2L 

3.34 

-1.66 

27.17 

-795.30 

762.42 

144.46 

L3L 

3.34 

-1.66 

27.17 

-795.30 

762.42 

144.46 

L1L 

3.34 

-1.66 

27.17 

-795.30 

762.42 

144.46 

L4U 

0.09 

76.30 

2.21 

1859.90 

94.21 

-2932.80 

L2U 

0.09 

76.30 


1859.90 

94.21 

-2932.80 

LIU 

0.09 

76.30 

2.21 

1859.90 

94.21 

-2932.80 

L4D 

40.19 

26.59 

225.66 

-2475.50 

8469.00 

-645.46 

L2D 

40.19 

26.59 

225.66 

-2475.50 

8469.00 

-645.46 

L3D 

40.19 

26.59 

225.66 

-2475.50 

8469.00 

-645.46 ; 

R4R 

3.34 

1.66 

27.17 

795.30 

762.42 

-144.46 | 

R2R 

3.34 

1 .66 

27.17 

795.30 

762.42 

-144.46 

R3R 

3.34 

1.66 

27.17 

795.30 

762.42 

-144.46 

R1R 

3.34 

1.66 

27.17 

795.30 

762.42 

-144.46 

R4U 

0.09 

-76.30 

2.21 

-1859.90 

94.21 

2932.80 

R2U 

0.09 

-76.30 

2.21 

-1859.90 

94.21 

2932.80 

R1U 

0.09 

-76.30 

2.21 

-1859.90 

94.21 

2932.80 

R4D 

40.19 

-26.59 

225.66 

2475.50 

8469.00 

645.46 

R2D 

40.19 

-26.59 

225.66 

2475.50 

8469.00 

645.46 

R3D 

40.19 

-26.59 

225.66 

2475.50 

8469.00 

645.46 


(Take', from Reference 13) 
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inhibited). Corresponding to each of these options is a jet-select table 
(Tables 4-6) which identifies the particular jet or combination of jets that 
is to be fired in response to each of the six RCS translational acceleration 
commands (+X, -X, +Y, -Y, +Z, -Z) and the six rotational acceleration commands 
(+R0L, -ROL, +PCH, -PCH, +YAW, -YAW). The jet select tables, which reside in 
the disk files named "$JSV", "$JSP", and "$JSPZI", are not routinely available 
for modification by the user. However, the HFRMP software system includes a 
jet-select editing processor that makes it possible to update the tables or to 
provide additional options with little difficulty. 

As indicated in Table 4 by the absence of any jet designations for the 
execution of translation commands, the V option (vernier jets) can be used only 
for rotational control. The P option (Table 5) and the PZI option (Table 6) 
can be used for translational and/or rotational control. 

In the PZI option, no jets are fired that would expel propellant directly 
upward with respect to the Orbiter body. Translational acceleration in the 
downward direction, if commanded, is achieved (at a comparatively high propel- 
lant cost) by firing the +X and -X thrusters simultaneously. The cant angles 
of the +X and -X jet thrust lines produce a small net acceleration in the +Z 
(downward) direction. This option normally is used only when the Orbiter is 
maneuvering in the near vicinity of a payload that must be protected from jet 
plume impingement. 

3.1.3 Shuttle Mass Properties 

* 1 * 

At %RMAT execution time, the user-defined mass properties of the Shuttle 


^Which remain constant during HFRMP trajectory integration. 
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Table 


Vernier (V) Jet Select Table, Disk 
File "$JSV" 













Table 5. Primary (P) Jet Select Table, Disk 
File "$JSP'' 


CMD 


THRUSTERS TO 

BE FIRED 

1 

2 

3 

■a 



7 

8 

+X 

R1A 

LI A 







-X 

F2F 

FI F 







+Y 

FI L 

L4L 







-Y 

F2R 

R4R 







+Z 

F3U 

L4U 

R4U 






-Z 


F2D 

L4D 

L2D 

R4D 

R2D 



+ROL 

L4D 

R4U 







-ROL 

L4U 

R4D 







+PCH 


F2D 

L4U 

R4U 





-PCH 

F3U 

L4D 

R4D 






+YAW 

FI L 

R4R 






■ 

-YAW 

F2R 

L4L 






1 











1 


: irirusters Inhibited (PZI) 
e. Disk File "$JSPZI" 


BE FIRED 












reside in the disk file named "1*", where they are stored by the data-base 
editing processor #DBED. The "1*" file contains 32 scalar quantities, of 
which %RMAT makes use of the following: 


Item No. 
1 
2 


Description 

Shuttle gross weight (lb) 

I, 


l XX 


3 

T YY 

4 

I zz- 

5 

V 

6 

! zx ' 

7 

J XY J 

8 

STA ' 

9 

BL • 

10 

WL 

32 

%RMAT 


moments of 
inertia 


products of 
inertia 


(slug-ft^) 


station coordinates 
of CG (in) 


Item 32 in the "1*" file is a flag that is tested immediately after entry into 
%RMAT to determine whether a re-computation of the response matrices is 
necessary. It is set equal to zero during initialization of the program disk, 
and thereafter (by the appropriate editing processor) when the contents of 

"$JFT", or any jet-select table are changed in any way. It is set equal 
to 9 by %RMAT upon completion of the computation and storage of the response 
matrices . 


3.2 OUTPUT DATA 

The output of the %RMAT link consists of eight 12x12 matrices, each corre- 
sponding to a particular jet select table and a particular mode of RCS cross- 
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coupling compensation, and each stored In a separate disk file as indicated in 
Table 7. The compensation modes available to the HFRMP user are designated 
NONE (no compensation), ROT (rotational compensation only), and FULL (rota- 
tional ard translational compensation). In any given flight profile segment, 
the HFRMP user may choose any jet-select option in combination with any com- 
pensation mode, except V with FULL. Full compensation is impossible with the 
vernier jets, and the compensation mode is internally defaulted to ROT when 
such a combination is specified by the user. 

Tables 8-10 show typical response matrices for the P jet-select option 
and increasing degrees of cross-coupling compensation. For purposes of illus- 
tration, each matrix was transposed and then partitioned into two 12x6 matrices 
so it could be printed conveniently on a single page. Each column of the 
(untransposed) response matrix corresponds to a particular translational or 
rotational command. Rows 1-3 contain the body-axis components of the steady- 
state linear acceleration, and rows 4-6 contain the components of angular accel- 
eration. Rows 7-12 contain RCS propellant consumption rates, broken down 
according to source (forward, aft left, or aft right tanks) and control function 
(translation or rotation). The computations that produce the results shown in 
Tables 8-10 are described by the flow chart shown in Figures 20a through 20g, 
wherein the symbols [U], [R], and [F] represent the uncompensated, rotationally 
compensated, and fully compensated response matrices. 

3.3 COMPUTATIONS 

It will be noted that six small rectangles, each enclosing a pair of 
acceleration components, run diagonally across the upper partition of the 
matrix shown in Figure 8. The components thus enclosed represent the uncompen- 
sated "principal response" of the Shuttle to each of the 12 translation and 
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Table 8. Uncompensated Response Matrix for Jet Select Option P, Disk File "*P 


E 

3 


V) 

c 1) 
u 


f — X3 
O C 
CJ ►— 


(/> 
C QJ 

F u 

3 * r ~ 
r — 

O C 
CO 


- vm to m un to o« no cn O r- cm 


• tM o 

LU 

Ol 

S. 

Q 




r- f- U i L»> Ul'lf-. f.tJQi 

o o co vo o o ‘ o ro o ol 

O o on co o r ;.r> O o Oi 

O O O O O - o CO O O O 1 

OQOOOOOOOO 

ooooooooool 


ro ro 

C\l CM 


SI 


ro 


u? 

cn 

, — 



1 — 

X 

0> 

cn 

CO 

CO 


1 — 

i — 

1 — ■ 

CO 

to 

to 

CM 

cn 

r— 

1 — 

m 

CM 

cn 

o 

1 — 

O 

, — 

1 — 

CM 


o 

o 

o 

O 

CO 

o 

o 

o 

o 

o 

CM 

r— 

o 

o 

o 

O 

CO 

o 

CO 

o 

CO 

o 

O 

o 

to 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

p 

a 


CO 

to 

«=£ 

Q_ 

c 

—I 

ID 

CD 


r — cm cn cn r - M 

O O CD CD O , 
oouj^od 

o o o o o o 
o o o o o o 


o o o o o o 


CO CO 

co co| 
cr > Lof 


co cn co co 
o o o 
o a ro ro 
o o o o 
o o o o 


o olo o o o 


r~ N Ol Ol 
CO CD CO CD 

co co ro on 
co O O' 
o o o o 


o o o o 


CM M 
LO 


*+ x i — cn CTi 

CO CO CO co CO CO 

N N O- CO n 

r-. o o 

o o o o o o 


5 

o 


o o o o o o 


o o 
o o 
o o 
o o 

O O' 


o o 


CO co 
I CM CM 

o ol 
CO co| 

CM CM 


o o 


On ' t^OONN 

Or-aicnooNn 

oococooooo 

oororooooo 

oooooooo 


oooooooo 


CJ 

o 

CC 

Q£ 

<c 

UJ 


ra 


M- 

co cn. 

to CM 

x coj 

CM CM 


CO o 


cocnOcncncTiCOCMcncn 
COCOr- ^ cc co cn cc CD CD 

rorouncorooncocuroco 

ooocMroonocooo 

OOOf—OOOOOO 


oooooooooo 


t "3 

5 ; 

C_) 


X >- >- 

+ 1+1 


NNM - JX 1 j 33 
+ t o O O co «a: cc 
o:a:a_a_>~>- 
+ 1+1+1 


I ' 1 

r - CM ■') V II ) UJ r •>, co o f - CM 




1 ■ 
31 


CO 

o 

o 

kO 

vO 

V. ) 

Cj 

CJ 

CJ 

o 

o o 



CO 



O 

CJ 

CJ 

CO 

o 

o 

1 — • 

r - 


r— 

*• •* o 




bT 

O 

o 

o o 

to 

LJ 

r • 

X 

r - 

r n, 

r ^ cj 




o 

o 

o 

o 

o 

CJ 

o 

O C 

- J 

C" r r, r™» 





o 

CO 

CO 

CO 

o o 

r— 

i— 

i — 

1 — 

i — o 



i — 




• 




» 

* 

♦ 





u. 


o to 

O CO 

o 

o 

to 

ro 

CO 

on 

to o 


o 

«i 













07 













































o 

h~ 














CO 

u_ 


o 

o 

o o 

o 

cn 

rn 

CJ 

to 

o 

Cj C . 1 



UJ 

X 

o o 

o 

o o 

o 

1 — 

i — 

r— 


(O- f~* 




mZ 

o 

o 


o 

o 

o 

r«. 

X 

r- 


O tx 


d. 



o 

o 

o 

CO 

o 

o 

o 

CO 

O 

CO 

o o 




J— 

o 

o 

o 

O O C,T 

1-“- 

r— • 



o 

* 

O 

LU 


• 







t 

< 




*— i 

<1. 


a o o 

o o 

o 

CO 

on 

CO 

ro 

o to 
















to 

< 













































o 





























Ol 




o 

o 

o 

o 

o 

o 

o o 

o 

o 

co o 

u; 




o 

C 5 

cn 

o o 

o 

o 

o 

CM 


H* 


O 


o 

o 

o 

o 

o 

o o o 



<; 


3: 

o 

o o 

o 

o 

cn 

o 

o 

i — 

o 

o o 

a 1 


u_ h- 

o 

o 

o 

o o 

o 

o 

o 

CM 



2, 




o 

o 

o o o 

o 

o 

o 

CO 

to 

t o to 

o 






























Cl. 















s: 















Z 3 















on 


37 


o 

o o o o 

cn 

< > 

C ’ 

to 

to 

O O’ 

z: 


to 

m: 

1 — 

o o 

r— 

f— 

CM 

cn 

CO O 

o 

c o 

o 


KH 


f" 

o 

o 


X 

M- 

c> o o 

o 

c; £n 

co 


q: 

o 

CO 

o o 

o 

( — 

o 

O 

o 

o 

cn o 


.J 


i — 

1— 

CO 

o 

1 — 

r— 

CM 

o o o 

o 

o o 


O 

F— 


• 

• 

♦ 




* 

• 





QC 

U_ 


CO 

o 

o 

on 

ro 

to 

o 

o o 

'O 

i'n 

c 

►— 





























-J 

o 














LO 

a. 

CO 

t 


O O CO 

o 

o 

cn 

o o 

o 

o o o 

o 

1 

LU 


r— - 

o 

f— 

o 

f— 

CM 

o 

o 

o 

o 

o o 

or 

<1 

—1 


x o 

X O X 

*3- 

o 

o o 

o 

o o 

a. 



ex. 

o o 

o 

o o 

r~~~ 

o o o 

cn 

o o 


O 

t 

h- 


o 

r ~; 

o 


CM 

cn 

o o 

o 

cn o 



< 


ro 

Q 

on 

o on 

to 

o 

o 

cn 

o 

o o 


<c 















-J 















to 






























oc 



o 

O O Co 

o 

m 

o 

C) 

o 

o o o 




m: 

o 

CM 

f — 

r — 

r~* 

CM 

o 

C n 

o 

o 

o o 



Q 


o 




r^. 

*^r 

o 

i i 

1 — 1 


o o 



dc *=C 

o 

1 — 

o o 

o 

t— 

o o 

o 

to 

o o 



Lu 

h— 

O CM 



1 — 

t^l 

o 

o 

c'3 

o 

o o 





o 

to 

ro 

on 

O') 


o o 

o o 

0 

0 


Q 



X 

X 


>- 

r-M m 

1 


X 

X 



i: 



+ 

1 

+ 

1 

-t- 

1 

o o 

CO 

CO 

«a; 


co 









o: 

a. 

a. 

a. 

>- >- 











+■ 

l 

+ 

t 

4 1 


o 

-M 

'a 

t, 

+j 

r> 

b 


87 






Transposed and partitioned for purposes o 



Table 9. Rotationally Compensated Response Matrix 7 for Jet Select Option P, Disk File "*PR 
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ansposed a^d partitioned for purposes of illustration. 




Table 10. Fully Compensated Response Matrix for Jet Select Option P, Disk File 
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Transposed and partitioned for purposes of illustration. 




rotation commands. All other components of acceleration represent "extraneous 
responses", commonly referred to as "cross-coupling effects", that arise 
from the canting of certain jet centerlines away from the body axes, uneven 
moment arms, and plume impingement on the exterior surfaces of the Orbiter. 

Thu details of the uncompensated response matrix calculations are defined in 
Figures 20b and ’Oc. 

The process of cross-coupling compensation can be explained in the follow- 
ing general terms. Assume that the jets activated by a given "primary command" 

P fire continuously for a long period of time At. The effects of extraneous 
accelerations are nullified by intermittent firings of jets activated by 
"compensating commands". Let the accumulated firing time of the jets activated 
intermittently by any particular compensating command K be represented by 6t K> 
The ratio = st^/At is referred to as the "duty cycle" of the (jets activated 
by the) compensating command. 

The uncompensated response to the primary command (a column in the uncom- 
pensated response matrix) we represent with the symbol [U] p e [(U-| p ) (U 2 p ) ... 
(Ui 2 p )] . The compensated response (a column in the compensated response 
matrix is defined by the expression. 

K 

[C] p = [U] p + z y k [u] k , 

where the summation includes all the compensating commands that are required 
to nullify the extraneous accelerations. 

Rotational compensation, which involves the nullification of only the 
extraneous angular accelerations (see Table 9), is carried out in two phases. 

The rotation commands themselves are rotationally compensated in the firs' 
phase (Figure 20d), which is of necessity an iterative process. Success! 


90 


compensations introduce new extraneous angular accelerations, whose magnitudes 
must diminish progressively if convergence is to be realized. Non-convergence 
of this process implies that (given the input mass properties and the jet- 
select table under consideration) the attitude of the Shuttle can not be 
control led , and an appropriate warning message is output. 

The translation commands are rotationally compensated in the second phase 
of rotational compensation (Figure 20e). This is a non-iterative process that 
is simplified by use of the results from the first phase. 

Full compensation is achieved by nullifying the extraneous linear acceler- 
ations that remain in the rotationally compensated matrix. Again, the process 
is carried out in two phases. Starting with tne rotationally compensated 
matrix as an input (which eliminates the need to worry about compensation 
commands producing extraneous angular accelerations) the first phase (Figure 
20f) is devoted to the translational compensation of the translation commands. 
This is an iterative process that is essentially identical to that defined in 
Figure 20d. Likewise, the second phase of translational compensation (Figure 
20g) is essentially identical to that shown in Figure 20e. 

The computations are repeated for each of the three jet select tables, 
except that full compensation is not attempted with the vernier jets. Appro- 
priate warning messages (if any) are stored in the desk files along with the 
response matrices. 

The last executable statement in the "RMAT link causes the code of the 
I’jTNIT link (Section 4) to be appended to that of the base link (#TRAJ), in the 
region of computer memory formerly occupied by the foRMAT code. The HPL get 
command is used for this purpose, which causes all data registers used by «RMAT 
to be de-al located (erased) before execution control is handed over to loTNIT. 
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Figure 20a. I5RMAT Logic Flow 
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Figure 20g. %RMAT Logic Flow (cont.) 
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4. TRAJECTORY INITIALIZATION LINK (35TNIT) 


The function of the %TNIT link is to initialize calculator memory in prep- 
aration for trajectory integration. Elements of the user-supplied data base 
(Appendix E) are read from appropriate disk files (where they will have been 
stored by a previous execution of the #DBED processor) and processed as neces- 
sary to define the constants and the initial state of the Shuttle/payload system, 

4.1 FUNCTION SUBPROGRAMS 


4.1.1 'JD' 

The 'JD 1 function subprogram computes the Julian day number corresponding 
to a given year, month, and day of the Gregorian (civil) calendar. 


4.1 .1 .1 Argument List 


pi = year 
p2 = month 
p3 = day 




> 


integers 


4. 1.1. 2 Example of Usage 

The instruction 'ID' (1980,4,2) -» D would cause D to be assigned the value 
2444332 (the number of the Julian day commencing at Greenwich noon on 2 April 
1980). 


4. 1.1. 3 Computations 

Figure 21 is a flow chart of the 'JD' computational procedure, which was 
taken from Reference 14. 
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Figure 21. JO Function Subprogram Logic Flow 







4.1.2 'PLOTYP ' 


Given the user-assigned mnemonic symbol which designates the desired type 
of data plot, 'PLOTYP' returns the appropriate numeric code for internal use. 

4.1 .2.1 Argument List 

None. The mnemonic symbol must be assigned to the string variable P$ 
before 'PLOTYP' is executed. 

4. 1.2. 2 Example of Usage 

If the character string "CPLV" resides in P$, the instruction 'P LOTYP' -> I 
will cause I to be assigned the value 2. 

4. 1.2. 3 Computations 
See Figure 22. 
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4.1.3 ‘UNIT 1 


Given the user-assigned mnemonic symbol which designates the unit of 
distance for data plots, the 'UNIT' subprogram re.turns the appropriate conver- 
sion factor (the number of feet in the designated unit). 

4 .1.3.1 Argument List 

None. The mnemonic symbol must be assigned to the string variable P$ 
before 'UNIT' is executed. 

4. 1.3. 2 Example of Usage 

The instructions 11 NMI" -*■ P$; 'UNIT' •> C would cause C to be assigned the 
value 6076.115. 

4. 1.3. 3 Computations 
See Figure 23. 
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4.1.4 'MONTH' 


Given the abbreviated name of a month, the 'MONTH' function subprogram 
returns the appropriate month number. 

4 .1.4.1 Argument List 

None. The abbreviated month name must be assigned to the string variable 
P$ before MONTH is executed. 

4. 1.4. 2 Example of Usage 

If the character string " AUG" resides in P$, the instruction 'MONTH' •* M 
will cause the number 8 to be assigned to the variable M. 

4. 1.4. 3 Computations 
See Figure 24. 
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4.1.5 'C' 


Given the size of a square matrix stored columnwise in a one-dimensional 
array, and the row and column indices of one of its elements, the 'C' function 
returns the relative address (in the one-dimensional array) of the matrix 
element. 


4.1 .5.1 Argument List 

pi = Row index of matrix element. 

p2 = Column index of matrix element. 

p3 = Matrix size (number of rows = number of columns). 

4. 1.5. 2 Example of Usage 

Assume the 3x3 matrix [T] is stored columnwise in registers r8 - rl6, as 
shown below: 


Relative 

Address 

Register 

No. 

Matrix 

Element 

0 

r8 

T l,l 

1 

r9 

T 2,l 

2 

rlO 

T 3,l 

3 

rll 

T l,2 

4 

rl2 

T 2 ,2 

5 

r!3 

T 3,2 

6 

r14 

T l,3 

7 

rl5 

T 2,3 

8 

r!6 

T 3 ,3 


The instruction r(8 + 'C 1 (1 ,2,3) ) A would cause the value of T-j ^ to be 
assigned to the variable A. 
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Figure 25 'C' Function Subprogram Logic Flow 
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4.1.6 IT 

Given the row and column indices of an element in a symmetric matrix 
whose upper triangular form is stored columnwise in a one-dimensional array, 
the ' T ' function returns the relative address (in the one-dimensional array) 
of the matrix element. 

4. 1.6.1 Argument List 

pi = Row index of matrix element. 
p2 = Column index of matrix element. 

4. 1.6. 2 Example of Usage 

Assume the upper triangular form 

*1,1 J 1 ,2 *1,3 

1 2,2 *2,3 

*3,3 - 

of the symmetric 3x3 matrix [I] is stored columnwise in registers rl0-r!5 as 
indicated below: 


Relative 

Address 

Register 

No. 

Matrix 

Element 

0 

rlO 

*U 

1 

rll 

1 1 ,2 =I 2,1 

2 

rl2 

*2,2 

3 

rl 3 

*1 ,3 = 1 3 ,1 

4 

rl4 

*2,3 = 1 3 , 2 

5 

rl5 

*3,3 


The instructions 10 + 1 T * (1 ,2) •>■ H; rH -> X would cause the value of 1^ 2 = ^ 
to be assigned to the variable X. 

4. 1.6. 3 Computations 
See Figure 26. 
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4.2 MATRIX DIAGONALIZATION SUBROUTINE ( 'DIAG' ) 


Given a symmetric matrix [M] of order C, the 'DIAG' subroutine solves the 
equation 

[R] [M] [R] T = [D] 

for [R] and [D], where [D] is a diagonal matrix (one whose off-diagonal elements 
are all zero, within some specified tolerance). 

4.2.1 Argument List 

pi = Order (size) of the input matrix [M]. 

p2 = Address of the input matrix [M], stored columnwise in upper tri- 
angular form. 

p3 = Address of the die onalized output matrix [D], stored columnwise 
in upper triangular form. 

p4 - Address of the square (pi x pi) output matrix [R], stored columnwise. 
p5 = Tolerance on the maximum squared value of D. . , where i/j (this input 

» 5 «J 

is optional; if not supplied by calling routine, it will be calculat- 
ed by 'DIAG'). 

4.2.2 Example of Usage 

Assume that the upper triangular form of the Shuttle's inertia tensor 
[I]g (referenced to the body coordinate system B) is stored columnwise in 
registers rlO - rl5 as illustrated in Section 4. 1.5. 2. Then the instruction 
cl£ 'DIAG' (3,10,10,1 ) would cause [I]g to be replaced by [I] p (the same inertia 
tensor but now referenced to the coordinate system P, whose axes are the 
principal axes of inertia). Upon return from 'DIAG 1 , the registers rl - r9 
would contain the 3x3 coordinate transformation matrix [R] that satisfies the 
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equations 


[R] A p = A b 
and 

[R]' 1 A b = [R] T Ag = Sp, 
where A Is any arbitrary vector. 

4.2.3 Computations 

The method used to find the matrices [D] and [R] is of the type known as 
the Jacobi method (Reference 15). This consists of a series of matrix rotations 
of the form 

[D] = ([R],[R ] 2 • - • ) CM] ( - - • [R] 2 T [R]/) 

where each [R]^ is selected so as to cause one pair of off-diagonal elements 

J. L_ 

to be zero after the k tn rotation. An initial trigger level, < 5 , is computed 
from 

6 = 'J a /C 

where a is the sum of the squares of all off-diagonal elements and C is the 
order of [M] . The square of each off-diagonal element of [M] is compared to 
6 . When an element is encountered whose square is larger than < 5 , the matrix 
[M] is rotated so as to cause that off-diagonal element to be zero. 

For example, suppose that the square of the (i,j)th element of [M] is 
greater than < 5 , then [M] is rotated by [R]^ to get [M]' where 

M'[i,j] = M'[j,i] = 0 
The matrix [R]^ has the form 
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R-j[i »i] = R-j[j .j] = cos <j> 

R-j[i .j] = - sin 4 . 

R-j Cj *i 3 = sir) 4 

and all other elements of [R]^ are identical to the unit matrix [1]. Rotation 
with [R]i , i .e. , 

CM]' - [R], CM] [R]{ 
gives rise to the elements 

2 

M'[i,i] = M[i,i] cos <p + 2 M[i,j] sin ^ cos + 

2 

M' [j ,j] = M[i,i] sin <f> - 2M[i,j] sin cos 4 , + 

M'[i,j] = M ' [ j , i ] = (M[j ,j] - M[i,i]) sin 4 , cos 

+ M[i ,j] (cos 2 (j) - sin 2 <p) 

from which it can be seen that 

M'[i ,j] = M’[j,i] = 0 

provided <(> is selected such that 

tan 2 (j> = 2M[i ,j]/(M[i ,i] - M[j,j]) 

After this rotation, scan of the off-diagonal elements continues using 
now, however, [M]'. It is true that any one rotation will cause other off- 
diagonal elements which were zero to become non-zero, but the trend is to re- 
duce all off-diagonal elements to small numbers. 

When the scan of all off-diagonal elements has been completed, the trigger 
level 6 is reduced by dividing it again by C, and the process is repeated 


2 

sin 4 , 
2 

M[j »j] cos $ 

<P 
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until s is less than a pre-set tolerance. The programmer may include this 
tolerance in the argument list or he may let ' DIAG * set the tolerance for him. 
This process is shown in detail in Figures 27(a) and 27(b). 


116 




« C • OPOfdOP INPUT NATHI* [h] 
pi * AOOMT ts OP INPUT nPTPI* 
pi or CHAooNAUiro output mptpi* to] 

pH»POOP*1* OP OUTPUT TPANSFOPMATION MPTHIK [P] 
P*r « • CONVtROCNCf TOLSRPNCt (OPTION*!. INPUT ) 



(W'CKfr.pl.KCplJ) ) 


T^TV 

D~- 


o * w ** 

••i 


I ♦ A 


"*■1 




Figure 27a. ' DI AG ' Subroutine Logic Flow 
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Figure 27b. 'DIAG' Subroutine Logic Flow (cont.) 
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4.3 MAIN LOGIC 


The aITNIT computations are described by the flow chart contained in 
Figures 28a through 28g. The memory allocation table in Appendix D will 
have to be consulted for the purpose of correlating the symbols appearing in 
the flow chart with rhe r-register numbers appearing in the HPL code, which 
is contained in Appendix C.3. 
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Figure 28b. %TNIT Logic Flow (cont.) 
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Figure 28d. 25TNIT Logic Flow (cont.) 
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Figure 28e. %TNIT Logic Flow (cont.) 
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Figure 28f. STNIT Logic Flow (cont.) 


125 





INITIALISE. PAYLOAD ROTATIONAL STATE 




-*t<0 



Figure 28g. %TNIT Logic Flow (cont.) 
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5. ORBITER GEOMETRY LINK (%SSVU) 


The function of this link of the #TRAJ processor is to plot an end and/or 
a side view of the Shuttle Orbiter profile. It is executed only when the user 
specifies RSBY (Rectangular Shuttle Body-fixed coordinates) for the plot type, 
in the graphics data file (see Appendix E.3). In addition, it draws dashed 
lines that represent the nominal limits of visibility (through the overhead 
and aft windows) from the on-orbit pilot's control position. The code (Appen- 
dix C.4) is straightforward, for the most part consisting of HPL pit commands 
followed by the station coordinates (measured in inches) of the points that 
define the end and side profiles of the Orbiter. Typical Orbiter profile plots 
are illustrated in Appendix H. 
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6. FLIGHT SEGMENT INITIALIZATION LINK (%SNIT) 


The %SNIT link is called into memory and executed at the beginning of 
every flight profile segment. Each time it is executed, it reads a new flight 
segment definition (see Appendix F) from file number 1, to which the name "7*" 
has been assigned in the %TNIT link. 

After each new flight segment definition is read from disk, appropriate 
processing of the new input data is performed in preparation for state variable 
propagation through the ensuing flight segment. The propagation itself (in- 
volving numerical integration of the equations of motion) is performed by the 
%PR0P link (Section 7). Upon completion of the segment initialization process, 
execution control is passed to %PR0P by means of a chain instruction that 
causes the %SNIT code to be replaced by that of %PR0P. 

Trajectory integration is terminated (i.e., %PR0P is not called into 
memory) when an end-of-file mark is encountered during the attempt to read a 
new flight segment definition. If a payload solid rocket motor (SRM) burn has 
been simulated, the Particle Impact Damage Integrator Processor (#PIDI) is 
called into memory immediately after the termination of trajectory integration. 
Otherwise, the user is asked whether he wants to start a new run. If the 
answer is yes ((C ONTINU E)) , the Data Base Editor Processor (#DBED) is called 
into memory. If the answer is no, ((2), CC0NTINU|)) , program execution is 
terminated by a stp command. 

6.1 FUNCTION SUBPROGRAMS 

6.1.1 'KAM 1 

Given the user-supplied mnemonic symbol that designates an attitude- 
maintenance option, 'KAM' returns the appropriate numeric code for internal use. 
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6.1 .1 .1 Argument List 


None. The mnemonic symbol must be assigned to the string variable P$ 
before 'KAM' is executed. 

6.1.1 .2 Example of Usage 

If the character string "LVRH" resides in P$, the instruction 'KAM 1 -> K 
will cause K to be assigned the value 2. 

6. 1.1. 3 Computations 
See Figure 29. 
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6.1.2 1 JSEL* 


Given the mnemonic symbol that designates a Shuttle RCS jet-select table, 
'JSEL' returns the appropriate numeric code for internal use. 

6. 1.2.1 Argument List 

None. The mnemonic symbol must be assigned to the string variable P$ 
before 'JSEL' is executed. 

6. 1.2. 2 Example of Usage 

If the character string " P" resides in P$, the instruction 'JSEL' -> J 
will cause J to be assigned the value 1. 

6. 1.2. 3 Computations 
See Figure 30. 
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Figure 30 'JSEL' Function Subprogram Logic Flow 
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6.1.3 'COMP 1 


Given the mnemonic symbol that designates a Shuttle RCS cross-coupling 
compensation option, 'COMP' returns the appropriate numeric code for internal 
use. 


6. 1.4.1 Argument List 

None. The mnemonic symbol must be assigned to the string variable P$ be- 
fore 'COMP' is executed. 

6. 1.3. 2 Example of Usage 

If the character string "NONE" resides in P$, the instruction 'COMP' ■» C 
will cause C to be assigned the value zero (0). 

6. 1.3. 3 Computations 
See Figure 31. 
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6.1.4 ‘ KTOMS * 


Given the mnemonic symbol that designates which (if any) of the Shuttle 
OMS engines are to be fired firing a given flight profile segment, "KTOMS' 
returns the appropriate numeric code for internal use. 

6. 1 .4.1 Argument List 

None. The mnemonic symbol must be assigned to the string variable P$ 
before "KTOMS" is executed. 

6. 1.4. 2 Example of Usage 

If the character string " L+R" resides in P$, the instruction "KTOMS" ->- K 
would cause K to be assigned the value 3. 

6. 1.4. 3 Computations 
See Figure 32. 
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6.1.5 ' KTRCS ' 


Given the mnemonic symbol that designates a Shuttle RCS thrust command, 
'KTRCS' returns the appropriate numeric code for internal use. 

6.1 .5.1 Argument List 

None. The mnemonic symbol must be assigned to the string variable P$ 
before 'KTRCS' is executed. 

6. 1.5. 2 Example of Usage 

If the character string "+R0L" resides in P$, the instruction 'KTRCS' -> K 
would cause K to be assigned the value 7. 

6. 1.5. 3 Computations 
See Figure 33. 
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nction Subprogram Logic Flow 
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6.1.6 ' DATYP ' 


HFRMP input data files contain three general types of data: non-integer 

numbers, integers, and mnemonic symbols (character strings). Given the item 
number of an entry in a flight segment definition (see Appendix F), 'DATYP' 
returns a numeric code which identifies the type of data so that it can be 
printed in the appropriate format. 'DATYP' also loads the appropriate 
identification text for the data item into the string variable B$, and the 
description of its unit of measurement into the string variable U$. 

6. 1.6.1 Argument List 

pi = Item number in the flight segment definition (see Appendix F). 

6. 1.6. 2 Example of Usage 

The instruction 'DATYP ' (7) -> I would cause I to be assigned the value 1, 

the character string "SS XB RATE OR INCR " to be loaded 

into B$ , and the character string "DEG/SEC " to be loaded into U$. 

6. 1.6. 3 Computations 

The 'DATYP' code, which is straightforward and has no effect on trajectory 
computations, can be found in Appendix C.5. 
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6.2 SEGMENT DEFINITION LISTING SUBROUTINE ('SLIST') 

The 'SLIST' subroutine is used to list each flight profile segment defi- 
nition on the output line printer, immediately after it is read from the disk, 
in the format that is illustrated in Appendix F. The 'SLIST' code, which is 
straightforward and has no effect on trajectory computations, can be found in 
Appendix C.5. 
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6.3 MAIN LOGIC 


The %SNIT computations are described by the flow chart contained in Figures 
34a through 34f. The memory allocation table in Appendix D will have to be 
consulted to correlate the logical symbols appearing in the flow chart with the 
r-register numbers that appear in the HPL code, which is contained in Appendix 
C.5. 
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Figure 34a 


%SNIT Logic Flow 








m 



3 











Figure 34c. %SNIT Logic Flow (cont.) 
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%SNIT Logic Flow (cont.) 
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Figure 34e. %SNIT Logic Flow (cont.) 
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Figure 34f. %SNIT Logic Flow (cont.) 




7. STATE PROPAGATION LINK (%PROP) 


The function of the %PROP link is to propagate the state of the Shuttle/ 
payload system through a single flight profile segment. Digital output data 
describing the state of the system are printed at the beginning and the end of 
the segment, and at user-specified intervals between these two points. If the 
user has specified "RSBY" (Rectangular Shuttle Body-fixed coordinates) or "CPLV" 
(Curvilinear Payload Local Vertical coordinates) for the PLOT TYPE (in the 
graphics data file, Appendix F-3), then the appropriate graphical data are 
also plotted at each data output point. Upon reaching the end of the flight 
profile segment, %PR0P passes execution control back to %SNIT. 

7.1 ROTATED-ELLIPSE PLOTTING SUBROUTINE ( ' REL IP ' ) 

When the user selects the "RSBY" plot type, %PR0P draws one or two pictures 
of the payload (depending on how many views are specified) at each data output 
point. These pictures represent orthogonal projections of the payload's 
cylindrical outline into the X-Z and/or the Y-Z planes of the Shuttle's body 
axes. The pictures are composed of straight lines (corresponding to the sides 
of the cylinder) and ellipses and elliptical arcs (corresponding to canted 
views of the circles that represent the ends of the cylinder). 

The ellipses and elliptical arcs are drawn by using the HPL ofs instruction 
to move the origin of plotter coordinates to the appropriate point on the 
HP-9872A plotting surface, and then by calling ' RELIP ' to draw an ellipse (or 
portion thereof) about the offset origin. The geometry of the ellipse, refer- 
enced to the offset origin of plotter coordinates, is shown in Figure 35. 


148 



149 


7.1.1 Argument List 


pi = a = Semi -major axis of ellipse. 

p2 = b = Semi -mi nor axis of ellipse. 

p3 = Ep = Eccentric anomaly of first point on elliptical arc, 

p4 = aE = E l - Ep, where E L = eccentric anomaly of last point on elliptical 

arc. 

p5 = Number of chords to be drawn between the first and last points, for 
the purpose of approximating the true arc. 
p6 = cos e , where e = ellipse rotation angle, measured from Xp^Qj axis 
p7 = sin 0 to semi -major axis. 




7.1.2 Examples of Usage 


Suppose that it is necessary to draw the upper half of the ellipse shown 
in Figure 35, centered on the Shuttle body-fixed coordinates Xg = 20, Zg = -50. 
We will assume that the relationship between the HP-9872A plot coordinates 
^ X PL0T’ Y PL0T^ anc * Shuttle coordinates (Xg, Zg) has already been defined 
by the execution of an HPL scl instruction. We will further assume that the 
lengths of the semi-major and semi-minor axes, respectively, reside in regis- 
ters A and B. The values of cose and sine we assume to reside in registers C 
and S, respectively. The necessary instructions for drawing the desired arc 
are then pen; ofs 20, -50; cU * RELIP 1 (A,B,0,tt,N,C,S) , where N is the number of 
chords that are to be used to approximate the true shape of the semi -ell ipse. 

T he 'RELIP 1 subroutine can be used to draw figures other than elliptical 
arcs. For instance, the instruction pen; c ll 1 RELIP 1 (R,R,0,27r,40,l ,0) would 
cause a circle of radius R (approximated by 40 chords of equal length) to be 
drawn about the origin of coordinates. If this were followed by the instruc- 
tions pen; c ll 1 RELIP * ( R ,R ,n/2 ,2-rr ,3,1 ,0) an equilateral triangle would then be 
inscribed within the circle. 
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7.1.3 Computations 


The coordinates of any particular point P, on the ellipse shown in Figure 
35, are given by the equations 

Xrlot = a cosE cose - b sinE sine 
and 

Ypioi = a cosE sine + b sinE cose. 

As indicated by the flow chart in Figure 36, these equations are embedded in a 
loop and evaluated at regular intervals of E. The elliptical arc is approxi- 
mated by drawing straight lines (chords) between the points thus defined. 

Since a great number of chords may be required to obtain an accurate 
approximatior. of the true arc, the trigonometric identities 

sin (E + 6 E ) = sinE cos fiE + cosE sin <$E 

and 

cos (E + 6 E ) = cosE cos sE - sinE sin <5E 

are used within the loop. This avoids repetitive references to the sin and cos 
functions, which could result in excessive execution time requirements. 

It should be noted that the 'RELIP ' subroutine makes use of the volatile 
simple variables (H,I,J,K and W,X,Y,Z); therefore, values assigned to those 
registers by the calling routine will be lost upon execution of ‘RELIP 1 . 
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Figure 36. 1 RELIP 1 Subroutine Logic Flow 
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7.2 FOURTH-ORDER RUNGE-KUTTA INTEGRATION SUBROUTINE ( 'RK4' ) 

The function of the'RK4' subroutine is to propagate the state of the 
Shuttle/payload system across one integration time step, using the fourth-order 
method of Runge-Kutta. 

7.2.1 Input Data 

7.2.1 .1 Argument List 

pi = h = Value of time step. 

7. 2. 1.2 Others 

The ' RK4 ' routine calls ' DERI VS ' ; therefore, all the 'DERI VS' input data 
listed in Sections 2. 3. 1.2 and 2. 3. 1.3 must have been defined before calling 
' RK4 ' . 




3 A 


7.2.2 Output Data 


' RK4 * updates the contents of the extended array of state variables (r25- 
r74), which will be referred to symbolically in this section as [X]. It should 
be noted that time is one of the state variables in the HFRMP, and is integrat- 
ed just like any other state variable. It should also be noted that, although 
the contents of the derivatives array ( rl 25-rl 74) will change as a result of 
executing 'RK4', the values that reside there upon return from ’ RK4 ' do not 
represent the true derivatives at the end of the time step. To obtain the 
true derivatives, it is necessary to call 'DERI VS' again after executing 
' RK4 ' . 


7.2.3 Example of Usage 

The instruction c£l ' RK4 ' ( 300) would cause the state of the system to be 
advanced 300 seconds. 
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7.2.4 Computations 


Let [X] represent the array of first derivatives (with respect to time) 
of the state variables in the array [X]. Let [X] n represent the state of the 
system at time t n , and [X] n+ ^ represent the state at t n+ ^ , where 

h = Vl " *n 

is some relatively small time interval. Since [X] = f([X]), then according to 
the fourth-order method of Runge-Kutta we can write 

[X] n+1 - [X] n + (DO, + 2[k] z + 2[k] 3 + k[4])/6, 

where 

[k], = h f([x] n ) 

[k] 2 = h f([X] n + [k] ] ) 

[k] 3 = h f([X] n + % [k] 2 ) 
and 

[k] 4 = h f([X] n + [k] 3 ). 

The error introduced into the system state by a single fourth-order RK 
integration step is on the order of 

<< 5 m hi _ 

d t 5 5! 

As a rule of thumb, it has been found in using the HFRMP that the integration 
stepsize should always satisfy the relationship 

h < 40 degrees/o) 

HlaX 
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The symbol u> mav represents the angular velocity magnitude, measured In degrees 

IMWA 

per second, of whichever one of the four state- variable reference coordinate 
systems (B,b,G, or g) is rotating most rapidly with respect to inertial space. 
That Is to say, no reference coordinate system should ever be allowed to rotate 
more than 40 degrees during a single integration step. In some cases it may be 
necessary to reduce the single-step rotation limit to 20 or even 10 degrees to 
achieve the desired integration accuracy. 

A flow chart of the 1 RK4 ' subroutine is shown In Figure 37. The [Y] array 
is stored in registers rl75-r224, and the [Z] array is stored in registers 
r225-r274. Each quaternion in the [X] array (only) is normalized, each time 
that array is updated, by means of the computational sequence 

n 2 2 i 

V^o + 1l + ^2 + q 3 m 

V m + % 

q^m -> q 1 
q 2 /m + q 2 
q 3 /m -* q 3 . 

The quaternion normalization procedure is never applied to the [Y] and the [Z] 
arrays; to do so would introduce a systematic integration error. 
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Figure 37. 1 RK4 ' Subroutine Logic Flow 
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7.3 MAIN LOGIC 


The %PROP cor mutations are described by the flow chart contained in Figures 
38 a through 38h. The memory allocation table in Appendix D will have to be 
consulted to correlate the logical symbols appearing in the flow chart with 
the r-register numbers that appear in the HPL code, which is contained in 
Appendix C.6. 
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Figure 38a. %PR0P Logic Flow 
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Figure 38b. %PR0P Logic Flow (cone.) 
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Figure 38c. SSPROP Logic Flow (cont.) 
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Figure 38d. %PR0P Logic Flow (cone.) 
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Figure 38e. "PROP Logic Flow (cont.) 
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Figure 38g. %PR0P Logic Flow (cont.) 
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APPENDIX A: 


BASIC QUATERNION OPERATIONS 


A. 1 QUATERNION ALGEBRA 


Reference A-l defines a quaternion to be a mathematical quantity of the 

form 

» Q 0 + + jQ 2 + kQ 3 (A-l) 


A A 

where Q q , Q-j , Q 2 , and Q^ are real numbers and where the products of i, j, and 

A 

k are governed by the following conventions: 


A A A A A A 

ioi = joj = kok = -1 , (A-2) 

A A A A A 

jok = -koj = i, (A-3) 


A A A A A 

koi = -iok = j, 


(A-4) 


and 

A A A A A 

ioj = - jo 1 = k. (A-5) 

A A 

In addition to having the properties of imaginary numbers, the quantities i, j, 

A 

and k also have the properties of unit vectors that are aligned with orthongo- 
nal coordinate axes in a three dimensional space. 

The circular symbol (o) which is used to denote the quaternion product is 
adopted from the notation of Reference A-2. This symbol should not be confused 
with the dot (•) that is used to denote the scalar product of vector algebra, 
which is governed by the following conventions: 


A A A A A A 

i*i = j* j = k*k - 1 , 

(A-6) 

A A A A 

j*k = -k*j = 0, 

(A-/) 

k*i = -i *k - 0, 

(A-8) 


A-2 


and 


i-j = -j*i = 0. (A-9) 

It is worthwhile to observe also the similarities and differences between the 
quaternion product and the vector cross product, which is governed by the 
conventions 



i x i = j x j = k x k = 0 , 

(A- 10) 


A A AAA 

jxk=-kxj=i. 

(A-ll) 


k x i = -i x k = j , 

( A- 12 ) 

and 


A A A A A 

i x j = -j x i = k. 

(A- 13) 


A quaternion can be thought of as having a scalar part Qq, 

and a vector 

part 


Q = i Q i jQ 2 + kQ 2 . 

(A-14) 

Sometimes, then, it is convenient to express Equation (1) in its 

equivalent 

form 


Q = Q 0 + Q. 

(A-l 5) 


The sum of two quaternions is defined by 

W * (P 0 +Q 0 ) + + J‘(P 2 +Q 2 ) + MP 3 +Q 3 ), 

(A- 16) 


A- 3 


and their product by 



^°Q = ^0% " p l^l ‘ P 2^2 " P 3^ 

+ + P ^0 + P 2% ' P 3 V 

+ J ( p o^2 + P 2% + P 3^1 " P l^ 3 ) 

+ k(P 0 Q 3 + P 3 Q 0 + P } Q Z - PgQ, ) . (A- 17) 

Equation (A-17) results from applying the distributive law of algebra along 
with the conventions defined by Equations (A-2) through (A-5). It sliould be 
noted that, in general, quaternion multiplication is not commutative {i.e., 

QoP / PoQ). Except for the commutative property of multiplication, quaternions 
satisfy all the requirements for the definition of a field. 

The quaternion product of a scalar and a quaternion i_s commutative, and is 
given by 

SoQ = !foS = SQ q + iSQ ] + jSQ 2 + kSQ 3 . fA-18) 

which follows from (A-17) when the scalar is treated as a quaternion whose 
vector part is zero. In a similar vein, the product of a quaternion and a 
vector if formed by treating the vector as a quaternion whose scalar part is 
zero. This results in 

SoV = ( -Q-j V -j - Q 2 V 2 - Q 3 V 3 ) 

+i(Q 0 V 1 +Q 2 V 3 -Q3 V 2) 

+ J’(Qq V 2 + ^3 V 1 " W 

+ k(Q 0 V 3 + - QgV ) ( A- 19) 


A-4 


and 


VoQ = (-V-jQ-j - V 2 Q 2 - V 3 Q 3 ) 

+ 1 (V-|Q 0 + v 2 q 3 - v 3 q ? ) 

+ 2 % + ¥l - W 

+ k(V 3 Q 0 + V 1 Q 2 - VgQ,). (A-20) 

By examining Equation (A-17), it is seen that the quaternion product can 
be expressed in the form 

M = PqQ 0 - M + P c ^ + Q 0 F + P x (A-21) 

which leads to a relationship, 

QoF = FoQ - 2 P x Q, (A-22) 

that is sometimes useful. Equation (A-22) shows that quaternion multiplication 
is commutative whenever the vector parts of the two quaternions are parallel 
to each other. 

AAA 

The conjugate of the quaternion Q = Qg + iQ^ + jQ 2 + kQ^ is defined by 
Q = Qg - iQ-| " jQ 2 “ kQg- ( A— 23 ) 

The norm of a quaternion is defined by the product 

QoQ - = Qg + Q* + Q 2 + Q 3 • (A24) 
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A. 2 COORDINATE TRANSFORMATION 


Quaternions that have a norm of unity exhibit properties that make them 
very useful for transforming coordinates from one Cartesian system to another 
having a different orientation. Although the distinction is not made in some 
of the relevant literature, a quaternion of this special class is more pre- 
cisely known as a versor. In any event, it is well to remember that when the 
word "quaternion" is used in relation to coordinate transformation or the 
orientation of rigid bodies, almost always it refers to a quaternion whose 
norm is equal to one (1.0). A convention of using lower-case alphabetic 
symbols to designate versors (unit quaternions) has been adopted in this 
report. In other words, use of the symbology 

q = q 0 + iq-j + Jq 2 + kq 3 (A-25) 

implies that 

^0 + q l + q 2 + q 3 “ 1 * (A-26) 

According to one of Euler's theorems, any two Cartesian coordinate systems 
F and G that have a common origin can be brought into coincidence by rotating 
F through some angle a about a single fixed axis. That is to say, the angular 
displacement of any Cartesian system G with respect to another Cartesian 
system F can be described in terms of a rotation about a single fixed axis 
which is usually referred to as the Euler axis. The orientation of the Euler 
axis can be defined by a unit vector 

d = dp = dg = id-j + jd^ + kd 3 (A-27) 

with components d-j , d 2 ? d^ which not only satisfy the equation 


A-6 


1 


2 2 2 

d? + d 2 + d 3 K 


(A-28) 


but which also have identical values in F and G . Assuming that the angle of 
rotation is restricted to lie in the range 0 £ a < 2ir, the unit vector d can 
be taken to define both the orientation of the Euler axis and the sense 
(direction) of rotation according to the customary right-hand screw convention. 

It can be shown (e.g., see Reference A-2 and A-3) that any vector V" 
having the form 


V F = lV Fl + jV F2 + kV F3 


(A-29) 


in coordinate system F can be transformed to its equivalent form 


V G lV Gl + jV G2 + kV G3 


(A-30) 


in system G by use of the equation 


V G " q FG 0 V F 0 q FG* 


(A-31 ) 


where 


q~Pg = cos (‘a a) + d sin 0 2 a). 


(A-32) 


It should be noted that in this system of notation the unit vectors i, j, and k 
are not associated with a particular set of reference axes. Instead, they are 
understood to be aligned with the axes of whatever coordinate system is desig- 
nated by the alphabetic subscript attached to the vector symbol. In this 
connection, a vector symbol having no subscript represents an intrinsic physi- 
cal value that exists independently of the system of reference. For instance, 
the symbol V .?nght_represent the inertial velocity of one body with respect _ 
to another, while Vp represents just one of many possible quantifications of V: 
namely, that resulting from the projection of V onto the axes of coordinate 
system F. 

A-7 




The double subscript FG in Equation (A-32) identifies the versor as being 
the orientation versor which defines the angular displacement of coordinate 
system G with respect to system F 1 . 

Th? inverse displacement (of F with respect to G) is defined by 

q GF = ^FG = cos ^ ' d sin ^ “)* (A-33) 

and the inverse transformation of coordinates by 

^F = q GF 0 ^6 0 q GF' (A— 34) 

It is easily seen that the equations 

^G " q GF 0 ^F 0 q GF (A-35) 

and 

^F ~ q FG 0 ^G 0 q FG ( A- 36 ) 

are equivalent to (A-31) and (A-34). 

We now consider a third coordinate system H whose angular displacement 
relative to G is defined by the orientation versor 

q GH = cos ( l 2 3) + e sin (' 2 3 ). (A-37) 

Application of the coordinate transformation law yields 


1 An alternate rotation that produces the same result is defined by 

, A 

qp G = - cos (' 2 a) - d sin ( l 2 a), which represents a rotation through an angle 

of 2tt - a in the opposite (-d) direction. However, there is no need to con- 
sider the alternate rotation in the present discussion. 


A-8 


V H = q GH 0 V G 0 q GH 

" q GH 0 ^FG 0 Vp o q pG ) o q GH 

(^GH 0 q FG^ 0 0 (^pq 0 » (A— 38 ) 

whence it follows that 

q FH " q FG 0 q GH' ( A-39 ) 

The above result can be extended to any number of successive rotations. For 
instance, suppose that the body-fixed frame B of a spacecraft is displaced from 
an inertial reference frame I by rotating the spacecraft first through a pitch 
angle a about Y-axis, then through a yaw angle 4> about its body-fixed Z-axis, 
and finally through a roll angle <f> about its X-axis. The total displacement is 
defined by the quaternion product 

q IB = [cos e) + j sin e)] o [cos (' 2 ip) + k sin ( l 2 ^)] o 


[cos (J* <f>) + 1 sin 0 2 <p) , 


(A-40) 


which, after carrying out the indicated multiplications, reduces to 


q IB = ^ C e S S " S 0 S 4> V 


+ i(C fl C. S + S. S, C.) 

0 ip <p 0 \p ip 


+ j(S n c, c A + C Q S. SJ 

0 \p <p 0 1 p <p 


+ k(C S. - S fi C, S.) 

G ip <f> 0 ip (J) 


(A-41) 


A-9 




where C 0 = cos 0 2 e), S 0 = sin (% e), = cos (* 2 \p) , etc. 

The angles 0 , \p, i represent one of many Euler angle sets that can be 
used to define the orientation of one Cartesian system relative to another. 

Any relative orientation can be described by three Euler angles representing 
successive rotations about particular coordinates axes, taken in a specified 
sequence; the only restriction being that the second axis of rotation must not 
coincide with the first or the third. Reference A-4 contains a useful 
compendium of the relationships between versors, transformation matrices, and 
all of the possible Euler angle sets. 


* 

Not to be confused with the single angle of rotation about the Euler axis 
which produces the same result. 
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A. 3 TIME DERIVATIVE OF AN ORIENTATION VERSOR 


If the coordinate system B is rotating with angular velocity relative 
to system I, the time derivative of qj B is given by 

q IB = ’ 2 ^ IB 0 fi B (A-42) 

or, al ternati vely , by 

^IB = ’ 2 ^1 0 ^IB" (A-43) 

The consistency of Equations (A-42) and (A-43) with each other can be readily 
verified by substituting the relationship 

^B = ^IB 0 0 *4 1 b (A-44) 

into (A-42), which yields 

Tib “ ^ ^ib 0 ^IB^ 0 ^“I 0 ^IB (A— 45 ) 

One of the major advantages of using a versor to define the orientation 
of a rotating coordinate system lies in the fact that it has a finite deriva- 
tive at every possible orientation (assuming of course that the angular velocity 
fi is finite), thereby facilitating numerical integration of the differential 
equations that govern the rotational motion of the system. Such is not the 
case when Euler angles are used. No matter what rotation sequence is chosen, 
it is possible for the derivatives of two of the angles to approach infinity 
in the vicinty of certain critical orientations. 

The problem of infinite derivatives can be avoided by defining the 
orientation with direction cosines (i.e., the elements of a coordinate trans- 
formation matrix); however, this requires the integration of nine real 


A-l 1 


variables as compared to four in the case of a versor. The integration of 
direction cosines is further complicated by the necessity of maintaining the 
normality and orthogonality of the transformation matrix, as defined by six 
different equations among the direction cosines. In the case of versor inte 
gration, only a single ancillary condition of this nature is of concern: thi 

maintenance or normality as defined by Equation (A-26). 
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APPENDIX B: 


COORDINATE SYSTEMS 


B-l 




B. 1 SHUTTLE STATION (STRUCTURAL) COORDINATES 

Shuttle Station coordinates are referenced to a Cartesian system that is 
fixed to the Orbiter structure. The origin of the system lies in the plane of 
symmetry at a point 400 inches below the cargo bay centerline and 576 inches 
forward of the aft face of the forward bulkhead of the cargo bay. The X-axis 
is parallel to the cargo bay centerline and positive in the aft direction, the 
Y-axis is normal to the Orbiter's plane of symmetry and positive to starboard, 
and the Z-axis completes a right handed orthogonal system. The X coordinate is 
referred to as the station (STA), the Y coordinate is referred to as the buttock 
line (BL), and the Z coordinate is referred to as the water line (WL). This 
system is shown in Figure Bl. 




B.2 PAYLOAD STATION (STRUCTURAL) COORDINATES 


The Pay 1 co Station reference system is defined such that its coordinate 
axes are parallel to the Shuttle Station coordinate axes when the payload is 
stowed in the Orbiter payload bay. The origin of the Payload Station coordinate 
system is located in the center of the front face of the payload cylinder. This 
system is shown in Figure B2. Coordinates in this system are always measured in 
inches . 
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B.3 SHUTTLE BODY COORDINATES 


Shuttle Body coordinates are referenced to a Cartesian system that is 
fixed relative to the Orbiter structure, with its origin at the Orbiter CG. 

The X-axis is parallel to the cargo bay centerline and positive in the forward 
direction, the Y-axis is normal to the Orbiter's plane of symmetry and positive 
to starboard, and the Z-axis completes a right handed orthogonal system. These 
axes are parallel to the Shuttle Station axes; however, the positive directions 
of the X and Z axes are reversed. This system is shown in Figure B3. 


I 
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B.4 PAYLOAD BODY COORDINATES 


Payload Body coordinates are referenced to a Cartesian system that is fixed 
relative to the payload structure, with its origin at the payload CG. The X- 
axis is parallel to the payload longitudinal axis and positive in the Orbiter's 
forward dirction when the payload is stowed in the Orbiter's cargo bay. The Y- 
axis is normal to the Orbiter's plane of symmetry when the Payload is stowed in 
the cargo bay, and the Z-axis completes a right-handed orthogonal system. This 
system is shown in Figure B4. 
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Figure B4. Payload Body Coordinates 


B-9 


B.5 RECTANGULAR LOCAL VERTICAL COORDINATES 


Rectangular Local Vertical coordinates are referenced to a rotating 
Cartesian system centered at the CG of an orbiting vehicle. The direction of 
the Y-axis is opposite to the orbital angular momentum vector of the vehicle's 
CG with respect to the center of the earth, the Z-axis points toward the center 
of the earth, and the X-axis completes a right hand orthogonal system. This 
system is shown in Figure B5. 
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B.6 CURVILINEAR LOCAL VERTICAL COORDINATES 


Curvilinear Local Vertical coordinates are referenced to the CG of an 
orbiting vehicle. As shown in Figure B6, the X and Y components of position 
are measured along the surface of an imaginary earth-centered sphere that 
passes through the CG of the vehicle. The Z component is measured normal to 
the sphere, positive in the direction of the earth's center. These coordinates 
are essentially identical with those defined in Section B.5 when the distance 
from the vehicle CG is small. 
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ORBITAL MOTION 

= rY 
= r6 
= r - P 


Figure B6. Curvilinear Local Vertical Coordinates 



B.7 MEAN OF 1950.0 COORDINATES 


Mean of 1950.0 coordinates are referenced to a non-rotating earth-centerea 
Cartesian system that is defined by the orientation of the earth's mean equator 
at the beginning of the Besselian year 1950.0 (i.e., the Julian Date 
2433282.423357). The Z-axis points in the direction of the earth's angular 
momentum vector, the X-axis is aligned with the intersection of the equatorial 
and ecliptic planes (positive in the direction of the sun as seen from the 
earth at the time of the vernal equinox), and the Y-axis completes a right hand 
orothogonal system. This system is shown in Figure B7. 
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X-Y plane is Earth's equator of epoch 

X is directed toward the mean vernal equinox of 
epoch 

Z is directed along Earth's mean rotational axis 
of epoch and is positive north 

Y completes a right handed system 

Epoch is the beginning of Bessel ian year 1950 
for the Mean of 1950.0 system or the date of launch 
for the Mean of Launch system 


Figure B7. Mean of 1950.0 and Mean of Launch Date Coordinate 


B.8 MEAN OF LAUNCH DATE COORDINATES 


The Mean of Launch Date system of reference is also a non-rotating earth- 
centered Cartesian system defined in exactly the same manner as the Mean of 
1950.0 system, except that the orientation of the equatorial plane is fixed at 
the date of launch. This system is shown in Figure B7. 
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PROCESSOR CODE 
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C.l Base Link (#TRAJ) 


0* “#TRflJ<l223/10MOV79>"J »to "RUN" 

l: utb 6» 11 Hist #6 

2: *et ";:RMRT",0.1 

3« "RMG1 " *• 2 rtf rc<Pl/2rt>+p2 

4: if p 2<8J 2rt+p2*P2 

5* re t p2 


6: 


O • 

i”< • 


10 

11 

12 

13 

14 

15 


16 


17 

18 


19 


20 
2 1 


24 



26 


" RHG2 " : 2 rf f r c < p 1 / 2 rt > + p2 
if p2>if? p2-2rt+p2 
ret p2 

"flTNl " : 2rt+p5! ?t o +2 
"fiTM2 ,, »0*p5 
Plt2*p35 p2T2+p 4 
if p3+p 4=0? 9 to +6 
if p3>p4» , 9to +3 

atn<abs<Pl ■•'p2> > + p 35 if p 2<05 rt-p?*p3 
?t,o +2 


rt / 2 - o. t n < p 2 / o. ba (p 1 > ) +p3 


if pi < 0 * p 5 - p 3 1 p 3 
ret p3 

"HM.S" *• obs.'ipl )-*p2 
i n t, f p 2 / 6 0 ) * p 3 i int < p 3 / 6 0 ) + p 4 
1 0 0 p 4 + p 3 - 6 0 p4+ p 2 - 6 0 p 3 ) . 1 0 0 4 p 3 


ret p3s ?n ■' p 1 > 






"SE C S '' i a b * p 1 .) / 1 0 0 + p 2 
i nt ''' p2) -*p 3! 1 00f rc-<P2> +p4 
6 0 < 6 0 p 3 + i n t p 4 ) ) + 1 8 0 f r c '■ p 4>+P 3 
ret P3»9n<pl ) 

" DO TP" : ret cp 1 rp2+r (p 1 + 1 > r <p2+1 ) + r p 1 + 2) r 02+2 > 

"SXV" • p 1 rp2-*rp35 Pi r 02+1 ) + r (p3+1 > 5 p 1 r<'p2+2) + r , 'p3+2' 1 5 ret 
'"v'RDD" : 1 +p 4 > ?to +2 
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1 


K 


30 


31 


'.;b 


3-4 


3 b 



■3"? 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 


54 


#TRAJ Cont. 

"VSUB":-1+p4 

rpl+p4rp2+rp3 


r(pl+l >+p4r<P2+l ) + Kp3+1 > 


r<Pl+2)+P4r<p2+2>+r<p3+2) ? ret 

"CRSP": n'pl+1 > r <>2+2>-r<Pl +2) r<p2+l )+ p4 

r(pl+2>rp2-rpl r(p2+2>‘p5 

rpl f'T.'2+l ) - r ( p 1 + 1 ) rp2 >r<p3+2) 

p5+r(p3+ 1 > 5 p4-m'p 35 ret 

"ROT" ’• rp2*H! ?to +2 

" I ROT" : -rp2+H 

rpl 5 rCpl + l > + Y5 r < p1+2>+Z 

r ( p 2+ 1 ) ■» I > r C p 2+2 >+J? r <p2+3> +K 

'• HH + I I-JJ-KK ) :-:+2 ( < I J+HK ) Y+ (. I K-HJ > Z > + rr-'3 

■: HH- 1 1 +J.J-KK > Y+2 < < JK+H I >Z+< I J-HK >X ) + r < p3+ 1 ) 

r HH-I I - J.J+KK >2+2 < < I K+HJ > >!+ ( .JK-H I > V > + r < p 3+2 ) 5 re t 


"PHOT 


1+H50+W; rp 


2 + ! r < p 2 + 1 > / 2 + Y 5 r ( p 2 + 2 > / 2+2? 


CHQ": 1+H+W! no +3 


QHQC": l+Hi-l+Wisto +2 


"Q0XC!"8-1 + H5 1+W 

0 r < p 2+1 > -*X5 W r p 2 + 2 > + Y 5 N r ‘ p 2+ 3 > + Z J r p 2 + W 
H r C p 1 + 1 > + 1 5 H r < p 1 +2 > * J ‘ H r <. p 1 + 3 > +K5 rpl +H 
HN- 1 JY-K2 - i'p3 


H' .'+ 1 W+J2 -KY+ r ' p3 + 1 ) 

HY+ JW+KJ!- I Z + r 1 p 3+2 ) 

HZ H IJ+ 1 Y- r '■ p 3+3 ) i ret 

"0231 "s 1 -H) r*'Pl + l > + I? r <pl + 3) +K? 3 to +2 

" 021 3 " : - 1 Ml j - r < p 1 +3 > + 1 5 r < p 1 + 1 > +K 

rpl +H5 r ( p1 +2>+J 

hh-i I + JJ-KK+:-:; 2< jk-hi >+y 

' H T H2 ’ ' 2 '' I J + HI . > ? r ( KX + YY ) > + r T P‘2 + 1 ) 
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J 


no +5 




#TRAJ Cont. 


60 i 

61 s 

62* 


6 3 s 
64: 
65: 
66 : 
67: 


h, : 


b y : 

70: 

71 s 


? 4: 


’ 6 : 


SO: 

81: 


34: 


Sb: 


8 ? : 


•RTN2*<J+IiH+K>*XPfiTN2'<J-I.H-K>*Y 
*fiNC2*<X+Y>-»rp21 '0X02* <M<X-Y>)4r<P2+2>! rev 
"Q313": rPl->H!r(PlM> + Hr<Pl*2>-»J! r<Pl+3)+K 
2»fiTH2* < r< 1 1 ♦ JJ) » r<HH+KK> !> + r<p2+l> 

»flTH2»<K»H)->X; *FlTN2»<J. I>+Y 
’ flNG2*<X+Y>+rp25 »flHG2» (X-Y> + r<p2+2) 5 ret 
"231Q": ?to +2 

"213Q":-1+W 

rPl/2-»XJcos<X>-»Ilsin<X>->X 
r < p 1 + 1 ) /2+V 5 cos. < Y>+J! s in<Y>+Y 
r<pl+2>/2-»Zlcosf Z)+K»sin<Z>->Z 
IJK-WXYZ+rp2 
IJZ+WXYK+r<p2+2-M> 

X JK+W I YZ+ r < p2+2 > 

I YK-WXJZ+r < p 2+2+W> 5 ret 
" 3 1 30 " : r ( p 1 + 1 > /2 +Y ? cos <Y>+J5siri<Y) *Y 
( rPl+r^pl+2>>/2->XJcos<X>*I* s i n<X> +X 
< rpl-r<pl+2> > '''2-> Z 5 cos < Z > +K 5 s i n ( Z ) +Z 
J I •» rp2 5 YK* r < p 2+ 1 ) 5 YZ+r < p 2+2 > 5 JX+ r < p 2+3 > 1 ret 
" ItlflTQ” * r<Pl+5> + I}r<Pl*6>*J» r<Pl+l>*K 
r C p 1 +7 ) + X J r ( p 1 +2 > +Y * r<Pl+3)+Z 
1 + rp 1 + r >' p 1 +4 > + r f p 1 + 8) ->H5 1 f H<’ 15 3»t- o +2 

l/2rH+US 1- 4W*rp2 \ <1-X 'M*n’.p2+1> 5 < J-Y >W+r*'p2+2> 5 •,K-Z>U->r<p2+3) 5 ret 
2-H*W$W+2rPl*HS if HClisto +2 

1 • 2 r H ■* 1 J 5 < I - X > U *rp2J 1 •••" 4 W * r •' p 2 + 1 ) 5 ( K + Z > W + r ! p2+2 > » '■ J + Y > W * r < p 2+3 > < r e t 
W+2r''pl+4>*H* if H<l» 9t.o +2 

1/2TH+N5 <J-Y>W+rP2» <K+Z)W*r < p 2+1 > J l/4U+r (p 2+2> » <I+X>W+r<p2+3.*? ret 
W+2r ( p 1 +8 > *H 

1 /2rH+W5 (K-Z>W + rp25 < J+YM4+ r<p2+l>5 <1 +XMOr (p2 + 2> 5 1 ■ 4W*r<p2+3> 5 ret 
"DERI VS" : 25pl+0 
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*TRAJ Cont. 


901 ‘GRflVlcll ’QCXQ'CC^SiO+ZgtO+BrMcll 'QXQ' <0+87*0+80*0+129) 

9H 0+rl+r2l l+r3lcl 1 ' ROT’ < 1 .0+25.0+133) 1 -1 + r3l ell * ROT* < 1 .0+129. 1 > 

92: r <0+33) >RJ r24 RR>Clr<Q+133)+U*r<0+134)+T* r<0+135)+S 
9 3: . 0032431 <r23-'R) f 2+BJ C( B<5SS-- 1 )/2- 1 ) + R5 C6S + B 
94 : BU > r* 7 i BT »rSt B3-R+r9 


95: r (0+77' +L* r<0+78) >M* r <0+79> +Ni 3C- R+C 


93: 


Cr2r3(N H> + r 1 35 C r 1 r3<L 40*rl7*Cr 1 r2(H-L> *r 13 


97: r 1 24 r< 0+75) tP5 l f P=0! for 1 = 10 to 1510+rIlnext i;?to "flMGRC” 

93: "hTM" : . 00339342+E* TT +UU+C* -ESrC+D? EC+CJ B/< 1-C)+D 

99 : 2 03 5 5 5 9 1 > ! + < 1 + 3C • 4 ) C/2 ) +B 5 BD'R+fl f < R-B ) < 1 -DR/ 2 > >H 

1 00 : . 0O2 3 77 e xP ( r ( <. H -394 1 34 ) /' 1 2335 > + 1 7 . 449 ) -*Q 

10l: 7. 29 2 1 e -5*U! R < r i 0+35 > +TN ) + r 1 5 -RUW+ r 2 5 - r >: 0 + 34 > ■* r 3 

102: ,-.ll * ROT’ ( 1 * 0+87* 1 ) * i f pl=l*?to "Pt.fiERO" 

103 : "REPO" : r l rl + r 3 r 3 +R* f<R+r 2 i‘ 2 ) + V* r 2 /V+S*rfl V+C 
104 : O + Y'.l+XSif C #05 r 3 /VC--*VJ rl/VC+X 
105 : * RTN 2 * <S*C)-+B 5 ’ RTN 2 * < Y. X)-R 5 nbs<S) + S* s im 2 B> +E 
103 : < . 7 S 3 + 2. 41 3 ob= ( V)M . 316 ) < 1 -S> + ! . 928 S + . 346 abs<EYY> *D 
107 : - 1 345 PDQV+D 5 c 1 1 * SXV* <D » 1 * ID) 

103: . O 3 1 < ; l n > 3B) -Fees <2fl) >-. 0 36B+I 
109: i. . 04 5co*>. 3fl ) - . 1 32Y)CC- . 045Xc os( 36 ) + J 


1 10 : 
111: 
1 12: 
113: 
114: 


. 037 £ l iv 2 . 0936 > - . 0 1 < Y+Xr o b*. 1 E> >E+K 

POWERS 5321 Cn-B5 1 0493SR-+R5 RI * r 1 35 B J* r 1 4 > Ht - r 1 5 5 ?to +4 
" PL REPO " : r 1 22n bi v r 1 ) + r 1 2 3 r ( r2 r2+ r 3 r 2 i -R 
rpo-r: c 1 1 • 2. ; v * ■ R » i» lO) 

-rl21. E'-’r 1 i 0-*r2 + r'3! ell ’CRSP'd* 10* 13) 


115: "ENDEPO" sell ’CRSP* <0+84. 10. 1 > 


113: ell ' VERB* 113*1*13)1 ell * I POT’ < 10*0 + '37< 10 > 

1 1 7: 32. 174. n. 0+ 39 ) - G Jell * SX V * <G* 10* 10) 

113: ell 'POT*- 13*0+80. l>Jcll ’ VRDD’ >' 16* 1 * 13 1 

119: "RNCRC "sell ’POT' (0+33.0 + 30* 1 >* 1. r 1 •* r4 5 Mr 2 + r 55 Nr 2 * r3 
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#TRAJ Cont. 


120« ell ’ CRSP’ ( 1 * 4* 4 ) l c 1 1 ' VSUB* (16*4*1) 

S 21 > rl/L+rllr2/M+r2J r3/N>r35cll * IROT' < 1 *0+80*0+136) 

122S if Pl = l» ?to "PLTHR" 

123s “THR‘ , $ r94 + J5 Jnod2+I * <. J-I >.• 2+K 
1 24 s 1 9 . 1 bk » r 1 4 1* 1 9 . 1 6 1 + r 1 40 * 0 + r 139 
1 25s bOOOGU+K) *R!c 1 1 ' SUV* ui* 91 * 13) 

12b s r95 + J* i f J -0« for 1-142 to 1 475 0 + r 1 5 ne a I f -at- o +11 

127$ for 1=13 to 155 rl+RC 1-12* JJ+rllnext I 

123$ for 1 = 1 3b to 1 38 1 r I +RC I - 1 32* J ]•» r U next I 

129$ for 1 = 142 to 1 47* R[ I - 1 35* J 3+ r 1 5 next I5+to +8 

130$ "RLTUR": 1 *r 1 74 ** 0+V + r 1 3+ r 14 + r 1 5 + r 164 1 if r 1 1 9-0 J o +7 

131$ r 7 4 + T * if "PR[2 3*5t o +6 

132$ 1+V; T RC 1 3 + Ti if T\=Rl73$*to +3 

133$ fl[ 7 ]•*>;; 8C8 3-»Y*i read 3* RC 7 3*RC 8 3* if TR[73*?to +0 
134$ i fit 8 3 - v 9 *; Ftr. 7 3-:;>+flt 6 3 * v-:;rc 6 3 >rc 5 3 
135$ -kRC 5 3+TRC b 3 1 RC 3 3+rlb4+G* -32. 174GRC43 r64+G 
136$ for 1 = 13 to 1 5* Gr ( I + 1 03 ) + r I J next I 
137$ "ENDTHR" sell ’ 1 ROT' t 1 3*0+87 *13) 

138$ for 1=7 to 9 * r I + n. I +3 > + r 1 1 +6 ' + r ‘ I +0 + 1 23 ' * next I 
139$ "RMC " $ r 1 0+76 1 -fi* if R-0*?to "EMDDER" 

140$ it fl=l*for I --lb to 1 u; r * 1+0 + 1 20' + r I * ne .t I * *t o -t? 

1 4 1 $ 0 •* r 1 5 ■ r 1 0 + 3 5 1 • t 2 * ■■ r 1 0 + 1 3 1 : Rr2+r i 

142$ ell ’ 1P0T’ * 0 + 36 * 0+87- 7,' * r 1 1 ’ V3UB' 1 7- 1 < 4 • 5 c 1 1 ’ CF':.P' ■ 4 . 1 

143$ 2r 5-R! if p 1 = 1 * it V# 0 • R RC 6 3 rl64+rl64 rb4 • ■•R 

1 44$ V 1 1 ' CROP* iQ + 133* 1*4 1 * r > 0 + 34 '•*':■* r > 0+ 35 -'-W 

145$ . 00 3248 1 r24 1 r23. RR ) 1 2t r 5 r< 0+ 1 35 > + rS r 1 0+ 1 34 > ' -43 rS- R + R-*R 

1 4 6 $ *, r 1 0+ 1 30 1 -26N > F'+k* 0*r45 -o-r 55 R. RW-t 9 -P + l! W'*r3+rb 

147$ cl 1 •' C P'3F ' i 1 < 7 • 7 • » c 1 1 ’ VRDLi ’ 4 » 7* 4 '• *• c 1 1 ’ PUT ’ •'.4*0+87*4) 

14 : 3 : for 1-16 to 1 J: * r ( I +0 ► 1 20 * - r ( I - 1 2 > > r I 5 r*e t I 

149$ if Pl=!'+to "PI. tiMI" 
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150* " AM I" ! f o r 1 = 7 to 9!0+rI8next I 

151 ! for K=16 to 18?K~12+H?2H+J? J-l+G? if rK<05 G+J? J+l +G 

1 52* ob*<rK>*fliabs<:R[Hf J3>-»E! if R>E5 E+R5 1 +r285 

153s . 0 1 4 F 5 r99+C> r<H+92)+B{F+4r<RD>/E+TJCint<T/C>+T 

154: if T305T-F+S5O.B+Y+Z? no +2 

155! C + T? T F >85 EES$/ 1 6B~fl+B? B/cbs <R[ H> G 1 > +2? 2T/S+Y 

15b! i .A+BVE+XJ XT S + l-l 

157! for 1=7 to 95 I-6+H? r I +XRI H» J )+ZRt H> G ]+r 1 5 next 1 

153! for 1 = 136 to 138? I -1'3'2+H? r I +KRC H> J ]+ZR[ H> G 3+r I 5 next I 

159! for I ='142 to 1 47? I - 1 35+H? r I +WRC H> J 3+VRC H> G 3+ r I ? next I 

160! next t?cll ’ IROT' 1 7> 37» 7) 5 c 1 1 » VADD' < 130»7> 130) »*to +5 

161! " F'Lfif'II " : f or 1 = 161 to 163? r I -r <1- 1 45 ) * r I ? next I 

162.* cl 1 ’PUT’ <16- 105,- 1 >? -Lrl + rl ? - Mr2+r2? -Nr3 + r3 

163! ell » IROT’ > 1. 105. 1>? for 1=1 to 35rI/2+Y 

164: obs<Y )■>'!? M+V-r< 1+1643? X-Y*r< 1+16?) 5 next I 

165: "ENDDER" : t or 1 = 7 t o 95 r < 1+0+ 123) + r I ? next I 

166: r< O +34) + 5 + r < 0 + 1 3 3 ) 5 r 1 0 + 3 5 ) + W ? R M W - r 9 + r < 0+134) 

167: < r7-2SW ) R + r < 0+ 1 35 ) ? 0 + r 1 ? - W+ r 2 ? r 8 "RW+ r3 
1 68 : cl 1 ’ QDOT * 1 0+25 > 1 . 0 + 1 25 > ? c 1 1 » QD0T ’ < 0 +29 > 0+36 >0+1 29 > 5 re t 
169! "RUM" : ?et "".RMAT" > 169 
170: end 
-16018 
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C.2 Response Matrix Computation Link U'RMAT) 

0t "*RHAT<0403/24SEP?9)"!-no "RUN" 

l: wtb 6, 11 5 list #6 

2: ?et "UNIT'S 0.1 

3: "RUN" : lkd ? rod? if flUlJ^to +5 

4! dirt fiC33,BC 32 3,FC 12,123» IC3»33i JC8, 121 

5: dirt LC 3 1 , UC 12, 12 3* RC 1 2 « 1 2 1 , Tt 7 , 44 3 

6: dirt BTC 15 3, CTC 12,43, PTC 25 ] » Jf C 3,43, Kit 5 ] 

7: dirt MTC 44,3 3, PTC 3 3, V$C 80 3, NTC §0 3* sf * 11 
8: "* WARNING: SS " + B* 

9: " CONTROL IMPOSSIBLE WITH " + D* 

10: " V" + JTC13i" P"4.J$[ 2 ]f " PZI" + JTC3 3 

1 1 s " JETS"+K* 

12: frtt 1 - / < c 30 > / 

13: n i ?n "TCMIiID" , 1 5 s-read 1 , CT 

14: oi?n "TJFTM", 1 1 for 1 = 1 to 44 5 s. read 1»M*C I 31 next I 
15: o = =irt "TJFT", ISsread 1 , TC * 3 

IS: ostm "l*", 1 Jsread l»BC*35it BC 32 3=9 5 no "RMEND" 
17: for 1 = 1 to 3! BC I + 1 3+IC I , I 35 next I 
13: -B£5]+IC2,3]*It3i2] 

19: -BC83 + IC1.33 + IC3,13 
20: -BC 7 3+ IC 1 , 2 3+ 1 C 2 > 1 3 
21: inu I + I 

22: tor T=1 to 3*d3P "’.RMAT" , JTC T 15 JnP T 
23: files T J S V , * V < *■ V F: 5 ■? t o + 3 
24: files TJSF'j +F', *F'R, *F'F5 n.-o +2 
25: files TJSPZI . *F'ZI , *PIIR» *P2IF 
23: i no LI, P, F5 " " -VTSJT5 sreod 1 , JC * 3 
27: for C = 1 to 125 tor N = 1 to 8 
28: JCN,C3m!if M=05?to +13 
2 9 : for 1 = 1 to 6 « T C I > M ] + r 1 5 r i e x t I 
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30: for 1=4 to 65 < rl -BC 1+4 J>/12+rI 5 next I 
31: -r4*r45-r6+r65cl 1 » CRSP» <4, 1 , 4> 

32: for 1=4 to 65 r I + TC 7, N ]r< I-3>+rI 5 next I 
33: for I = i to 6! UC I , C 3+r I +UC I , C 35 next I 
34: H$[ M 3+Pf 

35: if P*C 1, 1 3="R"5 9 + I5 ? to +3 
36: if F'fC 1 » 1 ]="L" 5 3+15 3to +2 
37: 7+1 

38: if 06? I +3+ 1 
39 : 3. 1071+W5 if M >33 5 . 0923+W 
40: UC I,0 3+U+IJC I f C 3 
41: next N 

42: 32. 174/BC 1 3+flH or 1 = 1 t o 3 5 flUC I , C 3+IJC I , C ]J next I 
43: for 1=4 to 65 UC I , C 3+LC I -3 35 next I 
44: nut IL+fl 

45: for 1=4 to 65 flC I -3 3+UC I > C 35 next I 
46: if T = 1 5 i f 0<7i=ito +4 

47: int<<C+l>/2>*P54P-2C-l+S5if SUC P> C 3 > 0 ? no +3 
4 8 : B$ & C «C3 & D * Sc J fC T 3 Sc K t + V $ 5 w ft 6.1, V $5"" + V $ 

49.' if 065 B$&s'’fiTT"&D4 : ScJ4 ; C T 3S,K$+W* 

50: next CJsprt 2, UC * 3, W$5 ara U+R 
51: if W*#""in,o "REND" 

52: for 0=7 to 12! int < <0+1 >/2>+P 
53: for H=1 to 10 
54: for H=0 to 1 5 4+ , iF'+H)i v iod3+E 
55 : -RC E , C 3+X 5 2E+K 5 i f X >0 5 K- 1 +K 

56: X/UCE>K3+Y5for 1=1 to 1 25 RC 1 , 0 3+YUC I , K ]+RC 1,0 35 next I 
57: next H5 0+X 

53 : for H=0 to 154 + ( F'+H > n od3+E 5 X+RC E , 0 3t2 +X 
5 9 : n e x t Hi X ' R C P , 0 312+ X 5 l f X < 1 e - 3 5 ■? t o +2 
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zm 




60 
61 
62 
6 3 
6-4 
65 


66 


6? 


6 '? 


70 


71 


74 


i b 

i i 


80 

81 


84 


o t* 





SRMAT Cont. 

next Hi^to +2 

4P-2C-1+S!if SRCP»C]>0?9to +2 

B« ,, fiTT ,, mW«T]iK$+U*Jwrt 6.1»Wltaro U+Rl^to "REND" 

for H-0 to l»4*-tP + H)rt.jd3+ESO+RCE»C]5neAt H 

next C 

for C = 1 to 6 

for E ■ 4 t o 6 

-PCEf C3+KJ 2E>KS if XWK-l+K 

M RC E < K ]-*V; f or 1 = 1 t o 1 25 PC I > C 3+YRC I , K ]+RC I , C ]? next I 
0-*R[ E» C ]! next E 
it T= 1 5 ?t. o +3 

l n t < 1 L + 1 *F'5 4P-iiC”l *S» i f SRC P * C ] >0 5 ? t- o +2 

BIT " ROT COMP "t: Ctt C J8 Dtl J*C T lt<kt-»V*l wrt 6.1,'/* 
n e x t. C 

"PEND" : iprt- 3.RC*]»U*5if T=l»3to "NEXTT" 


a r a 

P-*F 



it 

5 3t o ' FEND" 

if \ 

/*#" " 

! 

o "hBliRTF" 

for 

C = 1 

t o 

6 5 l n t (. f C + 1 1 - 

for N= 1 

t 0 

10 

for H=0 

t. 0 

15 1 4 ' P+H ) riO'd 

-FEE* 17 ]- 

' ! * 

■E*i5it ;::0;i 


X F'C E* 1. ]■*'.'? t or 1 = 1 to 12 5FC I«o ]tVRC l-M-f-C I , t ]5ne,,t I 
f "l t t H 5 0 "* X 

t O' r H - Q t O' 1*1 * t F' + H ) i'i o d 3 ♦ fc » !•> -t F C E » i. J ! 2 X 

i ’> r ■■■ t H < X F C P * i. ] T X ■*'.7 5 if M 1 e ■ 8 ’ ? t o + 2 

ne .--.t. N 5 ?t o "HPuPTF" 

4 P - 2'0- 1 -*:• 5 1 1 F [ P « C ] •• 0 5 i t 0 ' ft B U R T F " 

f O' r H - 0 to 1 • 1 * 1 P + H ■' n o< d 8 *• F. 5 0 ■* F L E » L ] 5 n e : ■ t H 

fi o t C 
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90: for C=7 to 12 

91: for E=l to 3 

92: -FCE'C3+Xi2E*K; if X>0iK-l+K 

93: X/FC E» K 3->V 5 f or 1 = 1 to 121 F[ I » C 3+YFC I » K 3-»FC I » C 35 next 
94: 0>FIE»C3!next E 

95: int < <C+1 >/2>*P5 4P-2C-1+S5 if SFC P»C 3<=05 •eto "RBORTF" 
96: next C5no "FEND” 

97: "flBORTF": B* &“ FULL COHP"SiDfiJ$C T 3&K$*l<i$ 

93: wrt 6.1jW*5ara R*F 
99: "FEND":sprt 4»FC*3»W* 

100: "NEXTT" : next T 
101: 9+BC32 3 

102: o.i?n :# 1*“ > 1 ? 5.prt 1 > BC * 3 
103: "RMEND" : ?et ,, JjTNIT“» 169 
104: end 
*24706 
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■ inn 


1 ti 


C .3 Trajectory Initialization Link (STNIT) 


Os " V.TNI T( 1238/1 60CT?9> " ! ^t o "RUN" 
l: utb 6* US list #6 
2: *et "\SSVU">0 p 1 


o' • 

" JD " : int ( ( p2- 1 4 > 

12>-tp4 

•»: 

P3-3207 

5+ int < 1 4b 1 ( p 1 + 48O0 + p4>/4)-t 

5: 

p5+int < 

36?<p2-2-1 

2p4>/12)-p5 

S: 

ret p5- 

int (.3 int 

p1 + 4900+p4 ) "1 00> 

?! 

" PLOT VP 

" • 
• 


o • 
• 

if Pt=" 

" 5 ret 

0 

Q ; 

if P*=" 

NONE"? ret 

0 

10! 

if P» = 

" R'8BY " 5 ret 

1 

11! 

if PJ= 

"CPLV"? ret- 

2 

12! 

ret -9 



13! 

"UNIT" 

• 

• 


14! 

if F* = 

" FT"? ret 

1 

15! 

if PT = 

" f FT " ? ret 

1000 

It'! 

if PT= 

" NMI " 5 ret 

60 7b. 115 

1?! 

if Pf= 

M " ? ret 

1/. 3043 

18! 

it PT = 

" KM"? ret 

1 000'" . 3048 

19! 

ret- -9 



20: 

"MONTH 

" : 


21 : 

if PT = 

" JAN " ! re t 

1 

22: 

it FT = 

" FEE'."? ret 

7i 

23: 

if PT = 

MAP"! ret 

3 

24: 

it PT = 

APR"? ret. 

4 

•7 cj • 

if PT = 

" MAY"? ref 

5 

2b! 

it F'f = 

" JUN" ! ret. 

6 


if PT = 

" JUL " ? re t 

7 

28! 

if F'T = 

" AUG"! ret. 

8 

2 9 : 

if PT = 

SEP"? ret 

9 


ORIGINAL PAGE fg 

0F P 00 * quality 
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%TNIT Cont. 



30 * if P *=“ OCT"; ret 10 
31 : if Pf = " NOV"! ret 11 
32 : if Pf =" DEC"! ret 12 
33 : ret -9 

34 : ” C" : ret p 3 < p 2-1 >+ p !-1 

35 : " T ":« ftx<pl > p 2 >* p 3 ! ret p 3 < p 3-1 > •' 2 + mn < Pl » p 2>-1 

36 : " DIflG ": 0 * M+X 

37 : for J =1 to PlJt'or 1 = 1 to J 

38 : » T * < T » J ) + K ? p 2 + K + H 5 p 3 + » i + K J r H ■» r K + Y 

39 : if I = J ? W + YY + W * 1 1 r < p 4 +' C * ( J * J > p 1 ) ) 5 <u o +2 

40 : X+ 2 VY-»X 10 ^r<P 4 + 'C' < I » J» pi >) + r<p 4 + , C* < J > I » p 1 ) ) 

41 : next 1 5 next J 

42 : if p 0 < 5 1 r < W + X > / p 1 / 1 e 6 * p 5 

43 : rX / pl + p 650 + p 7 ^ p 8 


44: 

for J 

= t 0 

pi ; 

f 0 

r 1 = 1 

to 

J-l 


45: 

p 3 + ' T 

' c. 1 » J ) 

fZ 5 

if 

o. b s ( 

rZ> 

< = p 65 

3 1- 0 +12 

46: 

1+P8J 

p3+'T» 

(Ip 

I) 

+H5p3+'T 

* <J» J>+W 

47: 

rRrrW+X; 2rZ+Y? 

V / 

r<XX+ 

YY>+K! if 

x<o;-k+k 

48: 

K>T(2 

<l+r(i 

-KK 

> > 

>+YJ r 

<1- 

YY>+X 

5 rH+K 

49: 

r +2 

XYrZ+Y 

YrW 

+ r 

HJKYY 

-2X 

YrZ+X 

XrW+rW;e+rZ 

50: 

f o r K 

= 1 to 

p1 ! 

if 

K=I 5 

3 1 0 

+4 


51 : 

if K= 

J < 3 1 o 

+ 3 







52: p 3 + » T * < K • I > - H i p 3 + * T * < K , J > + N 


53: 

rH 

v Z 5 

XZ+YrW+rHI X 

rW-YZ-rW 


54: 

p4 + » C 

’ IjpD+H 

1p4+»C» <K 

» J»p1)+W 

55: 

rH 

■» z ; 

XZ+YrW*rH5 X 

rW-YZ+ rW 


56: 

ne 

xt 

K 




57: 

ne 

xt 

I ? next 

j ; if 

p 8=1 j 0 +p 

85*»,o +3 

58: 

if 

p 6 

<p55 ret 




59: 

p 6 

/pi 

+ p6 





I 





i 


C - 1 3 
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VTNIT Cont. 


60: p7+1+p 7? if p7<10OS*to -16 

61: fmt “DIRG ERROR =“ » elO. 2» 2/5 urt 6.PlPb!fmt 5 ret 
62: "RUN" : wtb 6* 115 if fl?lll?to +5 
63: dim fit 8 3. BC 32 ]. GC 1 2 ]. RC 12* 1 2 3. TC 7 ] 

64 : d i rt Rf C ? » 2 3 * Bf C 4 0 3 > Cf C 2 . 4 1 » Df C 3 * 4 3 f Ef C 5 » 4 3 > Ff C 4 i 4 1 

65: dim GfC 3> 4 3* HfC 8» 3 h IfC6» 3 3* -JfC6» 1 3> K'fC SO 3> LfC 4 3* HfC 6 3, NfC 6 ] 


66: dim 0fC6< 3 3 > F'fC 4 h OfC 6» 3 3> RfC 6> 3 3i SfC 6< 3 3-UfC 10 MlfC 30 3 1 ZfC 30 3 5 sf i 1 1 


67: 

"SS"-RfC 1 35 "PL 

"-RfC 2 35 "TP" 

->RfC3 35 "DX ,, +fl«4 35 ' 

'DY"*RfC53i "DZ"-RfC63 

bO • 

"DR" -flic?]; "DS 

" *RJ-C S ] S "DT" 

-RfC 9] 


b’? • 

"0M50"->CfC 1 ]; " 

IMLD"+CfC2 3 



70: 

" M50"-DfC 1 35 " 

SLV"-DfC2 35 

" PLV"-DfC3 3 


71: 

"R3BY"-EfC 1 35 ” 

PSLV"-EfC2 35 

“CSLV“+E« 3 3! "ESLV* 

‘ EfC 4 3 5 "SM50"-EfC 5 ] 

1 £ • 

"PP6'Y"-FfC 13?" 

RPLV"*FfC2 35 

"CPLV"*F*C3 35 " EPLV' 

' -Ff C 4 3 

“* *i • 

i _< • 

" SENS " -GfC 135" 

PCOH"-GTC 2 35 

"R I HP" -GfC 3 3 



74: "SMR"-HfC 1 35 "ECC"*HfC 2 35 " INC" -HfC 3 3 5 "RRN* 


-Mf C 4 3 


75: 

" RRG " 

-H-TC 5 3 ; ' TRR' , -HfC6 35 "HR "-HfC735 

"HP " - H f C 3 i 




76: 

"F'CH" 

- 1 f C 1 3i"YflW"-IfC235 "ROL"-IfC335 

"FUB"-IfC4 35 

"RYB" -I f ( . 5 35 

"RZB" 

-Ifl63 

i i • 

"X"- JfC 1 35 " V " - J f C 2 ] 5 "Z"* JfC3 ]i "R"-»JfC 4 35 "S" + JfC 5 

35 " T"- JfC 6 3 



*7* £» i 

1 v • 

"+X " 

-OfC 1 3 5 " + Y "-OfC2 3i" + Z "-OfC 3 35 

"-X " -Of C 4 3i 

" -Y "-OfC 5 35 

"_7 " 

-OfC 6 3 

7?: 

"rf: 

-Of [ 13 5" PUT " -Of C 2 ] 5 " RFU " -OfC 3 3 5 

"P.RM" -OfC 4 3? 

" RX " - 0 f C 5 3 5 

"DHL' 

- 0 f C 6 3 

30: 

" RFR " 

-F'fC 1 3? "F'LP" -F'fC 2 3i " RPR " -Rf [ 3 3 ! 

"RHP" -F'fC 4 3 i 

"RF: " - F'fC 5 3i 

" OMR ' 

-F'fC 6 3 

3 1 : 

"RFT " 

-3fC 1 3 j "FLT"-SfC ~i 3 5 "F'RT"-3fC 3 35 

"PRT" -3fC 4 3! 

"RT " * 3 f C 5 3 i 

" OMT' 

-3fC 6 3 


o = =>n 

"-ID". 1 ! i re i.i d l»Zf 5 "#TRRJ"-ZfC 1 

. 5 3 




3 '• 

b >= € P < 

e rit "EUTEF TRRJECTORV ID TEXT" > 

f f : if f 1 =3 1 3 ; 

Jt o + 13 



34 : 

30-le 

n ‘ ! f 1 - 1 ? if 005 9 t o +2 

- 




i“. cr « 

for H 

= 1 to h 5 1 f & " "-f.fi next H 





3 6 : 

f Ht 2 

> •: 3 0 >2 <e 3 0? ."i i.j r t 6 * Z f > f f J f pi t 





C* i • 

lid : 

t o d > 0- rZ 35 5 ti 1 30-Q 





'ri{“( * 

20325 

722-r73> 1 . 407646'^e 16-r24 





S 9 : 

n i 9n 

" 1 * " < 1 ■ = t e o d 1 » Bf - 3 i BC 1 3 - r39 
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%TNIT Cont. 


BC 2 3+r 1 * BC 3 3+r35 BC 4 3+r6l -BC 5 3+r5i -BC 6 3+r4! -BC 7 3+r2 

cl 1 *DIflC* <3» 1* 1 • 7>i rl + r77J r3-»r78j r6-M-791c.il ' IMflTQ* < 7* 30> 

< 1076.7-BC8 3>/12M-84 5 BC9 3/12-M-85!<3?5-b£ IB 3>/12+r86 
U076.7-BC 1 1 ])•• 12+r2?5SBC 12 3/1 2+ r2 76 i t3?5-BC 133) ■ 12+r277 
for 1 = 1 to 3 » QBC I + 1 3 3 • r I * next IJcll » 2 3 1 Q* < 1 * 278) 

BC 17 3 U3-M-99! f or 1=40 to -4 7 ? 0 -> r 1 5 next I 
o Jin " 2* " * 1 5 3 r •= o d 1 » BC * 3 5 BC 1 3-* r64 
BC2 3+rl*BC 3 3 »r 35 BC 4 3 + r65 ~BC 5 3-M-5? -BC 6 3+r45 -BC 7 3+r2 
cl 1 ’ DlftG’ C 3* 1* 1 *735 r Hr 102 5 r3+r 103 5 r6+r 104 5 cl 1 » I MftTQ ’ ( 7 . 1 05) 
-BC 8 3/ 1 2 M- 109 SBC 9 3/12+rl 105 -BC 103/12+rl 1 1 
-BC 11 3 1 2-* r 232 5 BC 12 3/12-M-283 5 -BC 13 3/12+r284 
BC 14 3 24 »R-» r 1 20 5 BC 15 3/12+L*rl21 
flF!R-»rl22S 2PL-rl23 
♦ o r 1=85 t o 7 0 5 0 •* r 1 5 n 6 \t I 

0 £?n "3-*-"*l*for 1 = 1 to 1 1 j 3 . r <= 0 d l*GCI35n-=xt I 
fts 1 GC 1 3 > ^P-T S ’ PLOT VP ’ -*GC 1 3 * if GC 1 3=0! »»• o +13 
fts <GC2]:wP4S »UNIT*-»GC2 3 

bsiipSent "PER D Y TO PLOT '? " » I * if fl?13=0??to +0 
PclrSwrt 7G5* " 1 p540* 878* 1 4778* 1OO30" 5 P*n# 25fxd 0 
1 c 1 0 < 1 4 < 0 * 9 * ('■' e t’i 5 p It 0 * - . 3 * 1 b 1 Ki 
GC 3 3-C < GC 3 ]-r GC 8 3- , Y 5 1 f GC5 3 = 0*=ito +4 
£ c 1 U • - 1 4 L < 0 • - 9 C ! o f £ - ij C 8 ] > - G C 9 ] * G C 8 3 -- C 1 1 C 5 3 1 X 
0 0. GC 10]' 0< GC 8 ]« GC 11 3? 0 ' 0* -GC 10 3-0* M+C 10* GC 1 1 3 

, 0 ... 0 * GC 1 0 3 • 0 * GC 9 3 * GC 11 3 5 ■ 0 0 * -GC 1 0 3 < 0 * V * GC 1 1 3 
it GC 4 3=0 5 ?t o +4 

£ c 1 0- — 1 4 C » 0* -9C5 of £ -C GC 5 3-GC 7 3* -GC 9 35 GC7 3-CGC 4 3+X 
o < U * GC 1 u 3 * 0 > GC 7 3~L 1 0 * GC 11 3* >.o.x 0 * -GC 1 0 3* 0 > n * GC 11 3 
'■ 0 0 • GC 1 0 3 • 0* IjC 9 3 * GT 1 1 3 * ■* av O* -GC 1 O 3 * O* V • GC 1 1 35 p*ri# 1 
<:• £ ? n " 4 " * 1 5 £ r 1 0 d 1 • B C * 3 
BC ! 3-V! t t r 1 BC 2 3J+P4-; •MOUTH’ -115 BC 3 3 + D 
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1201 

12ls 

122: 

123: 

124: 

125: 

1 ii'b • 

12 ?: 

12:3: 
12'?: 
1 30 : 
131: 
132: 
133: 
134: 
135: 
136: 
137: 
133: 
133: 
140: 
141 : 
142: 
143: 
144: 
145: 
146: 
147: 
143: 
149: 


<’ Ji)' <Y» 11» D > -24 33232 >/’36524. 22+T 
<2304. 943+ < . 302+. 018T>T)T+X 
<2004. 256- <. 426+. 042T>T)T+Y 


<2304 . ?43 + < 1 . 093+ . 0 1 3T ) T > T+2 
rf/643000 + C; n/2-CX+rl i CY+r25 -if/2-CZ+r3 
ell '31313’ U» 19) 

■' SECS ’ t B! 5 ])-»r?4 5 BC 6 1+r 124 

o; in "?*" « 15 on end 1 > "G02" 5 sre-od 1 > BC * ] 

"G02" : BC 5 3 + r75; B[ 18 l+rlOOi 0*r76+r94+r95*rlQl+rl 19 
asm "5* M » lisread l»B£*3Jfts <BC13>+P» 

6076. 115+T5 r-24 + G? if Pf= " QM50 " 5 TBC 2 3+ft 5 BC 3 ]+E 5 it o +2 
TBC2]+r23*flf tbc 3]+r23*P5 <fl-P)/<fl+P) +E5 <fi+P)/2+fl 
fi a -EE ) + P 5 QBC 7 ] +F S p/ < 1 + ECOS < F > > + R 5 r < G/'P > Es i n < F > + S 


QBC 5 3+r 1 5 QBC 4 ]+r2!QBt 6 ]+F+r3 


if Pf="0H5Q"; 
r 1 +H 5 r 2 I 5 r 3+ 
eos<U)*K5 s-in< 


ell ’ 313Q* <1*25) Jell >QCXQ* <19* 25»25>J it o +3 
U 5 c os < I +C 5 s i n < I ) +D 
LO +L5 KK-LL+M5 2KL 4 N 


1 . 62405e-3+ Ji Jr'23t2+Vi V/P+B 


R+BDDM/6 

S-BWDDH/ 


a\ 

<P A 

PS* 


v 


R 5 fl + B < 1 -3DD/2 > /3*fl 5 r < G/fi >/fi+U 
+ 85 fl+2fiflV< l-3BHLL>/3RRR+fl 
I J - B <1 - 7 H D / 6 ) N ■' 2 P * r 3 5 I + V C D M / 2 R R + r 2 i H + B CM / 4P+rl 
ell •' 3 1 3Q ’ f 1 < 25 ) 

if^2+r4» -r4+r5i 0 + rb> el 1 ’3130.’ <4'Q'*f cU ’01 10’ '25«0*25.' 

r i < 2-'R- 1 /fl ' j - s 3 ) ■■ ' R+r 3 55 3+r 3 4 ! R*r33 

for 1 = 1 t o 3 i OBC I +3 1+ r I .• next I 

ell ’ 2 3 1 Q ’ <1* 87 ) 5 ell ’ QXO » < 25 > S7 * 29 > 

f o r I = 36 t o 33 J QBC 1 -23 ]- r I i next I 

ell 'HER I VS’ •:0)5if/lS0+Q 

ft-s < BC 12 ];^F'f ; if Pf = " S'LV" i c l 1 ’ ROT’ 1 1 » 87» 4) J el 1 ’VflDD’ <4,36,36) 
o z i n " 6 * " , 1 5 s r e o. d 1 > BC * 3 
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%TNIT Cont 


150l ' for 1=7 to 12! BC 1-5 3-»r I ! next I 
15H fts <BC 1 ]>4P*i if P$="RSLV" ! ^to +4 
152s ell ’ ROT* < 1 » 37* 4 > 5 cl 1 » VSUB» <36»4»4> 

1 53 s ell ’ CRSP* ( 4 > 7* 4 > ! c 1 1 ' VflDD' <10»4» 10) 

1 54 s ell 'IR0T'<7*87*7)!cll * IROT’ < 10* 87* 10) 

155: r9-r33*r9! rl2-r34*rl2 
1 56* ell 'CRSP'(.l*7,4>»cll * V0DD» <10»4» 10> 

157s ell ’ I ROT* < 7* 25 * 7) 5 c 1 1 * IROT* <10»25» 10) 

153! r’DOTP’ <7*7>*R+r58Scll ’ CRSP» <7* 10» 4> 

159: r'DOTP' (4>4>*H! H/RR*W+ r60 

ISO: for 1=4 to 6! -r I/H+r I ! -r< I+3>/R+r< I+3> ! next I 
loll - ’ BOTP ’ < 1 0 * 7 > -> r 59 5 c 1 1 * CRSP' <4i 7> 1 > 

1 62s ell ’ IMflTQ* ( 1)50) 

1 S3: for 1 = 1 to 3! OBC I +3 ]*r I ! next Ilcll * 231Q’ < 1 » 1 12 

164! fts. <B[8])-*P*J if Pf = " PLV" ! no +2 

165: ell >QCXQ*<5O*29»0>lell ' QMQ' C0» 112* 112) 

166! ell ’ 0;;Q’ 1 50* 1 12*54) 

167: for 1=61 to 63 * OBC I -43 3-* r I ! ne xt I 
163: fts. (BC 12 ]> ’Pf* if P f = " M50" > no +6 
169: if Pf = " PLV " 5 n o +3 

170! ell ■' UL44U. ’ 1 39 > 5 4 * Q ) ! c 1 1 ' PUT ’ ( 36 * O * 4 ) 

171: ell ' VhL'D J ‘ 4.61 < 6 1 '■ * no +3 
172: ell » DERI VS* •' « 

17 3! ell ’ROT’fl* 112. 4) Sell ' VRUD ’ ' 4 * 6 1 » 6 1 > 

174! o.s9ri "7*"*1 

175: if GC 1 3 = 1? chain "‘■iSSVU" * 169 
176: chain "VBNIT"* 169 
177: end 
*23109 
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C.4 Orbiter Geometry Link (%SSVU) 

0: ,, 5iSSVU<1108/i6SEP?9>":Vto "RUN' 


l: 

utb 

6* m list #8 


2: 

get 

"*SNIT"»0, i 


o • 
• 

"SVM 

1" : obs( * nt (p6> )-H 

5 p 5/H"* J 5 c 

4: 

cos < p4 > +X 5 s i n< p4 > * Y » 

0*K 

5 s 

pit 

p!+X'p3»p2+Vp 3* if 

K=H » ret 

6: 

•JX+I 

Y*WJ IX-JY+XJW+YS 

K + 1 -> K j ■) t. 0 

"? • 

♦ • 

"RUN 

":deg» 1 2GC 2 3-> D ? DGC 3 1-»C 

8 s 

pen# 

2! if GC5]=05gto 

"SVU" 

8s 

sol 

0 * ~ 1 4 C j 0 . 9 C 5 o f s 

-DGCSl.EiGI 

1 9 i 

' 1 i rl 

DGC S3* DGC 8 1-CGC 5 ], -DG[ 9 j. 

1 1 s 

Of s 

-12r85* 1 2r8S~375 

125 

-us 


13s 

pen 

Spit 0 , 2 6 1 


14: 

pit 

350. 262 


15: 

Pit 

S 1 80 .263 


16: 

pit 

3468 » 300 


17: 

p 1 1 

3468 » 312 


18s 

pit 

3200 j 327 


19: 

p 1 1 

3 1 3 0 . 3 3 6 


28: 

pit 

31 25 » 403 


21: 

Pit 

3 1 50 - 399 


22: 

p 1 1 

3175. 399 


< *! S 

p 1 1 

3200,. 403 


ro 

4^ 

p 1 i 

3220.412 


25 • 

pH 

3241 . 424 


26 : 

Pit 

3232.424 


27: 

p 1 1 

3220. 422 


l2 !“t * 

pit 

3280 >413 


29: 

Pit : 

3175. 409 
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30: 

pit 

S150I 4|*9 "" SSVU C ° nt - 

31: 

pit 

5125*413 

32: 

pit 

8120*415 

33: 

pit 

81 30*4 30 

34: 

pit 

'3139* 450 

35: 

pit 

SI 43* 482 

38 : 

p 1 1 

8145*475 

37: 

pit 

8144*490 

38: 

pit 

8140*509 

3? s 

p 1 1 

SI 30*515 

40: 

p 1 1 

SI 15* 524 

41: 

p 1 1 

8100*528 

42: 

pit 

880* 528 

43: 

pit 

880.524 

44: 

p 1 1 

840.518 

45: 

p 1 1 

830.512 

48: 

p 1 t 

820.513 

47: 

pit 

8 4 * 8 1 8 

48: 

if 8 

- 0! 1 *85 ■sit o -35 

49: 

P 1 1 

-84,818 

50: 

1 1 fi€ 

2, 1 * 8 1 5 -* 2 • 4 34. 8 » V 5 2 4 C -* 2 

51 : 

PUT* 

pit ;■! ■ V ! p 1 1 ! 7ts2c o 5 1 52 ' * Y f-2i i n i. 52 1 

LT • 

p t n * 

pit ! 1 * V ! p 1 1 7 c i i = 1 1 1 . 5 ' ■ Y + 7 = i t'i ' 

Cj » 

"SVU 

" : i * 8[ 4 1 -O! jt ci "l.'tlli" 

54: 

= c 1 

0 « 1 4C < 0, 9C 5 o ♦ B CGI 5 ]+ [n,L 7 ] < ItijT 9 ] 

55: 

1 l m 

-D8C 7 ]* C IjC 4 3-IiGC 7 ], DGC 9 ]* r. GC 6 ]-DG 

58: 

of s 

1 2 r04 -1078. 7* 12r38-375 

S 7 j 

pen! pit- 1813*291 

tr • 

p 1 1 

1534* 293 

59: 

pit 

1 5 2 8 * 3 3 8 


1 1 n. 


[9] 


f. > ; 1 ins 


C - 1 9 


",SSVU Cent. 


60t 

pit 

1550*335 

611 

pit 

1575.334 

62« 

pit 

16O0.334 

6 3 s 

pit 

1623 

.3 36 

64: 

pit 

1689 

.412 

65» 

pit 

16O0 

.413 

66! 

pit 

1575.402 

6?s 

pit 

1550 

.392 

68! 

pit 

1525 

. 382 

6'?! 

pit 

1517 

» 376 

ro! 

pit 

1502 

• 441 

71 ! 

p i t 

1525 

-441 

3 • 

) kl • 

pit 

1550 

<442 

73! 

pi t 

1 575 

• 445 

74! 

pit 

1598 

453 

75! 

pit 

1 5 78 

• 526 

76! 

pit 

1550 

517 

77! 

pit 

1525 

506 

78! 

pit 

1500 

489 

79! 

pit 

1492 

483 

SO! 

pit 

1482 

519 

81 ! 

Pit 

1676 

5p2 

82: 

Pit 

1 700 

3 1 6 

C* s,* • 

Pit 

1 59 3 

8 1 6 

84! 

Pit 

1325- 

548 

85: 

Pit 

1 31 7. 

540 

86: 

Pit 

1 21 1 • 

5 3 3 

87: 

p 1 1 

1 808* 

524 

• 

p 1 1 

1 207. 

516 

'p • 

pH 

! : >.i 7 . 

420 
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90: 

Pit 

576* 420 

91: 

Pit 

576*500 

92: 

Pit 

500*500 

93: 

Pit 

491*499 

94: 

Pit 

480*493 

95: 

Pit 

472*489 

96: 

Pit 

434*450 

97 : 

Pit 

300*394 

98: 

pit 

283*385 

99: 

Pit 

270,373 

1 00: 

: Pit 

259* 372 

101 i 

' Pit 

250*364 

102: 

ell 

'SYM* < Si 

103: 

Pit 

263, 308 

104: 

pit 

286* 308 

105: 

Pit 

312*294 

106: 

p 1 1 

337*289 

107: 

pit 

362* 286 

108: 

Pit 

387* 283 

109: 

Pit. 

414*280 

HO: 

p 1 1 

450* 278 

111: 

pit 

4 9 2 * 275 

112: 

Pit. 

535*273 

113: 

Pit 

575*272 

114: 

p 1 1 

625* 271 

115: 

p 1 1 

675, 269 

116: 

Pit 

725*268 

117: 

Pit 

775, 267 

118: 

pit 

900* 265 

119: 

Pit 

1 00O *264 
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1201 pH 1100*263 
1211 pH 1150*262 
1 22 * pH 1250*261 


123 : 
124 : 
125 : 
126 : 
127 : 
1 23 : 
129: 
130 : 
131 : 
132 : 
133 : 
134 : 
135 : 
136 : 
137 : 
133 : 
139 : 


Pit 1300*261 
Pit 1350*264 
Pit 1400*267 
Pit 1450*271 
Pit 1500*275 
Pit 1550*231 
Pit 1600*238 
Pit 1613*290 

line' 2*1* 545. 2-»X* 434. 3*YS 24C+Z 
pen * p 1 1 . X * V * p 1 1 X+Zcos < 1 00 ) , y+Zs- x n< 1 00 > 
p e n > p 1 1 X * V * p 1 1 . X +Zc o s < 3 3 > , y +Zsin(38) 
560->X ? 480*Y 

p * n 5 p 1 x x ’ V * p 1 1 X + Z c o s < 3 0 > , Y + Z s. i n < 3 0 > 
Pen* Pit X , Y > p 1 1 X+395cos(-32) * Y+395sin<- 
E N D : p € n * p e n # 1 
chain "XSNIT" * 169 
e nd 


1 " 


3 


32)1 line 
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C.5 Flight Segment Initialization Link (%SNIT ) 


0: " , /.SNIT<1340/27OCT?9)"S'jto "RUN” 
1: wtb 6 . 1 1 ! 1 1 st #6 

2: ?et '".PROP" » 0» 1 
3: "KflM" s 


4: 

if 

P!=" D"J 

re t 

0 

Cj a 

if 

Pl=" IRH"; 

ret 

1 

rfi : 

if 

Pf-"LVRH" 5 

ret 

1— 

r • 

ret 

-9 



o * 

" .12 

EL": 



9: 

if 

M 

u_ 

ret 

0 

10: 

if 

Pf=" p" 

! ret 

1 

11: 

if 

Pf=" PZI" 

5 re t 

2 

12: 

r € 

t -9 



13: 

" C 

OMP": 



14: 

if 

PX = " 

5 ret 

0 

15: 

it 

P t ="M0ME" 

! ret 

0 

16: 

it 

Pf = " POT" 

5 ret 

1 

1 r : 

if 

PJ-="FULL" 

5 ret 

2 

IS: 

re 

t -9 



19: 

"t 

TONS": 



20: 

i t 

P -f = " 

5 re t 

0 

21 : 

i t 

F T - " HONE " 

: ret 

0 


it 

F‘ J = " L" 

! re t 

1 


it 

R|=" R" 

i re t 

o 

24: 

it 

Pf=" L+R" 

! ret 

3 

■ tr i 
C J * 

re 

t -9 



2 b*. 

"1 

TftCS" : 



I 1 • 

1 1 

PS = " 

i t e t 

0 

IS: 

it 

Ft--" DONE " 

'< r e t. 

0 

29: 

1 1 

p4-=" 

! ret 

1 


S»* 
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301 

if 

Pt*“ 

-X' 

J ret 

2 

311 

if 

P*»" 

♦y 

5 ret 

3 

32: 

if 

P*»" 

-Y* 

» ret 

4 

33: 

if 

P*»“ 

+2' 

p ret 

5 

34: 

if 

Pt = " 

-2' 

5 ret 

6 

35: 

if 

P* = " 

+R0L' 

! ret 

7 

36: 

if 

P*=" 

-ROL* 

5 ret 

3 

37: 

if 

P* = " 

+PCH' 

» ret 

9 

33: 

if 

P* = " 

-PCH* 

5 ret 

10 

39: 

if 

P*=" 

+ V R W 

J ret 

11 

43: 

if 

P* = " 

-YAW 

'» ret 

12 

41: 

re 

t -9 




42: 

"DATYP* 

: .ji'iP 

Pi 



43: 

"FLT 

SEGMENT TYPE <PR0P> 


" 

"+U»; ret 

44: 

"FLT 

SEGMENT LENGTH 


"HHMM.SS 

"+U#I ret 

45: 

"PRINT /'PLOT INTERVAL 


"HHMM.SS 

"-U#! ret 

46: 

"MAX 

1NTEG STEPS 1 2E 


" HHMM . SS 

"-UJfp ret 

47: 

" O' c 
o 

AERO FACTOR 


•• 

" +IJ! 5 ret 

CO 

" c*c 

■Z> o 

RATE OPT <INCR»IR»LVR> 


it 

"Hlf’p ret 

49: 

" C;C 

XB RATE OR I NCR 


" DEG-'SEC 

' +UI 5 ret 

53: 

" c c 

YB RATE OP I NCR 


"DEG/SEC 

ret 

51 : 

" C; g 

ZB RATE OP I NCR 


"DEG-'SEC 

" -pijf : ret 

52: 

" gg 

ATT MAINT OPT (Dp IRHp LVRH> 



" -tUf ! ret 

CO t 

•J * 

" 

XB DEADBAND 

"*B$J 

"DEG 

"*U*J ret 

54: 

»• c c 

•DO 

YB DEADBAND 


"DEG 

“^•Uf 5 ret 

55: 

" g C; 

ZB DEADBAND 


"DEG 

“ +UT p ret 

56: 

M C; C 

RCS OPT < V p P p PZ I > 


•• 

" -j U 4 5 ret 

57: 

"S3 

XC COMPENSAT I ON < NONE > ROT p FULL ) . 

"-»Bf 5 


ret 

53: 

“SS 

CMS CMD <NONEiLpRiL+R> 



"•»Uf 5 ret 

59: 

"S3 

PCS CMD •' NONE p +X < -X p+Yp...p + YAW « 

-YAW ) " •»BI' 5 

•• 

" •> U 1 ! ret 
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1 


\ 

I 

\ 

I 

K 


'f.'SNIT Cont. 

ret 1 
ret 3 
ret 1 

ret 1 
ret ! 
ret 3 
ret 3 

67« "SLIST": 

63! fnt 1.2/VFLT PROFILE SEGMENT" , f 3. 0 

6'?: fnt 2 * 2/ » "ITEM" » 2x» " DESCRIPTION" » 33x» "VALUE" , 2/ 


70! 

frit 

3, f 4 . O, c 42 

,f 14.4,cll 





71 ! 

f rlt 

4 • f 4 . 0 * c-42 

,f 1 4 . 0 » c. 1 1 





■J? “i • 

frit 

5 > f 4 . O , c 4 2 

, c 1 4 » c 1 1 





-vs • 
i z< • 

f (It 

9 > 2 " 5 f fit 






74! 

If 

pi M? i,i t b 61 

1 1 





“* C, « 
I >J * 

t.ift 

6 . 1 » p 1 * i,.i r t 

6 . 2 





?6i 

f o r 

1=1 to P2i 

jfiP ’ DATYP’ 

a> 




t c • 

i,i rt 

6.3. I < BT > BC I MJTi no 

+ 3 




78: 

i.i r t 

6.4, I , BT < B[ I ],UT; no 

+2 




79: 

1 1 ; 

(BC I 3>+PT! 

wrt 6.5- I,BT,PT,UI 




80: 

n€xt Ili'irt 6.9 

• ret 





8 1 : 

"PUN": rod! n IS 

0 ■» 0 





y d : 

TC 6 

] + 1 ■-* T C 6 1 « 0 

G C 12]« if TC 

6 3 = 1 ; 1 ->GC 12] 

" " H. T 



* 

24-» 

L ? o n e ri d 1 « 

"FIN" 





84: 

f >:■ r 

1 = 1 »o L : 

read 1 . BC I ] 

j next 1 5 c 1 1 

SLIST 

’ ( TC 6 

3 * L > 

y 5 • 

f o r 

I - 1 to 

SECS’ < BC 1 + 1 

] j f T C I ] * n e x t 

I! TC 1 

] + r 7 4 

-.TC 1 3 

86: 

6C5 

3-r75 






C; ^ • 

f t i 

■ BC 10 ]'.'-PT 

! ’ KAM'-Hn7 

6 




yy : 

f t s 

1 BC 1 4 ]> -*FT 

! ’ JSEL ' 





89 : 

l f 

J - o ; " V " - M T 

< j t o +3 






60! "PL AERO FACTOR "+BTS" 

61! "PL RATE OPT < INCR. IR.LVR) "-*6*1“ "+U# 

62! "PL XB RATE OR INCR ">BT! "DEG/SEC “+UT 

63! "PL VB RATE OR INCR "*BT! "DEG/SEC "-NJT 

64! "PL 28 RATE OR INCR "*Bt> "DEG/SEC "+UT 

65! "PL ATT MAINT OPT CD, IRH, LVRH> "-MJT 


66! "PL THR TABLE (NONE. IL> SA» SD) "*BT! " "HIT 


L 
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90: if J=t* "*P"->M$5*to +2 
91i • , *PZI ,, *M* 

92: fts <BC 15]>*P*J ’COMPOC 
93: if lOO 5 i v i'ix ( 1 ? C OC 
94: if J=0Jnin<l»C>-»C 
95: if C=0 5 Ilf 5 sto +3 

96: if C = 1 j -*Hf i no +2 

97: M«e H F“-»MI-»N* 

93: o =■ ?n Mf* 2isread 2>Rt*3»W* 

99: if Mf# " " 5 f lit ■'> c00> /'5 wrt 6« iJf-5 f nt 

lOO: fts <B[ 1710P4 5 ’ KTRCS* ->r95 

1011 fts (B[ 16 3 > -> F' -f i ’ KT0M3’ *r945 if r94=05 9to +6 

1 S2: < 1513-1076. 7>/12+rl 5 0*r25 < 492-375)/ 1 2+ r3 

103: if r94=l I 83/12-M-2 

104: if r94=25 -88 / 12->r2 

105: oil ’ VRDB ' < 1 ? 84 > 9 1 > 5 r ’ DGTP ’ < 9 1 » 9 1 ) *X 

106: for 1-91 to 93 i r I/OiT 5 next I 

107: for 1=96 to 98 5 QBE 1-85 ]*r 1 5 next I 

108: oil "DERIVS* 03> Stivl80-»Q 

109: 0* r ?•* r8 5 — r 3 3 r9 5 c J. 1 ’ CRSP’ < 1 1 7> 10) 

110: rl2-r34*rl25 cl 1 ’ ROT’ <10» 87* 10> 

111: ft s 1 ' B [ 6 3 ) * P $ 

112: if Pf = " I NCR” 5 for 1=4 to 65 0*r 1 5 next 1 5 n o +3 
113: it F‘f = " IR"5f or 1=4 to 65 r< I+32>*r 1 5 next- 15 no + 
114: ell ’ROT’ (l t 87 >4 >5 ell »VSUB* <36»4>4) 

115: for 1=4 to 6 5 QBC 1+3 3— r I r 1 5 next I 
116: o.= ?n Ml* 25 sread 2>R[*]»Wf 
117: if Wf#" " 5 f nt / j c80» /i wrt 6> Wf 5 frit 
118: for 04 to 65 if rO0 5n, o +6 
119: 200 5 if rf OiC-OC 
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120: nax<0»rK/RCK.C JW 

121: for 1 = 10 to 12! rl+TRt I-9»C 3+rI 1 next I 
122: for 1=36 to 385 r I+TRI I-32f C 3*rl 5 next I 
123: for 1=42 to 475 r I+TRC 1-35, C 3 + r 1 5 next I 
124: next K 

125: ell ' IR0T» (7, 25. 7 > » c. 1 1 » IROT* < 10, 29, 10) 

126: ell * CRSP* <7, 10, 4) 5 r33+R 
127: r’DOTP' <4 , 4) ->H5 H/RR-»N+ r35 

128: for 1=4 to 65 -r I/H+rl 5 -r< I+3)/R+r(I+3) Jnext I 
129: D0TP» <18, 7)->r345 cl 1 » CRSP' <4, 7. 1 ) 

130: ell » IMRTQ» (1 >25) 

131: os -an M*,25sreod 2,RC* J 

132: BC 1 8 3->rl 0Q5 ft s < B[ 23 ] ) +Pf 5 ’ KAM * ■» r 1 0 1 

133: fts kBC 24 3>*P* 

134: if Pf="N0NE" 5 a to +4 
135: if Pf = " "Sato +3 
136: if Lf=“ " 5 ato +4 

137: frit "SRM “*c4» " IGNITION COMMAND IGNORED" > - 5 wrt 6, P* 5 fin 
138: if rl 1 9#Q5 f' mt "SRM ",c4»" BURN CONTINUES" > / 5 wrt 6>Lf5frit 
139: ato +8 
140: pf*Lf 

141: tor 1 = 1 to 45 if P*Cr»I] = " "5 next I 
142: "00"S F'f [ 1,4 3-N* SN-fC 3»6 3-P* 

143: oi jn Nf , 35 = read 3, AC2 3, AC 'J 3« AC 4 3, AC 8 3 
144: f74-AC 1 35 r74+AC 2 3 =*RE 2 3 5 AC 8 ?.-AC 5 HO* AC 6 3-AC7 3 
145: 1 -» f 1 1 9 f r 1 1 6 5 0 -* r 1 1 7 •> r 1 1 8 
146: 0 3 an "8*" » 55 sr-rt 5, P* > K$ 

147: ell ’ DER I VS’ 1 1 > 5 tf / 1 80*i! 

148: 0*r7->r85 ~r58-r '95 c 1 1 ’ CROP’ < 1 > 7, 1 0) 

149: r 1 2 - r 5 9 •» r 1 2 5 c 1 1 ’ R 0 T ’ ( 1 6 , 1 1 2 , 1 0 ) 
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1501 ft* <BC 19 ]>*P$ 

151« if P$=" INCR" * f or 1=4 to 6S0+rI5next I5<»to +3 
152s if Pt = " I R “ 5 f o r 1=4 to 61 r< I+57>+rIlnext I5?to +2 
153: ell 'ROT' < 1* 1 12»4>5c 1 1 » VSUB* <61 » 4* 4> 

154: for 1=4 to 65 QBC 1 + 16 ]-r I*r I 5 next I 
155: for 1=61 to 635 rl + r< I-57> >r I J next I 
156: ell ’ROT’ <4, 105» 1) 

15?: for 1 = 1 to 35 rl r< 1 + 101 >*rl 5 next I 
153: oil »IR0T* <1» 105, 1> 

159: for 1=1 to 35 rI/2+YJ abs<Y)+X 

1 60s r( 1+64 )+X+Y*r < I +64 > 5 r ( I +6?) +X-Y+ r I +6? ) ? no xt I 

161: chain "‘;F'R0P" » 169 

162: "FIN": if rl 19=05 no +3 

163: fivit "SRM "»c4>" STILL BURNING AT TRAJECTORY TERMINATION" .2/5 wrt 6.L* 
1 64 * i p r t 5* "end" 

165: f fit 2 / » c y w » 2 ■ ' > c s 0 . 2 / 5 w r t 6. Zf. K$ 

166: for 1 = 1 to 35becp5uo.it 5005 next I 
16?: if L*#""5=iet "#PIDI" 

168: ent "START MEN RUN ? " , 1 5 i f fl 913=05 no +2 
169: 9*t "#DBED" 

1 70: stP 
1 ? 1 : end 
■*- ^ y y 1 9 
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8« M V.PROP<1103/16OCT79>":<jto "RUN" 

i: utb S.USlist #6 
2s <»et M *PIDI"»0»1 

3s "RELIP” s obs( int <p5> ) +H5 p4/H+ J5 cos<J>+ I ? sin< J>+J 

4s cos<p3>+X? sin<p3)+YS0+K 

5s pit Xp 1 p6-Yp2p7 > Xp 1 p7+Yp2p6 ? i f K=HJ ret 

6: JX+IY+WSIX-JY+XSW+YSK+l+KSeto -1 

7s "RK4" : p1+p25 p2/2+p3? p2/6+p4! 1 *p5 

8s for 1=25 to 745 r I + r < I + 1 50> + r < I +200) ? next I 

95 ell * DERI VS* <0> Sell * DERI VS* ( 1 ) 


10! 

f o r 

1=25 to 74! rC I + 150)+p3rn + 100>->rI5r.ext I 




1 IS 

for 

K=0 to 25 by 255 for J=25 to 29 by 4? J+K+HS 0+W 




12! 

for 

I=0 to 3? W+r t I+H>t2*NS next I?rW+U 




13! 

for 

1=0 t o 35 r ( I+H)/W* r < I +H> 5 next I S next JS next- K 




14! 

if 

p 5=4 1 ret 




15! 

for 

1=225 to 2745 rl+p4r<' 1-180 )*r I ? next I 




16! 

if 

p5=1 I p2' 3-'p 4? no +3 




17! 

if 

p5=2! p2*p 35 no +2 




18! 

p2/ 

6*p35for 1 = 175 to 224? r< I+50)s>rI ? next I 




19! 

p5+1*p 5? no -10 




20! 

"RUN": rud! 0+TC 7 3 




21! 

f fit 

1 ) 2 ) " 'j E T = " « t 1 0 . 3 ) " HHMM . S S " ) y 




22: 

t nt 

2«o 2' " : " • c 4 > t 9 . 1 ) c 4 ) f 8 . 1 > c 4 > t :3 . 2 < e 4 • t 8 . 2 < c 4 • t c 

. 2 - c 4 ) t 


2) c4 

23: 

t fit 

3 ) ■: 2 > " : "*c4» t9. 1 ) c-4 ) f 8 . 5 • c-4 > f 8 . 2 » c 4 « f 8 . 2 j c 4 > ♦ 6 

. 2 • c 4 ) f 


2 ) c 4 

24! 

f fit 

4 > c 2 ) " : " ) c 4 > t 9 . 2 ) c 4 ) t' 8 . 2 » 1 4 ) f 8 . 2 * c 4 ) f 8 . 3 ) c 4 ) f E 

.3< c4) f 

ft 

3 ) c 4 

25: 

f fit 

5 ) c 2 ) ": " ) c 4 ) f 1 0 . 2 ) c2 ) f 9 . 2 ) c2 ) f 9 . 2 ) o2 ) f 10.2»c3) 

f 1 0 . 2 . 

‘j> * 

f 1 0 . 2 - c 

26: 

f fit 

6 ) c2 ) " : " ) c 4 ) f 1 0 . 0 ) c 2 ) f 9 . 0 ) c2 ) f 9 . 0 ) c2 ) t 1 0 . 2 ) c 3 ) 

f 1 0 . 2 > c 

3 ? 

♦ 1 0 . 2 ) c 

•i* 1 • 

f fit 

7 ) c 2 ) " : " ) c 4 ) f 1 0 . 8 ) c 2 ) f 9 . 2 > c 2 ) f 9 . 2 ) c 2 ) f 1 0 . 2 ) c '3 ) 

+ 10. 3 ) c 

3) 

f 10. 3* c 

28: 

f fit. 

8 . c-2 « 0 5 ) f 9 . O > e 4 > f 8 . 0 ) c 4 ) f 3 . 0 > o4 ) f 8 . 0 ) c 4 ) f 8 . 0 ) c 

4 1 f 8 . 0 > 

c4 


29 : 

f fit 
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36 : 
31 * 
32 : 

33 : 

34 : 

35 : 

36 : 

37 : 

38 :' 


33 


40 

41 

42 


42 


44 

45 

46 

47 

48 


49 


50 

51 


54 


■jb 


59 


"PRINT": if TC73nocl2=0Jwtb 6*11 
wrt 6. 1. »HH.S»<r74> 
ell »DERIVS'<0)5 130/if4Q 

if r285#0» f "*** RCS CRN* T CONTROL SS RTT"*/*wrt 6ifnt *04r285 

for 1 = 1 to 3* r 1 4 r < I + 1 74 > > next I 

for 1=7 to 95 rI+r<I+178>Jnext I 

"SSTRRN" s r244G* r33*R5 r34*S; r35*W 

i:,R/<2G-RRRNW-R8S)4R;RR<2/R-l/fl-SS/G)4P 

Sr<P/G>4Y*P/R-l4XJ r<XX+YY>-»E! 'RTN1* <Y»X>*F 

fl/24 r 4 » - r-44 r5 ! 0-* r6 ! c 1 1 » 3 1 3fl ’ < 4 * 7 > ! c 1 1 » QXQC » < 25 * 7* 1 1 ) 

ell »QXQ’<19» llt7>lcll *G313»<7*4> 

Q ’ RNG 1 » < r 4 > 4H i Q r 5* 1 5 Q ’ RNG 1 • < r6~F ) 4 J 
6076. 1154T5QF4F 

wrt 6. 3* RfC 1 3*C$[ 1 3»R/T»H« 1 ]* E* Hit 2 3*1* Hf-C 3 3 » N » HTC 4 3 j J j HtC 5 ]> F> HIT 6 ] 

1 . 62405e-3+ j; Jr23T2-»V5 V/P*B 

oil ’ 03 1 3 ’ ( 1 1 > 4 ) 5 cos < r5 ) -*C j s i n < r5 ) -*D 

cos ( r6 ) *K ! s i n < r6 > 4L * KK-LL-*t1 * 2KL4N 

R-2RRVa-3riDLL>/3RRR4RJr<G/fl>/R4W 

fl-B < 1 -3DD/2 > / 34fi 5 R-BDDM/'6-»R j S+BWDDN/34S 

RR <2/R- 1 /fl-SS/G ) *P i Sr< P/G > 4 V » P/R- 1 ->X 

r(XX+YY)4E5 > RTNl ’ <YiX>+F 

r 6 + B (1 - 7 D D ■ 6 > N/2P4IJ5 r5- V C It M ■ ' 2 R R -> 1 5 r 4 - B C N ■ 4 P 4 H 
Q 1 4 1 j Q> RNG1 ’ (H)4N* O’ RNG1 ’ ( U - F ) - . J * Q F 4 F 
Cfl ( 1 -E ) - r23 ) /T4P i < R ( 1 +E > - r23 > /T4fl 

wrt 6.2*flfC 1 3* C$C 2 ]* fl.« H$[ 7 3* P* H$C 8 3* I * Hf[ 3 ] * N> HIT 4 ]* J< Hit 5 3* F* Hit 6 3 
wrt 6 

" 38R0T "sell ’ O; !Q ’ r 1 9 * 29 . 7 ) • c 1 1 »G231» <7*4 ) 

Q r 4 4 1 5 Qr54 J* u’ AN G 1 * ( r 6 ) 4 1 Q r 3 6 4 K 5 Q r 3 7 4 V * 0 r 3 8 4 Z 

wrt 6.4* flit 1 3*Dlt 1 ]» I* lit 1 3* J. 1 * r 2 3* K* I ft 3 ]> X » ITT 4 3* Y* lit 5 3* Z* IIC63 
ell ' ROT ’ < 1 * 87* 1 > • c 1 1 » VSIJB’ <36* 1. 7'> >cl 1 »Q231» <87*43 
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%PROP Cont 


60« Qr4-»I}Qr5->JJQ , flHGl'<r6>+K5Qr7-»XIQr8-»YIQr9-»Z 

61i wrt 6. 4* flit 1 ]i Dlt 2 3* I » lit 1 3*J* Wt2 3»K* lit 3 3* X* III 4 3* Y* lit 5 3* Z* lit 6] 
62: wrt 6 


63: "PLR0T":cl 1 » DER IVS* <1 ) 5 180/rf+Q 
64: for 1 = 1 to 35 rl + r< 1 + 177) 5 next I 
65: for 1=7 to 95 r 1 + r< 1 + 181 ) 5 next I 
66: ell ’ QXQ* < 1 9. 54* 7) 5 cl 1 » Q231 * <7» 4> 

67: Qr4 + I!Qr5+J!Q !, flMGi* <r6)+K5Qr61+X5Qr62+Y5Qr63+Z 

63: wrt 6. 4 .Plft 2 3. Dlt 1 1* I > lit 1 ]* J. lit 2 3* K» lit 3 ]. X. lit 4 3. Y. lit 5 ]. Z> I ft 6 3 
69: ell * ROT' <1 f 112> l>Icll • VSUB* <61 * 1 • 7) 5 cl 1 » Q231 * < 1 12» 4) 

70: Qr4 + I > Qr5+ J* Q* fiNGl * < r6)+K! Qr7+X5 Qr3+Y. Qr9+Z 

71« wrt 6. 4. flit 2 3. Bit 3 3. 1. lit 1 3. J. lit 2 3. K. I It 3 2. X. I It4 3. Y. I It 5 3. 2. I It 6 3 
72: wrt 6 

73: ell *QCXQ* <50*25*0)5 if GC13#2!no +9 
74: if GC 123=0! l+Gt 1 2 35-no +3 


< *J • 

76*. 

77: 

73 : 

79 

SO 

81 


84 i 
85: 
86 : 


“CPLV" : GC 2 3Gt 3 3 + 85 GC 2 3GC 9 3+T 

ell 1 Q213’ <0. 4) 5 -r4r53+r75 r5r53+r8! r5S-r33+r9 

for L=3 t o 15 if Gt 5-L 3=E 5 no +4 

sc 1 0 » — 1 4 S > 0 » -98 5 GC 2 3GC 8-L 3+X+M! if L = 1 5 W+SGt 5 3+W 

of s — i'l » - T 5 1 in X.X~SGt5-L3.T.T-SGC6 3 

ofs r <8-L> . r 9 5 Pit O. 0? pen? cl 1 ’ R E L I P ’ <$/20* $'20* 0* 2rf* 30* 1*0) 
pen! next L 

0+ r 7 + r 8 5 - r33+r9 5 c 1 1 * CRSP » <175*7* 10) 5 r 12-r34+r 1 2 
0 + r 1 3+r 145 -r58+r 15! cl 1 *CRSP* <173* 13* 16)5 rlS-r59+rl8 
ell * ROT* < 1 3* 0. 1 3) 5 cl 1 ’ ROT ’<16*0*16) 

c 1 1 ’ VSUB’ (.13* 7* 13) 5 cl 1 » VSUB* <16. 10. 16)5 if r 1 1 9=05 no + 1 2 
"°LSM50": cl 1 ’ROT’ <188.0. 1 S3 > 5 c 1 1 ’VSUB’ ''188.185. 188) 
c 1 1 ’ IROT’ <1 16.54* 1 > 5 cl 1 * I ROT’ <1.19.1) 
sprt 5. r74-flt 1 3. r 1 64 * r 1 . r2. r3 
c 1 1 ’ IPOT’ 18.25* l)!cl 1 * IROT* <1. 19. 1) 
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%PROP Cont 


901 ell ' IR0T'<16.25»4>5ell » IROT' <4» 19» 4> 

91* sprt 5» rl » r2i r3i r4i r3* r6 

92* wrt 6.6.fl$C2 3»EfC53, rli J*C 1 3. r2» J«21ir3» J4C3 3. r4»fi«41i r5ifl«5 J» r6»fl#C6 1 
93* ell »IR0T*<188»25>l>»cll ' IROT' <1 » 19 » 1 > 

94* sprt 5 » r 1 » r 2» r 3 


95* if r74<flC 2 3! no +2 
96* sprt 5> "end" 5 0+r 1 19 

97* "PLTRflN" : f or 1=1 to 35 r< 1 + 1 2>*r 1 5 next I 
98: ell ' CRSP' < 1 75 » 1 , 4 > ! c 1 1 ' VSIJB’ <16* 4> 4> 

99: wrt 6. 6i fifC 2 3» ETC 2 3> r 1 > JTC 1 3> r2i JTC 2 3* r3> JIC 3 3» r4> fl*C 4 3>r5» R*C 5 3i r6>fi$C6 3 
ICO: cl 1 ’ IROT’ <36.87. 7> Jell ' VSIJB* < 1 75. 7i 7) 

101: ell ' CRSP’ <7> 1 > 7> ! c 1 1 ’ VRDD' <4. 7> 4> 

102: ell ' ROT’ <! 1 > 3?» 1 > 5 el 1 ’ ROT' <4 , 87. 4) 

103: wrt 6» 6.< FlfC 2 3* ETC 1 3* r 1 > Jit 1 3» r2* J#C 2 3> r3> JTC 3 ]> r4 > RfC 4 3* r5» RTC 5 3> rS> fiTC 6 3 
104: wrt 6 


105: ell * ClCJiQ * (.54 < 29 > 1 3 1 > ! if GC13#l!no +15 
106: if GC 12 ]=0i 1 + GC 1 2 35 n o +14 

107: "R.SBY":GC2 3GC3 3+SSGC2 3GC9 3+TS r 120+R! -r 121 *r 105 0+r 1 l*r 12 
108: ell ’ ROT ’ '.'189* 1 8 1 > 7 > 5 e 1 1 ’ ROT’ < 10. 1 31 * 10) 

109: for L=0 to 15 if GC 5-L 3=05 no +10 


110 : 
111: 
112 : 
113: 
114: 
115: 
1 16: 
117 : 
118 : 
119: 


sc 1 9 f - 1 43 > 0 j -98 5 GC 2 3GC 8-L 3+X+W 5 i f L= 1 i W+8GC 5 3+W 
of s -W t -TS 1 i n '/, < K-SGC 5-L 3 » T » T-SGC 6 3 

of i r < 2-L) > r-3 5 p 1 1 0 > O 5 pe n Sell ' PEL I P * < £•• 20 j S-'20» ti •' 2 > 2fl » 3> 1 > 0 > 5 pe r> 
r< 1 1-L)+U5 rl2+V 5 -r (. 8-L >+115 -r9+Y 
if ( 1 -2L) r 1 10+L ) > 6 5 X + U + X 5 Y + V + Y 5 -IJ+U5 - V + V 
r(UU+VV)+E5fimax<0* 1 - < E r 1 2 1 ) T 2 > + B 5 if E=05 1-»C5 0*D5 no +2 
V/E+C5-U/E+D Ofy 

of s X > Y Sell 'PELIP' <fi» B» 0* 2if? 30> C» D> 5 if E = 05 no +2 ^ 

of e IJ > V Sell ’ PELIP’ <fl> B»0> it* 15) C> ID S iPlt -IJ,-V 
pen? next L 
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''■PROP Cont. 


"TPSENS'tcll » VSUB* <282. 109. 7> J ell * ROT 1 <7. 181 * 7> 

ell ’R0T’(61il81'13)lcll ' VSUB* <13*36* 13)lcll *CRSP' < 13. 7. 10) 

ell ’VflDD' (4. 1 0 > 4 > I c 1 1 ' VRDD* < 1 • 7» 1 > 

ell ’ VRDD’ ( 1 > S4 * 1 ) 5 c 1 1 ’VSUB’< 1.275. 1) 

ell ’ROT'a.ZYS. DJcll ’ROT' <4.278.4) 

r 1>X» r2 » V# r'3+2. r4nll r 5 »BI r6+C 

NK+VY»«lr<Q + 22)->R5ra->Q 

»flTNl'<Y»X)-»Sl»flTN2»<2»Q>+T 

C8N+BY+CZ) R+U 

if Q R< 1 e -t*» 0 »V > r (Rfl ♦■BB+CC-ULO/'R + Nl *t o +2 
tBX-RYVQQ+VJ <CR-UZ)rQR+W 
1 SO-' n «Q. QS*Si QT -»T» QV*V J QM+M 

wrt 6. 7. flf[ 3 3.G*t 1 3. R. Jft 4 3. S» Jft 5 3. T. Jft 6 3. U. Rft 7 3* V* Rft 8 3. N, Rft 9 ] 
wrt 6 

"R I PC ON"! for l-l to 6 * r < 1 + S4 ) ♦ r 1 5 no \t I 

wrt b. 8> FlfC 2 ]>Gft 3 3. r 1 . Oft 1 3. r2> Oft 2 3. r3. Oft 3 3. r4. Oft 4 3. i 5. Oft 5 3. rO. Oft t. 3 
w r t 6 

for 1=8 to 7* i t 1 +40>» r I » tie >;t I 

r3+r4-Ri R+r2«B 

re + r7 -*C i C + r5-D 

r2 + r 5- 1 • r Or 1 1 ?-» J > r 4 t r 7 ’K 

0 t- C ■* ! i > B *■ D ♦ Y » r 0 + r 1 -» Z 

wrt b.S*Hfr 1 3> Oft 2 3. rs.f'ft 1 3. r3» C'«T 2 3« r4-i'fl' 3 3.R. i'll 4 3. fc. Oft ‘S 3. rO.Oft 0 3 

wrt b . 0 « Oft 1 1- Oft 21. rS.Pft 1 1* rt'.Pft 2 3- rf.Rfl 3 3>i. < F-fC 4 3. Ii.Rft t>3» rl.Rft 3 

Mil b.e.ftft 1 1« Oft 2 3.1. Sft 1 3* J.Sft 2].| - Oft i: Ml- 3ft 4 l.Y.Sft 5 3.2. 3ft b 3 


if o b 3 . 174 -11 13" . OO01 . chciiri " , .3MlTMt"'? 

"RDVRtlb": Tt 2 3'. 1 I ret t*74 Tt 2 3> >-D5 if > 10 . UU01 i Tt 2 3-D 
D + r 7 4 - E S H Tt 1 3-E . OOOl.Tt 1 3-E!E-r74-D 


ir.t ■ . Tl 3 3-> + l »TL5 3!D. Tf *5 3 *TC 4 3 
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‘vPROP Cont. 


1501 if TC 5 3= 1 ; if absCTC 1 ]-T[43-r?4X.00O15Tt 1 ]-r?4 + TC 4 1 
15H fxd 2! dsp ":.PROP T = " i ’ HM. S’ < r?4 ) » c 1 1 * RK4 * < TC 4 3> 
1521 TC51-1+TC5]; if Tt 5 3=01 TC 7 3M+TC 7 31 dsp ""iilo "PRINT 
153: ?to -3 
154: end 
*290? 
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APPENDIX D: 
MEMORY ALLOCATION 


D-l 



14076469 x 10 = gravitational parameter 
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V 


SS STATE VARIABLES 


DESCRIPTION/NOTES 


SS local-vertical system ' 


orientation versor, with respect 


to MLD ECI system. 




SS CG distance from center of earth. 


dR/dt 


SS Keplerian angular velocity magnitude. 


t SS body ang ula r ve locity with 
respect to inertial space. 


.resolved onto SS body axes. 


SS gross weight (constant) 


QMS propellant consumption 


! W PRRX 


lb 

* 

forward 

lb 

aft left 

i lb 

t 

aft right 

i 1b 

forward 

! lb 

aft left 

l lb 

aft right 


control 



RCS propellant 


consumption 



control 


(RESERVED) 


(RESERVED) 












































PL STATE VARIABLES 


REG. 

SYMBOL 

UNITS 

DESCRIPTION/NOTES 

50 

q 0 


% 


/ PL local vertical system ' 


51 

q l 



* q_ l9 - 

orientation versor, with 


52 

q 2 




( respect to MLD ECI system 


53 

q 3 



J 

54 


q 0 



/ PL body orientation versor, . \ 

55 


q l 


- q Ib 

with respect to MLD ECI 


56 


q 2 



\ system. / 


57 


q 3 

.. 


58 

r 

ft 

PL CG distance from center of earth. 

59 

• 

r 

ft sec" 1 

dr/dt 

60 

“k 

sec" 1 

PL Keplerian angular veolcity magnitude. 

61 

“1 

sec" 1 



/ PL body angular velocity with ' 


62 

- co ^ 

sec" 1 


. -Ib 
“h 

respect to inertial space. 


63 


sec" 1 



( resolved onto PL body axes. 


64 

w 

lb 

PL gross weight 

1 

65 


lb ft se< 


Integrals of control torque applied 

67 

/ a 2+ 

lb ft se< 


in the positive direction about each 

68 

n 3+ 

lb ft se< 


body axis. 

69 

/£-j _ 

lb ft se< 


Integrals of the control torgue magnitude 

70 

fl 2- 

lb ft se< 

applied in the negative direction about 

71 

n 3- 

lb ft se< 

J 

each body axis. 

71 



(RESERVED) 

72 



(RESERVED) 

73 



(RESERVED) 

74 ! t (sec 

— X— * — - 

State time 
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SYMBOL UNITS 



PL CONSTANTS 


DESCRIPTION/NOTES 


PL aero factor 


PL att maint flag (0:D, 1:IRH, 2:LVRH) 


PL principal moments of inertia. 


PL principal axes of inertia 


orientation versor, with respect 



Position of PL actual CG with 


respect to "nominal CG" (front face 


of cylinder), in PL body axes system. 


value computed in 'DERI VS') 


U[!j (PL thrust unit vector, in PL body axes system. 


PL thrust flag (o:no, l:yes) 


PL radius 


PL length 


u d*V 4 


Atmospheric den si ty factor (applies to SS & P L] 


D-6 


























































-8 
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REG. SYMBOL UNITS 


MISCELLANEOUS DATA 


DESCRIPTION/NOTES 


/ SS sensor position with 


respect to nominal CG, in 


\ SS body axes system. j 


/ SS sensor orientation ) 


versor, with respect 


\ to SS body axes. j 




/ PL target point position with respect^ 



* pfit 

Lh_. 

to "nominal CG" , front fact of PL 

- 



cylinder), in PL body axes system. / 


Shuttle attitude control failure flag 


D-13 



















APPENDIX E 

USER-SUPPLIED INPUT DATA BASE 


E-l 


E.l SHUTTLE DATA FILE 


SHUTTLE DATA FILE * 


ITEM DESCRIPTION 


VALUE 


1 HEIGHT 

2 IXX 

3 IVY 

4 IZZ 

5 IYZ 

6 IZX 

7 IXY ... 

8 CG STA 

9 CG BL 

19 CG WL 

11 SENSOR STA 

12 SENSOR BL 

13 SENSOR ML 

14 SENSOR AXES PITCH 

15 SENSOR AXES VAN.. 

1 6 SENSOR AXES ROLL. 

17 DAP CYCLE 


200017. 0000 

LB 

887302. 0000 

SLUG • 

0386877. 0000 

SLUG- 

00 94 30 7. OOOO 

SU.JG-- 

-971 . 0000 

SLUG- 

247376.0000 

SLUG - 

5622. 0000 

SLUG- 

1095.3000 

IN 

O. 3000 

IN 

377. 4OO0 

IN 

1095. 3000 

IN 

0. 3000 

IN 

377. 4000 

IN 

O . 0000 

DEG 

O. 0000 

DEG 

0 , 0000 

DEG 

SO. OOOO 

MILL I 


FTT2 

FTt 

f r t 

FT t 
FTt 
FTt 


To i\j To To To 



E.2 PAYLOAD DATA FILE 


PAYLOAD DATA FILE * 


ITEM DESCRIPTION 


VALUE 


1 HEIGHT 

2 I XX 

3 IYV 

4 IZZ 

5 IYZ 

6 IZX 

7 IXY 

8 CG ST A 

9 CG BL 

19 CG ML 

11 TARGETED I NT ST A 

12 TARGETPOINT BL. 

13 TARGETPOINT ML. 

14 DIAMETER 

15 LENGTH 


37253. DODD 

7439. 0000 

45463. 0000 

45433. 0000 
0. 0000 
O . 0000 
0 . 0000 

300 . 6300 
0. 0000 
0. 2000 
300 . 6300 
O . 0000 
0 . 2000 
112 . 0000 

428 . 0000 


LB 

SLUG -FT I 
SLUG -FT t 
SLUG-FTt 
SLUG -FT t 
SLUG-FTT 
SLUG -FT f 
IN 
IN 
IN 
IN 
IN 
IN 
IN 
IN 
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r>j r-j ro ro ro r*j 


E.3 GRAPHICS DATA FILE 


GRAPHICS DATA FILE $ 


ITEM DESCRIPTION VALUE 


1 PLOT TYPE ( NONE * RSBY > CF'LV > RSBY 

2 PLOT UNIT> PU <FT» KFT> NMI > M> KM> FT 

3 PLOT SCALE. 120 PU/IN 

4 X PANEL WIDTH 14.0000 IN 

5 Y PANEL WIDTH 0.0000 IN 

6 PANEL HEIGHT 9.0000 IN 

7 X MAX 1500.0000 PU 

8 V MAX 0.0000 PU 

9 Z MAX 30.0000 PU 

10 TIC INTERVAL 20 PU 

11 LABEL INTERVAL 5 TICS 
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E.4 TIME/ATM DATA FILE 


TIHE/flTM DATA FILE $ 

ITEM DESCRIPTION VALUE 


1 LAUNCH YEAR 1931 

2 LAUNCH MONTH < JAN * FEB » MAR? NOV» DEC) FEB 

3 LAUNCH DAY 27 

4 LAUNCH GMT 1935.0000 HHMM.S 

5 INITIAL GET 1005. OOO0 HHMM.S 

6 ATMOSPHERIC DENSITY FACTOR 1.0000 


CO CO 


E.5 SHUTTLE I STATE DATA FILE 


SHUTTLE ISTflTE $ 

ITEM DESCRIPTION VALUE 


1 STATE TYPE <0(150* IMLD) OH50 

2 SMA 359S. 1726 MM I 

3 ECC 9.0012 

4 INC 28.3014 DEG 

5 RAN -89.5443 DEG 

6 ARG. -29.4152 DEG 

7 TRA 3.3863 DEG 

3 ATT REF <SLV> SLV 

9 PITCH -100.9000 DEG 

10 YAH -33.9000 DEG 

1 1 ROLL -160.O000 DEG 

12 RATE REF <M50»SL¥>,« M50 

13 XB RATE 0.0000 DEG- "SEC 

14 YB RATE O.0O00 DEG- SEC 

15 ZB RATE 0.0O00 DEG- SEC 
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E.6 PAYLOAD ISTATE DATA FILE 


PAYLOAD ISTATE $ 

ITEM DESCRIPTION VALUE 


1 STATE TYPE tPSLVi PS8Y> RS8Y 

2 X t.3775 FT 

A Y -0.0250 FT 

4 7 3.1292 FT 

5 XDOT 0.7120 FT/ SEC 

6 YDOT 0.0O00 FT/'SEC 

7 ZDOT -0.3392 FT - SEC 

3 ATT REF (PLWSBYT SBY 

9 PITCH 53.0000 TEC 

10 YAM 0.0000 DEG 

11 PALL 0.0000 DEG 

12 RATE REF <. M50 » PLV • SBY > M50 

13 :-iB PATE 0.0000 DEG SEC 

14 YB RATE 0.0000 DEG SEC 

15 ZB RATE 0.0000 DEG- SEC 
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APPENDIX F 

USER-SUPPLIED FLIGHT PROFILE DEFINITION 


F-l 


F.l FLIGHT PROFILE SEGMENT 


FLT PROFILE SEGMENT 20 


ITEM DESCRIPTION VALUE 


1 FLT SEGMENT TYPE (PROP) PROP 

2 FLT SEGMENT LENGTH 2. 4OO0 HHMH.SS 

3 PRINT/PLOT INTERVAL O.O5O0 HHMH.SS 

4 MAX INTEG STEPS I ZE 10.0000 HHMH.SS 

5 SS AERO FACTOR 1.0000 

b SS RATE OPT CINCRj IRjLVR) IR 

7 SS XB RATE OR INCR 0.0000 DEG'" SEC 

3 SS YB RATE OR INCR 0.0000 DEG/SEC 

9 SS ZB RATE OR INCR 0.OOO0 DEG -SEC 

10 SS ATT MAI NT OPT <D»IRH*LVRH> IRH 

11 SS XB DEADBAND 2.OOO0 DEG 

12 SS YB DEADBAND 2.0000 DEG 

13 SS ZB DEADBAND 2.0000 DEG 

14 SS RCS OPT ( V > P * P Z I > P 

15 SS XC COMPENSATION (NONE* ROT » FULL) ROT 

1 6 SS OHS CMB (NONE * L > R* L + R) 

1 7 SS RCS CMD ( NONE * +X , -X * (-YAW > -YAW > 

13 PL AERO FACTOR 1.0000 

19 PL RATE OPT CINCRp IRpLVR) IR 

20 PL XB RATE OR INCR 0.0000 DEG- SEC 

21 PL YB RATE OR INCR 0.0000 DEG- SEC 

22 PL ZB RATE OR INCR 0.0000 DEG-'SEC 

23 PL ATT MAI NT OPT (DpIPHpLVRH) IRH 

24 PL THR TABLE (NONE* ILj SR« SD) IL 
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APPENDIX G 
DIGITAL OUTPUT DATA 


G-l 




G.l DIGITAL OUTPUT DATA 


GET * 1056.350 HHMM.SS 


SSSOM50 

3617.3 

SMfl 

0.O0407 

ECC 

28.31 INC 

279.22 

RAN 

26.O0 

ARC 

153.42 

TRA 

S3: IMLD 

187.5 

HA 

149.8 

HP 

28.46 INC 

279.57 

RAN 

19.90 

ARG 

159.58 

TRA 

S3: M50 

-21.13 

PCH 

-57.48 

YAW 

93.13 ROL 

-0.000 

RX8 

0.000 

RYB 

0. 000 

PZB 

S3: 3LV 

-ii4.ro 

F'CH 

0.00 

YAW 

130.00 ROL 

0.000 

RXB 

-0.065 

RYB 

0. 800 

RZB 

PL: M50 

161.78 

PCH 

-12.72 

YAW 

37.86 ROL 

0.000 

RXB 

0.000 

RYB 

0.000 

RZB 

PL: PLV 

-3.51 

PCH 

-10. 14 

YAW 

0.63 ROL 

-0.011 

RXB 

0.068 

RYB 

0.013 

PZB 

PL : f M50 

-325994 

X 

-277410 

Y - 

•196164 Z 

-868.89 

DX 

-716.93 

DY 

-347.56 

DZ 

PL:RSLV 

413200 

X 

-992 

Y 

225774 Z 

1273.34 

DX 

-155.27 

DY 

106.44 

DZ 

PL: RSBY 

32454 

X 

1811 

Y 

469739 Z 

181.56 

DX 

155.31 

DY 

1164.21 

DZ 

TP: SENS 

47OS60 

R 

1.73 

S 

36.05 T 

1168.77 

DR 

0.263 

DS 

-0.003 

DT 

PL PIMP 

173 

+X 

205 

+ Y 

272 +Z 

173 

-X 

331 

-Y 

261 

-Z 

S3 PC ON 

62 

RFX 

0 

RLX 

0 R'RX 

0 

RAX 

62 

RX 

671 

Of'lL 

S3 PC ON 

46 

RFR 

30 

RLR 

31 RRR 

61 

RAR 

107 

RR 

671 

OMR 

S3 PC ON 

108 

RFT 

30 

RLT 

31 RRT 

61 

RAT 

169 

RT 

1341 

OMT 


GET = 10 

56. 

40O HHMM. 

, S3 













-i* 

i OM50 

36 

17.3 

SMA 

0. 00407 

ECC 

28. 

31 

INC 

279 

.22 

RAN 

25 

. 90 

ARC 

153.85 

TRA 

S3 

: IMLD 

1 

37.5 

HA 

149.8 

HP 

28. 

46 

INC 

279 

.57 

RAN 

19 

.90 

ARG 

159.91 

TRA 

S3 

: M50 

-2 

1.13 

PCH 

-57.48 

YAW 

98. 

13 

ROL 

— 0 . 

000 

RXB 

0. 

000 

RYB 

0. 000 

RZB 

C; C; 

: SLV 

-114. 37 

PCH 

0.00 

YAW 

188. 

00 

ROL 

0. 

000 

RXB 

-0. 

065 

RYB 

8. 000 

PZB 

F'L 

: M58 

16 

1.78 

PCH 

-12.72 

YAW 

87. 

86 

ROL 

0. 

000 

RXB 

0. 

000 

RYB 

0. 000 

RZB 

PL 

: PLV 

- 

3. 16 

PCH 

-10. 07 

YAW 

0. 

03 

ROL 

-0. 

011 

RXB 

0. 

068 

RYB 

0.014 

RZB 

PL 

: SM50 

-3 

30744 

X 

-281099 

Y - 

•19803 

4 Z 

- 

1033. 

87 

DX 

-759 

.41 

DY 

-401.51 

DZ 

PL 

: R3LV 

4 

20086 

X 

-1849 

Y 

226314 Z 


1452. 

11 
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Figure HI. Illustration of Graphical Plot for CPLV Plot Option 
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Figure H2. Illustration of Graphical Plot for RSBY Plot Option 



